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ABSTRACT 


A  practical  lumped  response  watershed  model  has  been 
developed  to  study  lake-watershed  systems  in  a  Canadian 
prairie  setting.  The  model  is  designed  to  utilize  basic 
climatic,  hydrological  and  physical  data,  such  as  is 
normally  available  for  lake  studies,  in  order  to  simulate 
variations  in  lake  water  quality  and  quantity. 

The  model  has  been  used  to  examine  the  role  that 
climatic,  hydrological  and  physical  variables  play  in 
controlling  the  hydrology  and  chemistry  of  prairie  lakes. 
These  analyses  indicated  that  climatic  parameters,  such  as 
precipitation,  evaporation  and  temperature  play  a  dominant 
role  in  influencing  lake  levels  and  salinity  of  shallow 
lakes  and  lakes  situated  in  small  watersheds.  Large  deep 
lakes  or  those  located  in  a  large  watershed  are  primarily 
controlled  by  groundwater  discharge  ana  surface  runoff. 
Physically  altering  the  watershed  may  also  produce 
significant  changes  to  the  lake  chemistry  and  surface 
elevations . 

The  practical  a pp iicab il ity  of  the  model  is 
demonstrated  by  applying  it  to  a  study  of  Baptiste  Lake, 
Alberta,  Canada,  where  a  favourable  compar ision  is  achieved 
between  the  real  and  simulated  h ydr ochemical  and  lake  level 
patterns.  The  model  was  primarily  used  to  illustrate  the 
importance  of  groundwater  outflow  to  Baptiste  Lake; 
indicating  that  groundwater  constitutes  a  significant 
proportion  of  all  recharge  to  the  lake. 
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I.  INTRODUCTION 

Statement  of  Tnes is 

Hydrological  studies  of  lake-wat ershed  systems  have 
concentrated  typically  on  quantifying  the  balance  between 
precipitation,  evaporation,  surface  runoff,  groundwater 
inflow  and  outflow  and  lake  discharge.  Although  the 
contribution  cr  components  other  than  groundwater  to  the 
water  balance  are  obtained  through  standard  measurement 
methods,  in  most  cases  groundwater  inflow  and  outflow  have 
been  unknown  quantities  and  their  significance  in  the  water 
budget  has  been  speculative. 

In  practice,  the  groundwater  contribution  to  lakes  is 
often  assumed  to  be  negligible  because  accurate  methods  to 
calculate  groundwater  recharge  and  discharge  are  not 
available.  However,  there  are  increasing  indications  that 
suggest  that  groundwater  contributions  to  stream  and  lake 
systems  may  be  considerable  (summarized  in  Schwartz,  1979). 
An  important  research  priority  is,  therefore,  the 
development  of  methods  that  consider  the  groundwater 
component  and  that  are  capable  of  determining  how 
significant  groundwater  is  in  the  hydrologic  balance  of  a 
lake- watershed. 

The  first  objective  of  this  thesis  is  to  develop  a 
practical  hydrologic  model  that  is  oriented  primarily  to 
lakes.  It  is  designed  to  utilize  basic  climatic  and 
hydrological  data  in  order  to  predict  variations  in  lake 
water  quality  and  quantity.  A  second  objective  is  to  use  the 
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model  to  examine  the  role  that  various  climatic  and 
watershed  variables,  such  as  precipitation,  evaporation, 
surface  runoff,  groundwater  and  seasonal  climate  changes, 
play  in  controlling  the  physical  and  chemical  hydrology  of 
lakes  in  a  prairie  setting.  The  final  objective  is  to 
demonstrate  the  practical  applicability  of  the  model  in 
evaluating  an  actual  lake-watershed  system.  The  model  will 
be  used  specifically  to  estimate  the  groundwater 
contribution  to  Baptiste  Lake,  Alberta,  and  to  determine  the 
impact  that  further  cottage  development  may  have  on  the 
quality  of  the  lake, 

Baptiste  Lake  has  neen  selected  for  study  because  there 
are  ongoing  government  programs  concerned  with  this  lake  and 
because  there  is  a  large  existing  data  base.  The  Water 
Quality  Control  Branch,  Pollution  Control  Division  of 
Alberta  Environment  are  establishing  background  limnological 
data  for  Baptiste  Lake,  are  defining  the  extent  and  sources 
of  pollution  in  the  lake  and  are  assessing  the  potential 
impact  of  a  proposed  cottage  development  (Water  Quality 
Control  Branch,  1977).  The  Earth  Sciences  Division  of 
Alberta  Environment  was  approached  by  the  Pollution  Control 
Division  to  assist  in  determining  the  guantity  and  quality 
of  groundwater  entering  the  lake.  The  study  of  Baptiste  Lake 
is  an  extension  of  the  field  work  and  the  preliminary 
analysis  of  those  data  carried  out  by  the  author  and 
incorporated  in  a  report  by  Trew  et  al.  (1979), 

Each  of  the  three  objectives  will  be  addressed  in 
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successive  chapters  of  the  thesis.  The  de velopaent  a 
theory  of  the  watershed  model  is  presented  in  Chapte 
Chapter  III  examines,  via  sensitivity  analyses,  the 
that  the  previously  mentioned  physical  parameters  ha 
lake  hydrology.  Chapter  IV  discusses  the  practical 
application  of  this  model  to  Baptiste  Lake. 
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The  initial  developments  in  hydroloqic  modelling  were 
water  budget  calculations  based  on  a  simple  accounting  of 
moisture  input,  outflow  and  storage.  These  calculations 
primarily  require  precipitation  and  temperature  as  input. 
Occasionaly  other  climatic  data  sue h  as  dew  point , 
temperature,  humidity,  solar  radiation,  wind  speed  and 
geographical  data  suen  as  soils  and  vegetations  are  used. 
Examples  of  the  application  of  these  methods  to  study 
water-balances  of  prairie  lakes  are  given  by  Laycock  {1968, 
1971,  1973),  Saskatchewan  Nelson  Basin  Board  (1972)  and 

Prairie  Province  Water  Board  {  1952)  .  The  drawback  of  these 
models  is  that  they  represent  only  a  crude  long  term 
accounting  of  water  in  the  hydroloqic  cycle.  Without  routing 
water  through  a  more  realistic  representation  of  a 
watershed,  the  temporal  variability  in  processes  such  as 
evaporation  and  groundwater  discharge,  cannot  accurately  be 
assessed . 

Commencing  in  the  early  1960‘s  modelling  progressed 
towards  simulating  the  hydrologic  cycle  by  routing  water 
through  storage  units  which  represent  components  of  the 
watershed.  The  watershed  is  simulated  either  as  a  single 
lumped  unit  or  as  a  set  of  lumped  models  which  has  the 
watershed  divided  into  smaller  sub-areas.  Bepresentation  of 
a  watershed  by  a  lumped  model  is  illustrated  on  Big  are  1  and 
Figure  2 . 

The  first  total  watershed  response  model  that  simulated 


EVAPORATION 


5 


■o 

0) 

-C 

in 

<1) 

4-J 

(u 

2 

fD 


C. 

c 


<u 


0) 

X 

CD 

CL 

E 

D 


a) 


D 

CD 


Impervious  Boundary 


6 


the  hydrological  cycle  was  the  Stanford  Watershed  Model 
(Crawford  and  Linsley,  1966),  While  this  model  was  designed 
for  streams,  the  requirements  of  a  hydrological  model  that 
they  advanced  can  readily  be  applied  to  a  lane-wa tershed : 

1.  "The  model  should  represent  the  hydrologic  regimes  of  a 
wide  variety  of  streams  and  rivers  with  a  high  order  of 
accuracy . 

2.  "It  should  be  easily  applied  to  different  watersheds 
with  existing  hydrologic  data. 

3.  "The  model  should  be  physically  relevant  so  that 
estimates  of  other  useful  data  in  addition  to 
stream! low,  such  as  overland  flow  or  actual 
evapo t ranspir at  ion ,  can  be  obtained." 

The  Stanford  Watershed  Model  simulates  portions  of  the 
hydrological  cycle  by  routing  and  accounting  for  all 
moisture  in  the  basin  over  a  period  of  time.  The  model 
consists  of  water  storage  units  which  represent  those 
components  of  the  real  system  such  as  the  soil  zone  or 
groundwater  zone,  which  act  to  store  water  for  various  time 
periods.  A  set  of  empirical  algebraic  equations  control  the 
movement  of  water  among  the  storage  units.  Parameters  within 
these  equations  are  used  to  characterize  features  of  the 
watershed  such  as  land  surface  and  soil  profile.  These 
parameters  are  calibrated  in  order  that  the  model 
predictions  of  streamfiow  will  simulate  recorded  streamflow. 
The  model  simulates  only  the  hydrologic  components  or  a 
watershed,  having  no  provisions  for  handling  water 
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chemistry.  The  model  has  undergone  many  revisions  and 
improvements  since  its  original  conception  and  these  are 
briefly  discussed  by  Linsley  et  ai.  (1975). 

Huff  (1968)  extended  the  Stamord  Watershed  Model  to 
simulate  the  hydrologic  transport  of  selected  radioactive 
aerosols.  This  was  later  reprogrammed  and  is  known  as  the 
Wisconsin  Hydrologic  Transport  Model  (Patterson  et  al. , 

19  74)  . 

Crawrord  and  Donigiac  (1973)  also  developed  a  transport 
model  based  on  the  Stanford  Watershed  Model.  The  model, 
known  as  the  Pesticide  Transport  and  Eunoff  Model,  simulates 
the  loss  of  pesticides  from  agricultural  lands  by 
"p iggy -backing”  the  pesticides  onto  the  movement  of  water 
and,  accounts  for  alteration,  volitization  and  degradation 
of  the  pesticides. 

Fields  and  Watson  (1975)  used  an  optimizing  procedure 
with  the  Wisconsin  Hydrologic  Transport  Model  to  determine 
input  paramaters  which  lead  to  the  best  fit  for  simulated 
and  observed  contamination  transport  and  hydrologic  values. 

Generally  there  is  little  or  no  indication  in  the 
literature  of  application  of  these  three  transport  models  to 
lake-watershed  systems.  In  each  case,  the  model  was  tested 
on  a  small  stream  watershed. 

Up  to  now,  hydrological  studies  of  lake  systems  have 
relied  on  more  conventional  water-balance  methods  (Laycock, 
1968,  1971,  1973;  Nursali  et  al. ,  1971).  Fritz  and  Krouse 

(1973),  using  dissolved  sodium  and  natural  stable  isotopes 
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of  oxygen  and  hydrogen  (oxygen-18  and  deuterium),  evaluated 
the  pattern  of  groundwater  inflow  to  a  prairie  lake  and 
characterized  the  evapcritic  regime.  Schwartz  and  Gallup 
(1978)  demonstrated  how  major  ions  (Ca,  Mg,  Na,  K,  HC03  ,  S04 
,  Cl),  can  be  utilized  to  determine  the  relative  importance 
of  groundwater  recharge  and  discharge  to  prairie  lakes. 
However,  these  studies  did  not  take  into  account  the 
time-varying  changes  that  occur  in  a  lake- watershed  system, 
such  as  the  flow  of  water  or  chemical  changes,  but  simply 
predicted  an  average  lake  chemistry  from  the  contributing 
components,  independent  of  the  time  at  which  the  data  was 
collec ted . 

A  logical  extension  of  the  chemical  balance  technique 
lies  in  the  routing  of  major  ions  through  a  prairie  lake 
system  over  a  period  of  time,  Sucn  a  model  could  be  a  useful 
tool  in  understanding  the  hydrology  of  a  prairie 
iak e- wa t ershed .  However,  as  the  above  discussions  explain,  a 
model  need  not  be  complex  in  order  to  adequately  represent  a 
lake- water shed  system.  Hence,  the  lumped  response  model 
developed  in  this  thesis  is  adequate  for  examining  features 
of  lake- watershed  systems  in  a  prairie  setting. 


. 


II.  THE  MODEL 


The  model  developed  in  this  thesis  is  a  lumped  response 
model,  similar  in  purpose  and  structure  to  the  Stanford 
Watershed  Model  and  the  Pesticide  Transport  Eunoff  Model.  It 
is  capable  of  simulating  the  hydrologic  responses  of  a 
lake-watershed  system  to  variations  in  precipitation, 
potential  evaporation  and  air  temperature.  The  uniqueness  of 
this  model  exists  in  the  fact  that  both  the  hydrologic  and 
mass  transport  functions  are  much  simpler  and  much  less 
input  data  are  required.  This  model  is  designed  to  operate 
with  only  a  basic  set  of  input  data.  In  this  respect,  the 
data  requirements  for  the  model  are  compatible  with  the  type 
of  information  commonly  available  for  lake  systems  in 
Alberta. 

The  model  operates  by  routing  water  throuqh  an 
idealized  lake-wa tershed  over  a  period  of  time.  At  each  time 
step,  water  flow  through  the  various  storage  units  and  the 
lake  is  determined.  Derived  values  such  as  lake  stage,  are 
calculated  from  given  sets  of  precipitation  and  potential 
evaporation  data.  Information  on  air  temperature  is  used  to 
simulate  hydrologic  conditions  during  winter,  such  as 
restricted  infiltration,  ice  formation  on  the  lake,  and  snow 
melt  in  the  spring. 

Mass,  in  the  form  of  ions,  is  also  cycled.  Each  major 
source  of  water  that  flows  to  the  lake  is  character ized  by  a 
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chemical  parameter  such  as  total  dissolved  solids  or  ionic 
composition.  These  sources  are  precipitation  failing 
directly  on  the  lake,  runoff  from  impervious  areas,  surface- 
runoff  or  overland  flow  and  groundwater  discharge.  The 
chemistry  of  the  lake  water  is  a  mixture  of  all  the  water 
sources  modified  by  processes  such  as  evaporation  or  ice 
formation. 

The  hydrologic  components  of  the  model  ana  tne  flow  of 
water  among  them  are  shown  in  detail  in  Figure  2. 

Precip itat io n  is  assumed  to  be  rain  if  the  temperature  is 
above  0°C,  or  snow  at,  or  oelow  0°C.  Snow  does  not  enter  the 
system  immediately  during  the  winter,  but  is  stored  until 
spring  when  it  is  allowed  to  melt.  The  difference  between 
precipitation  and  potential  evaporation  represents  water 
that  may  infiltrate  into  the  unsaturated  zone,  directly  to 
the  groundwater  zone,  cr  move  overland  to  the  lake  via 
rivers . 

Helted  snow  or  rain  which  falls  on  the  impervious  area 
is  added  directly  to  the  lake.  Impervious  areas  could,  for 
example,  be  comprised  cf  streams,  rivers,  lakes  or  urban 
areas  with  storm  sewers  and  roads,  where  precipitation  is 
transported  directly  tc  surface  waters,  fiunoff  from 
impervious  areas  must  be  discharged  directly  to  the  lake  or 
a  stream  before  the  land  can  fce  defined  as  an  impervious 
area.  As  is  usually  the  case,  all  runoff  from  some  ouildiags 
and  roads,  etc.  must  flow  over  pervious  areas  before 
reaching  surface  waters  and  hence,  in  this  thesis,  these  are 
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not  included  as  impervious  areas. 

The  unsaturated  zone  is  a  moisture  storage  unit 
overlying  the  water  table.  It  is  characterized  by  a  definite 
moisture  capacity  which  must  be  exceeded  before  recharge  can 
take  place  to  the  groundwater  zone.  Any  water  in  excess  of 
the  unsaturated  zone  moisture  capacity  is  availarle  for 
recharge  to  groundwater  storage  units.  However,  the  rate  of 
recharge  is  fixed.  Potential  recharge  in  excess  of  this 
amount  will  flow  to  surface  waters  as  overland  flow. 
Following  the  procedure  descriced  ry  Ayers  (1970),  surface 
runoff  is  defined  as  the  amount  of  water  availarle  for 
infiltration  that  is  in  excess  of  the  water  retained  m  the 
upper  and  lower  unsaturated  zones  and  the  maximum  allowacle 
recharge  to  the  groundwater  zone. 

The  groundwater  zone  is  composed  of  a  set  of 
groundwater  storage  elements  or  units  which  are  defined  with 
characteristics  resembling  actual  hydrostratigraphic  units. 
Waters  from  these  units  may  flow  directly  to  tne  surface 
water  system  or  to  other  groundwater  storage  elements. 
Groundwater  recharge  to,  and  discharge  from  the  storage 
elements  is  calculated  according  to  a  set  of  routing 
equations  defined  by  Dooge  (I960). 

Fivers  are  not  as  fully  constructed  by  tne  model  as 
they  would  be  in  the  real  hydrologic  environment.  Only  the 
groundwater  discharge  component  of  stream  flow  is  computed 
with  the  model.  Groundwater,  discharged  to  a  river,  is  a 
separate  parameter  because  one  of  the  objectives  of  the 
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model  is  to  determine  how  much  water  directly  enters  the 
lake  and  how  much  indirectly  enters  the  lake  via  rivers.  It 
is  assumed  that  all  surface  runoft  ana  impervious  flow  will 
reach  the  lake  whether  or  not  it  enters  a  river ,  and  hence, 
these  two  components  are  routed  directly  to  the  laxe.  Kiver 
flow  can  be  estimated  by  summing  the  calculations  of  surface 
runoff,  runoff  from  impervious  areas  and  groundwater 
discharged  to  rivers.  'Water  can  leave  lake  storage  as 
surface  outflow,  groundwater  outflow,  evaporation  or 
temporarily,  as  ice  in  the  winter. 


Theory  of  the  do  del 

The  Uosaturated  Zone 

The  unsaturated  zone  is  a  vital  component  of  a 
watershed  system  because  it  controls  the  amount  of  moisture 
that  will  recharge  the  groundwater  system.  The  unsaturated 
zone  also  acts  as  a  reservoir,  gaining  water  during  periods 
of  precipitation  and  being  depleted  during  periods  of 
evaporation. 

Water  is  routed  through  the  unsaturated  zone  to  the 
groundwater  units  where  the  water  table  is  not  near  or  at 
the  surface.  Otherwise  infiltration  is  assumed  to  enter 
groundwater  storage  directly.  The  unsaturated  zone  controis 
the  amount  of  recharge  to  the  groundwater  system  in  a  given 
period  of  time.  Therefore,  the  quantity  of  surface  runoff  is 
controlled  by  charac teristics  of  this  zone. 
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The  unsaturated  zone  (Fiqure  3),  is  divided  into  two 
pacts  (Ayers,  IS70);  an  upper  unsaturated  zone  which 
receives  direct  infiltration  and  underqoes  evaporation,  and 
a  lower  unsaturated  zone  which  receives  infiltration  from 
the  upper  unsaturated  zone  and  underqoes  evaporation  at  a 
reduced  rate.  Before  water  flows  no  the  lower  unsaturated 
zone,  the  moisture  storaqe  of  the  upper  zone  must  be 
satisfied.  Also,  before  evaporation  from  the  lower  zone 
occurs,  all  moisture  storaqe  in  the  upper  unsaturated  zone 
must  be  depleted. 

Linear  units  of  measurements  are  utilized  by  the  model 
for  precipit ation,  potential  evaporation,  recharge,  water 
storage  and  discharge.  The  units  are  determined  as  deptn 
over  the  area  of  the  watershed  (eg.  cm,  m) . 

Recharge  to  the  unsaturated  zone  is  determined  in  the 
following  way: 

.....(  2) 

R«=P-Ep 


where;  R  =  total  available  recharqe  to  the 

upper  unsaturated  zone 
P  =  precipitation 
Ep  -  potential  evaporation 
Water  storage  in  the  upper  unsaturated  zone  is 
determined  by  a  parameter  called  the  maximum  upper 
unsaturated  zone  storage  capacity.  Only  when  this  storage 
limit  is  exceeded  by  infiltration,  will  there  be  recharge  to 
the  lower  zone.  Any  infiltration,  or  infiltration  combined 
with  moisture  already  in  storage  that  is  in  excess  of  the 
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Figure  3-  Schematic  diagram  of  the  unsaturated  zone  element. 
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maximum  storaqe  capacity  is  routed  to  the  lower  unsaturated 
zone  according  to  the  tela tionship : 

K=  P-Ep-R',  . (  3) 

where:  Rk  =  total  available  recharge  to  the 

lower  unsaturatea  zone 

Ru2  =  actual  amount  of  water  recharged  to  the 
upper  zone  that  remains  there  in  storage 
Any  moisture  in  excess  of  the  lower  unsaturated  zone 
requirements  is  routed  to  the  groundwater  units  using  an 
expression  similar  to  (  2)  or  {  3); 

Rgw^P-Ep-  R'u2-R  u) 

where:  Rgw  =  total  available  recharge  to  the 

groundwater  zone  from  the  unsaturated 
zone 

R|z  =  actual  amount  of  water  recharged  to  the 
lower  zone  that  remains  there  in  storage 


or: 

Rgw'O  ,  if  (4)  <  0  0) 

From  {  2)  it  is  evident  that  evaporation  may  exceed 
infiltration.  In  this  case,  the  present  water  storage  in  the 
upper  unsaturated  zone  is  reduced  by  the  amount  of  excess 
evaporation.  If  the  evaporation  requirements  are  not 
satisfied  (i.e.  Rlz  <  0)  ,  then  evaporation,  defined  by  the 
following  equation,  taxes  place  at  a  reduced  rate  from  the 
lower  unsaturated  zone: 


■ 


where:  E  -  actual  amount  of  water  evaporated  from 

the  lower  unsaturated  zone 
S,z  -  maximum  lower  unsaturated  zone  moisture 
sto  ra  ge 

Sl2m  =  present  moisture  content  cf  lower 
un saturated  zone 

No  evaporation  is  allowed  from  the  groundwater  units  which 
are  overlain  ly  an  unsaturated  zone,  even  if  evaporation 
requirements  are  not  met  in  the  unsaturated  zone  (i.e.  Rgw 
cannot  be  less  than  zero) , 

After  the  moisture  requirements  of  the  unsaturated  zone 
have  been  satisfied  and  the  maximum  recharge  allocated  to 
the  groundwater  zone,  any  water  surplus  is  discharged  as 
surface  runoff.  This  calculation  is  made  using  the  following 
expression: 

Qs  -  Rgw  “  Rgwm  .....(7) 

where:  Qs  -  surface  runoff 

Rgwm  -  maximum  allowable  recharge  to  the 
groundwater  zone 

or: 

Qs  =  0  ,  if  (  7)  <  0  . (8) 

Runoff  from  impervious  areas  is  routed  directly  to  the 
lake.  The  amount  of  runoff  from  impervious  areas  is  defined 


by : 

Q,  --  P  •  r  .....(  9 ) 

where:  Q,  =  runoff  from  impervious  areas 
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y  =  proportion  of  basin  which  is  impervious 
If  such  areas  are  present  in  the  basin,  then  infiltration  to 
the  unsaturated  zone  is  reduced  by: 

I  =  P  -  Q,  . (  ,0) 

where:  I  =  infiltration  to  the  unsaturated  zone. 

The  term  I  replaces  P  in  equations  (  2)  ,  (3)  ana  (  4)  . 

The  model  was  designed  to  simulate  hydrologic 
conditions  during  the  winter,  as  well  as  during  the  summer. 
Before  the  unsaturated  zone  section  of  the  model  is  reached, 
a  check  is  made  in  the  program  to  determine  if  it  is  winter. 
During  winter,  there  will  be  no  mriltration  to  the 
unsaturated  zone  and  precipitation  is  stored  on  the  surface 
as  snow.  During  the  first  month  in  which  the  average  monthly 
temperature  is  above  0°C,  the  snow  melts  and  this  moisture 
is  added  to  precipitation  in  (10). 

Groundwater  Storage  Elements 
Dooge  (I960)  proposes  that  a  groundwater  system  can  be 
represented  by  one  or  more  linear  storage  units  or  elements. 
When  these  are  properl y  linked  together,  they  provide  an 
adequate  basis  for  predicting  groundwater  outtiow  when 
groundwater  recharge  is  known.  Practice  has  demonstrated 
that  this  simple  system  of  linear  storage  elements,  can 
successfully  model  a  real  groundwater  system  (Dooge,  I960; 

Ay  ers,  1970).  Input  or  recharge  to  an  element  can  occur  from 
three  possible  sources:  infiltrating  precipitation, 
groundwater  routed  from  another  element,  or  ftom  a  surface 
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source  such  as  a  stream  or  a  lake.  Groundwater  outflow  is 
routed  to  other  storage  elements,  to  a  surface  water  body, 
or  to  the  atmosphere  by  evaporation. 

Dooge  (I960)  descrines  four  types  of  groundwater 
storage  elements.  The  behavior  of  each  element  is 
charac terized  by  the  following  set  of  parameters;  storage 
delay  time,  recharge  tame  and  recharge  characteristics.  The 
constraints  placed  on  the  groundwater  storage  elements  and 
the  terminology  used  in  this  thesis  has  been  modified  and 
simplified  from  Dooge  {I960).  In  this  thesis  elements  are 
defined  as  fellows; 

1.  a  shallow  groundwater  element,  which  underlies  an 
unsaturated  zone, 

2.  a  shallow  groundwater  element  which  has  no  overlying 
unsaturated  zone, 

3.  a  deep  groundwater  element, 

4.  a  composite  groundwater  element. 

Shallow  groundwater  storage  elements  gain  water  from 
precipitation.  Type  I  elements  are  overlain  by  an 
unsaturated  zone  through  which  precipitation  infiltates  to 
the  groundwater  system  and  from  which  water  evaporates. 
Hence,  no  water  will  evaporate  from  the  Type  I  elements. 

Before  there  is  any  groundwater  recharge,  the  water 
storage  requirement  of  the  unsatuated  zone  must  first  be 
fulfilled.  In  a  case  where  the  water  table  is  located  at  or 
very  near  to  the  ground  surface,  precipitation  and 
evaporation  is  assumed  to  occur  directly  from  the  Type  2 
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groundwater  storage  element. 

Groundwater  discharge  fccm  the  shallow  storage  elements 
may  be  routed  either  to  other  elements  or  discharged 
directly  to  surface  water  bodies.  It  is  also  possible  for 
water  to  flow  from  a  surface  water  bod y  into  the  groundwater 
storage  element.  For  example,  during  a  period  of  sudden  rise 
in  the  level  of  a  river  or  a  lake  surface,  water  will  flow 
into  the  bank.  This  situation  is  assumed  to  apply  only  to 
the  two  shallow  groundwater  storage  elements. 

Fecharge  to  the  deep  groundwater  storage  element. 

Type  3,  can  only  come  from  other  elements  (either  shallow, 
deef  or  composite).  In  addition,  it  is  assumed  that 
evaporation  cannot  take  place  from  these  deeper  elements. 
Discharge  from  the  deep  groundwater  element  is  routed  to 
other  storage  elements  or  a  body  of  surface  water. 

‘The  Type  4  element  has  the  properties  of  both  the 
shallow  and  deep  groundwater  elements.  It  can  act  like  a 
shallow  element  with  discharge  by  outflow  and  evaporation. 
Once  tue  storage  is  depleted  to  a  minimum  limit,  the  element 
assumes  properties  of  a  deep  element,  discharging  by  outflow 
only.  The  physical  counterpart  of  this  type  of  storage 
element  is  a  hydrostra tigraphic  unit  where  groundwater 
undergoes  evaporation  until  the  water  table  is  lowered  to  a 
level  where  evaporation  no  longer  removes  water.  Further 
depletion  to  lower  the  water  table  will  only  occur  via 
groundwater  outflow. 

Figure  4  schematically  illustrates  a  typical 


‘ 

' 


2  I 

groundwater  storage  element.  Water  discharging  from  an 
element  is  allocated  to  various  recipients.  Ihe  figure  shows 
all  possible  sources  of  recharge  and  all  possible  recipients 
of  discharge.  Each  of  the  four  types  of  elements  will 
actually  use  only  part  of  these  components. 

As  previously  discussed,  only  shallow  groundwater 
storage  elements  can  be  recharged  by  water  leaving  the 
unsaturated  zone.  Ihe  quantity  of  recharge  defends  on  the 
proportion  of  the  basin  surface  area  occupied  by  the  storage 
element  as  follows: 

Rer,5-Rgw  if  an  unsaturated  zone  is  present  .....(II) 

or : 

Re=<5-I  with  no  unsaturated  zone  .....{1 2) 

where:  Re  =  proportion  of  recharge  allocated 
to  an  element 

<5  -  proportion  of  area  represented  by  a 

groundwater  storage  element 
I  =  percolation  to  the  groundwater  zone 
Groundwater  discharges  to  other  elements,  rivers  or 
directly  to  the  lake.  Ihe  outflow  is  proportionally  divided 


according  to: 

Q  g  w  “  Q  gw  * 

.....  (1  3) 

where:  Qgw 

=  proportion  of 

outflow 

from  an  element 

discharging  to 

a  given 

place 

Q'9w 

=  total  outflow 

from  an 

element 

* 


proportionality  factor 
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LEGEND: 

CC  ELEMENT  NUMBER 


C 

P-i  ’  * ‘  P  n 


© 

© 

© 
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GROUNDWATER  SALINITY  (mg/*) 
PROPORTION  OF  OUTFLOW  FROM 
LAKE 

ANOTHER  STORAGE  ELEMENT 
RIVER 

OUT  OF  BASIN 


ELEMENT 


Fi gu re  A 


Schematic  diagram  of  a  typical  groundwater  storage  element. 
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As  previously  noted,  evaporation  may  occur  from  a  Type 
2  shallow  groundwater  element.  The  amount  of  evaporation  is 
determined  by  a  negative  value  of  infiltration  as  determined 
by  equation  (10).  The  negative  value  of  infiltration  is  then 
inserted  into  equation  (12),  indicating  a  net  less  of  water. 

Groundwater  Routing  Equations 

Storage  in  an  element  is  not  cescnbed  by  a  simple 
balance  relationship  (for  example,  in  equation  {  1) ) ,  nut  is 

controlled  cy  a  parameter  termed  the  storage  delay  time. 
Hence,  storage  in  an  element  is  aeiined  by: 

S  =  K  •  q  .....(14) 

where:  S  =  storage  in  an  element 

q  =  outflow  from  an  element 
K=  storage  delay  time 

This  relationship  assumes  that  the  time  periods  in 
which  recharge  occurs  are  equal  and  that  the  rate  of 
recharge  in  the  time  period  is  constant. 

Dooge  {  I960)  developed  mathematical  expressions  to 
describe  the  discharge  of  groundwater  from  a  linear  storage- 
element  and  to  calculate  the  amount  of  storaqe  remaining  in 
an  element.  Once  the  volume  of  recharge  to  the  storage 
element  and  the  routing  coefficients  are  known,  the 
following  equations  hold: 


Qt  =  c0  Rt  +  c,Rm  +  c2Qt_,  - - (15) 

St  -  c3Rt  +  c4Qt 
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where:  Q,  =  volume  or  outflow  during  time  step 

Qt_,-  volume  of  outflow  during  the 
previous  time  step 

Rt  =  volume  of  recharge  over  time  step 
R,_,  =  volume  of  recharge  during  the 
previous  time  step 

Sf  -  total  storage  at  the  end  of  time 
cQ,  c, ,  c2 ,  c3 ,  c4  =  routing  coef  ficients 
Ihe  routing  coeffecients  are  dependent  only  on  values 
of  storage  delay  time  and  time  for  recharge  to  an  element. 
The  routing  coefficients  are  defined  as: 


c 


o 


K 


C|=  T 


.  -H\  -  .-hr 


Co  =  e 


zH 


c 


3 


C4  = 


1 


-1 


(17) 


(Id) 


(2d) 


(21) 


where:  K  =  storage  delay  time  of  an  element 
T  =  duration  of  uniform  recharge  to 
an  element 

Discharge  from  an  element  declines  exponentially  and 
recharge  increases  exponentially  (figure  5).  The  reader  is 


RELATIVE  STORAGE  ^ ,  RELATIVE  STORAGE 


25 


RELATIVE  TIME 


Figure  5-  Groundwater  outflow  from,  and  storage  in  a 
groundwater  element. 


referred  to  Dooje  (  196  0)  for  a  complete  mathematical 
development  cf  these  routing  coefficients.  Also,  because  of 
the  time  delay  factor,  recharge  to  an  element  may  be 
occurring  but  the  discharge  continues  to  be  dominant, 
indicating  a  net  less  of  storage. 

Should  a  shallow  element  that  is  adjacent  to  a  body  of 
surface  water  experience  excessive  evaporation,  the  storage 
in  the  element  will  drep  to  zero,  after  which  further  water 
requirements  for  evaporation  will  be  satisfied  by  blow  from 
surface  waters  to  the  storage  element.  During  this  time  only 
enough  water  to  fulfill  the  evaporation  requirements  will  be 
obtained  from  the  surface  water  body  and  hence,  storage  in 
the  element  remains  at  zero. 

If  recharge  from  a  body  of  surface  water  is  not  allowed 
but  evaporation  takes  place,  groundwater  discharge  will 
continue  until  storage  is  zero,  at  which  time  both  the 
storage  and  outflow  will  remain  at  zero  until  the  storage 
element  is  recharged.  However,  zero  storage  is  usually 
attained  during  the  time  interval  rather  than  at  the  end  of 
the  time  interval.  The  following  equation  (Dooge,  I960) 
provides  an  adjustment  to  the  discharge  calculated  in 
equation  (15)  to  allow  for  outflow  during  only  part  of  the 
time  period: 


<Vc5S 


t-l 


(22) 


where:  St_,  =  storage  in  element  during  the 

previous  time  interval 
Qt  =  adjusted  outflow  from  the  storage 


' 

element  for  fart  of  time  step 
when  evaporation  is  occurring 
routing  coefficient  defined  as; 


c 
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C5=  1  -  lQ9e(l~A) 

A 


where ; 


123) 


(24) 


The  La he  System 

The  climate  of  the  Canadian  prairies  is  characterized 
by  a  large  vanaticn  in  temperature  and  the  nature  of 
precipitation,  which  affects  the  hydrology  of  a  lahe.  It  is 
important  for  a  hydrologic  model  to  account  for  these 
climatic  conditions  in  order  that  water  is  routed  through 
the  hydrologic  elements  properly. 

In  this  model,  a  year  is  divided  into  three  climatic 
divisions.  Defined  primarily  on  the  basis  of  average  monthly 
temperature,  they  are: 

winter  months;  average  monthly  temperatures  ate  less 

than  0  °  C . 

spring  month:  first  month  following  winter  which  has  a 

temperature  greater  than  0°C. 

summer  months:  months  between  spring  with  temperatures 

above  0°C. 

During  summer  months,  the  sources  of  water  for  the  lake 
are  surface  runoff,  runoff  from  impervious  areas,  river 
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flow,  groundwater  discharge  and  precipitation  tailing  on  the 
lake,  water  is  lost  from  the  lake  Ly  evaporation,  surface 
outflow  and  discharge  via  ground wa ter  flow. 

In  winter,  only  groundwater  and  river  flow  are  assumed 
to  add  water  to  the  lake.  Snow  accumulates  and  is  not  free 
to  infiltrate  or  runoff  over  the  land  surface.  It  is  assumed 
that,  during  the  wintei,  water  does  not  evaporate  from  the 
lake,  the  ice  or  the  snow  pack  ana  that  all  snow  falling  on 
the  the  ice  remains  until  spring.  Water  temperature  is 
assumed  to  decline  to  G°C  during  tne  initial  half  of  the 
first  winter  month  because  air  temperatures  are  usually 
colder  than  water  temperature  during  the  fall.  During  this 
time,  evaporation  takes  place  from  the  lake  tut  there  is  no 
surface  runoff  or  infiltration. 

All  components  of  recharge  and  discharge  that  exist  in 
a  summer  month  are  assumed  to  te  active  in  the  spring.  At 
this  time,  the  snow  on  the  watershed,  and  ice  and  snow  on 
the  lake  are  assumed  to  melt.  Two  assumptions  are  made  tor 
the  spring  month:  (1)  all  ice  and  snow  has  melted  by  the  end 
of  this  month  and,  (2)  the  ground  thaws  rapidly  during  the 
spring,  allowing  water  to  infiltrate  to  the  unsaturated  and 
groundwater  zones.  These  assumptions  are  stated  primarily  to 
simplify  the  model  and  are  justified  because  the  large  time 
step  that  is  used  by  the  model  would  cover  the  actual  date 
of  snow  melt  and  thawinq  of  the  ground. 

Two  pathways  of  outflow  from  the  lake  are  via  the 
surface  and  subsurface.  Surface  water  outflow  is  calculated 
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in  the  model  from  the  volume  of  lake  water  above  the  surface 
water  outlet  elevation.  A  capaci ty-discharqe  curve  function 
for  the  lake  should  normally  re  provided  by  the  model  user 
(discussed  in  the  following  Operation  or  tfodel  section). 
Outflow  via  groundwater  is  simply  removed  from  the 
calculations  rather  than  ty  movinq  through  groundwater 
storage  elements.  The  quantity  of  groundwater  is  calculated 
by  a  defined  depth  of  water  loss  from  the  lake  times  the 
area  of  the  lake. 

Ice  Formation 

The  formation  of  ice  is  an  important  process  in  lakes 
of  the  Canadian  prairies.  Its  importance  is  evident  when  one 
considers  that  lakes  are  often  covered  with  ice  from 
November  to  April,  and  that  a  large,  proportion  of  the  total 
water  volume  of  the  lake  may  be  frozen.  One  major  result:  of 
ice  formation  is  partitioning  of  ions  between  the  unfrozen 
lake  water  and  ice.  The  net  effect  is  for  the  ions  to  become 
concentrated  in  the  unfrozen  water.  Therefore,  an  ice 
forming  component  in  the  model  is  essential  for  a  rigorous 
description  of  the  behaviour  ct  lakes  in  a  prairie  setting 
during  winter. 

A  simple  expression  describing  the  formation  of  ice, 
assuming  that  the  temperature  gradient  is  uniform  through 
the  ice,  is  the  Stefan  Equation  (Fertuck  et  ai. ,  1970; 

Ingersoll  et  al. ,  1948;  and  Pounder,  1965): 


. 
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D2T^/(W')dt  . <25) 

*0 

where:  D  ■=  ice  thickness 

k,  =  thermal  conductivity  of  ice 
Pi  =  ice  density 
Li  =  heat  of  fusion  of  ice 
Tf  -  temperature  at  freezing  point 
T.  =  temperature  of  ice  surface 
t0  -  time  at  which  ice  formation  commences 
t  =  time  at  which  ice  thickness  is  calculated 
Equation  (25)  cannot  me  directly  applied  to  formation 
of  ice  on  lakes  because  most  often  air  temperature  is  only 
available,  while  the  temperature  of  the  surface  of  the  ice 
is  actually  required.  Direct  substitution  of  air 
temperatures  m  equation  (25)  will  not  produce  accurate 
results  because  snow  on  the  ice  is  an  insulator,  resulting 
in  a  considerable  difference  between  air  and  ice  surrace 
tern  {.e  ra  t  ures. 

Fertuck  et  al.  (  19  7  0)  discussed  a  one-dimensional 
finite  difference  model  for  the  freezing  of  ice  under  a 
cover  of  snow  given  only  air  temperature,  wind  speed  and  the 
depth  of  snow  on  the  ice.  In  their  representation  of  the 
Stefan  Equation,  the  temperature  at  any  point  in  the  ice  or 
snow  is  calculated  by  (  p,  ,  H  ,  k  ,  are  assumed  to  be 


constan  t )  : 


3  1 


"n.t  +1 


Tn-U  +  (M  -  2)  -  Tn,t  +  Tn  +  l,t 

M 


(26) 


where:  Tnt  =  temperature  at  a  finite  difference  node 
n  =  subscript  used  to  designate  the  finite 
difference  node 


M  - 


p  H  ax' 
k  At 


(27) 


H  =  specific  heat  of  ice  or  snow 
k  =  thermal  conductivity  of  ice  or  snow 
Ax  =  vertical  distance  between  nodes 


At  — 


time  step 


The  temperature  of  the  snow  at  the  snow-ice  interface  i: 
aetermined  by  choosing  Ax  of  the  snow  so  that: 


AX, 


1  n  -  1  ,  t 


T„.t 


AX; 


ki.f  Tn  +  l,t 


(28) 


where:  ks  -  thermal  conductivity  of  snow 

A xs  =  vertical  spacing  of  nodes  in  the  snow 
AX;  =  vertical  spacing  of  nodes  in  the  ice 
and  the  temperature  at  the  snow-air  interface  (set  At  in 
eguation  (27)  so  that  M  ■=  2)  is: 


t,  +  N-T„, 

1  t 


1  +  N 


(29) 
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where:  Ta  t  =  dir  temperature 

ki  -  hAX  . (33) 

k 

s 

h  =  heat  transfer  coeincieut 
Ih is  f or muia tion  has  been  simplified  because  the 
watershed  mcdcl  requires  cnly  an  approximation  solution  to 
ice  growth.  Cnly  3  nodal  points  were  defined  (as  compared  to 
6  used  by  Fertuck  et  al.  ),  with  n- I  at  the  snow-air 
interface,  n  =  2  at  the  snow-ice  interface  and  n — 3  at 
ice-water  interface.  The  influence  of  wind  is  also 
neglected.  With  these  simplif ications,  equation  {23) 


becomes: 

T,  t  =  Ta.t  •  • 

and  the  temperature  of  the  ether  two  nodes  are: 

f  ,  =  ft 

T3tt  Tf.t 

where:  Tiit  =  temperature  at  the  ice-snow  interface 
Tf  t  =  temperature  of  the  freezing  point 
of  water  (0°C) 

Therefore,  equation  (28)  can  fe  rewritten  as: 


(3!) 


(32) 

(33) 


A  Xs 


_K 

AX; 


A_ 

AX, 


(34) 


Equation  (34)  is  substituted  in  equation  (26)  to  give  the 
ice  formation  equation: 
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/-2Ti  t  k,  t  1 

D  =  \  LiPi  . (35) 

Ice  begins  forming  during  the  first  month  in  which  the 
average  nonthly  temperature  is  below  0°C,  and  continues  as 
long  as  the  temperature  remains  below  0°C.  Ice  thicknesses 
are  calculated  during  this  month  according  to  equation  (35), 
assuming  there  is  no  snow  on  the  ice;Tit  =  TQt.  Ihe 
calculations  for  the  following  winter  months  use  equation 
(35)  with  Ti  t  calculated  from  equation  (34). 

Snow  is  allowed  to  accumulate  on  the  ice.  The  depth  of 
snow  is  taken  to  be  10  times  the  liquid  equivalent  of  the 
precipitation.  Normally,  as  time  progresses  snow  will 
compact  and  become  denser.  However,  in  this  model  its 
density  remains  equivalent  to  that  of  fresh  snow.  It  is 
assumed  that  a  thick  layer  of  snow  with  a  low  density  has 
the  same  insulating  effect  as  a  thin  layer  of  dense  snow. 
This  assumption  is  based  on  the  fact  that  the  neat 
conductance  is  a  function  of  the  depth  of  the  scow  times  its 
density.  As  the  snow  compacts,  its  depth  decreases  and  its 
de  ns  lty  increases.  While  the  relationship  is  not  linear,  it 
is  adequate  for  a  simple  ice  forming  model. 

Variables  used  in  the  ice  forming  equations  were 
obtained  from  Fertuck  et  al.  (IS7u)  because  their  study 
area;  Saskatoon,  Saskatchewan  and  this  one,  have  similar 
climatic  and  geographical  settings.  Pertinent  values  are 
listed  on  Table  1.  Substituting  the  values  shewn  on  Table  1 


. 


■ 


into  equation  (35)  will  simplify  tne  ice  terming  equation 
to : 


D, 


V 


2  25  •  L 


(36) 


(37) 


where:  D,  -  ice  formation  during  the  first 
winter  month 

Dt  -  ice  formation  during  the  following 
winter  months 


fable  i. 

VALUL3  OF  TIE  ICE  FOBBING  MC  DLL  PAR  AM  Ei  hnC 


ks  =  1  .  729*  10-*  cal/sec-cm-°c 
ki  =  5 .  3 5  0 ^  1  0~  3  cal/sec-cm-°C 
ps=  0.  130  gra/cm3 

Li  Pi=  6.276*  10-=  cm  2/sec-°C 

(ice  with  a  snow  cover) 

_k 

LjPi=  4.88!*I3~S  cm2/sec-°C 

(ice  without  a  snow  cover) 

All  values  are  from  Fertuck  et  al.  (1970) 
except  ps  which  is  from  Pounder  (1965). 
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The  ice  melts  daring  the  first  month  which  has  an 
average  monthly  temperature  alove  0°C.  All  ice  and  snow  is 
melted  by  the  end  of  this  month.  This  assumption  has  been 
substantiated  by  Pounder  (1965)  and  by  observation  at 
Baptiste  Lake.  Snow  is  converted  to  water  at  a  rate  of  10  cm 
of  snow  being  equivalent  to  1  cm  of  water  and  ice  at  a  ratio 
of  1:1. 

Chemical  Transport  and  Balance 
The  transport  of  dissolved  constituents  through  the 
watershed  is  based  on  a  simple  routing  and  accounting 
procedure.  Although  the  model  will  now  only  consider  one 
chemical  parameter,  it  could  be  extended  to  handle  several. 
Chemical  reactions,  dissolution,  precipitation  or 
stratifying  the  lake  are  net  included  in  the  model  directly. 
However,  because  the  chemistry  of  water  rrom  each  watershed 
component  can  be  specified,  all  important  chemical  processes 
operating  in  the  basin  are  directly  considered.  The  chetricai 
quality  of  the  lake  water  is  calculated  according  to  the 
dilution  equation: 


C  = 


C,Q|+  C;Q2  +  C  Q 
Q,  +  Q2+ . . . .  +  Qi 


(38) 


where:  C  =  concentration  of  the  lake 

C,  =  concentration  of  the  water  sources 
Qj  =  volume  of  flow  of  the  water  sources 
The  chenical  inputs  to  the  laxe  originate  from 
precipitation,  impervious  flow,  surface  runoff  and 


’ 
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groundva ter.  Calculations  of  the  cont ributions  from  the 
first  three  are  straight  forward.  The  total  inflow  is  the 
chemical  quality  of  these  sources  rimes  their  respective 
volumes  of  water.  However,  groundwater  discharge  may 
originate  from  several  storage  elements,  each  with  a 
different  water  quality.  Groundwater  can  discharge  directly 
to  the  lake  or  to  rivers.  Rivers  also  add  chemical 
components  to  trie  lake.  However,  as  previously  discussed 
with  respect  to  the  discharge  cf  water  to  rivers,  tne 
chemistry  of  the  rivers  are  not  defined  but  are  the  same  as 
the  groundwater  that  discharges  to  them  cecause  no  surface 
runoff  or  precipitation  is  mi  ied  with  tnis  water. 

The  model  can  account  for  the  concentration  or  ions  in 
the  lake  water  by  evaporation.  It  is  assumed  that  no 
dissolved  constituents  are  lost  through  the  evaporation 
process  or  no  salts  ate  precipitated  as  the  salinity  of  the 
lake  increases.  This  assumption  rs  valid  only  if  the  volume 
of  water  ir  the  lake  rs  sufficiently  large  so  that 
evaporation  removes  only  a  small  portion  of  lake  water. 

During  the  winter,  ice  will  form  on  the  surface  of  the 
lake  from  lake  water.  As  the  ice  forms,  ions  are  generally 
excluded  from  the  ice  but  remain  in  the  lake  water,  thus 
increasing  the  salinity  of  the  lake.  Snow,  accumulating  on 
the  ice,  will  have  the  same  chemistry  as  precrpitation ,  rut 
this  is  not  added  to  the  lake  until  the  ice  and  snow  melt 
during  the  spring. 

Ions  are  removed  from  the  lake  only  through  natural 


flushing  due  to  the  outflow  of  surface  water  and  discharge 
via  groundwater. 

Operation  of  toe  Com pu  te  r  Program 

The  model  was  programed  in  Fortran  IV  and  run  on  the 
AMDAHL  470V/6  computer  at  the  University  of  Alberta.  The 
program  consists  of  a  main  program  and  four  subroutines  that 
must  be  supplied  by  the  model  user.  A  complete  listing  of 
the  program  and  user’s  manual  are  is  provided  rn  Appendices. 

The  basic  sequence  of  operations  involved  in  the  model 
is  shown  in  figure  6.  The  program  is  designed  to 
continuously  loop  through  the  simulated  watershed  processes 
using  a  time  steps  of  one  month.  5. I,  units  are  used 
throughout  the  program. 

Some  of  the  input  parameters  may  not  te  readily 
attainable  from  actual  field  data,  but  require  a  trial  and 
error  process  to  determine  a  set  of  values  for  these 
parameters  which  will  produce  an  acceptable  fit  between 
recorded  and  simulated  watershed  responses.  A  detailed 
description  of  the  calibration  procedure  tor  the  model 
follows. 

It  should  be  noted  that  this  model  produces  an 
approximation  to  actual  field  setting,  and  an  exact 
correspondence  between  measured  and  predicted  responses 
should  net  be  expected.  Simulation  and  calibration  runs 
should  continue  only  until  the  user  is  satisfied  that  a 
physically  reasonable  set  of  model  parameters  yield 
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6.  General  flow  chart  of  the  lake-watershed  model. 


Figure 
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predictions  which  are  adequate  for  the  particular  reeds  of 
the  study. 

User  Supplied  Functions 

Four  functions,  which  determine  the  physical  hydroloqic 
characteristics  of  the  lake.  Bust  be  supplied  by  the  model 
user.  Ihe  data  required  to  determine  the  functions  should  be 
obtained  from  hydrographic  surveys  of  the  lake  and  rivers. 
Should  the  data  not  exist,  reasonable  estimates  can  Ce  made 
from  published  information,  The  functions  are: 

AEEALK  :  Function  used  to  calculate  the  area  of  the  laxe 
in  km2,  given  a  lake  stage  in  &. 

ELEVLK  :  Function  which  finds  the  lake  stage,  in  meters 
above  datum,  given  the  volume  of  the  lake  in  m3. 

VOLLK  :  Function  required  to  calculate  the  volume  of 

water  m  the  lake,  in  m3,  given  the  laxe  stage  in  m. 

OIJTFLO  :  Function  which  gives  the  surface  water  outflow 
from  the  lake,  m  m3,  given  the  lake  stage. 

Input  parameters  and  the  calling  sequence  for  the 
functions  are  supplied  in  Appendix  I.  Examples  of  the 
functions  are  provided  at  the  end  the  program  listing. 

Input 

Two  types  of  input  are  required  for  the  program:  firsx, 
watershed  parameters  tnat  define  tne  basin  and  second, 
climatic  variables.  The  watershed  parameters  consist  of 
information  required  to  define  the  physical  characteristics 
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of  the  lake  and  watershed,  the  groundwater  storage  elements, 
chemical  composition  of  the  various  waters  ana  the  route 
taken  by  the  groundwater  to  reach  the  lake.  The  actual 
parameters  arising  from  this  information  are  defined  only 
once  in  the  model  and  do  not  change  with  tame.  Watershed 
parameters  can  be  sub-divided  into  two  groups  on  the  basis 
of  information  which  can  be  attained  from,  field  measurements 
and  those  which  cannot.  Lake- watershed  parameters  not 
obtainable  by  field  measurements  require  calibration  t raa Is 
to  select  a  value  which  appear  to  give  the  best  fit  with 
actual  field  records.  The  problems  associated  with  choosing 
the  appropriate  values  for  these  parameters  is  considered  in 
a  following  section  called  Calibration. 

The  second  group  of  input  comprise  the  climatic 
variables.  These  are  input  at  each  successive  time  steps. 

The  climatic  variables  are,  precipitation,  potential 
evaporation  and  air  temperature. 

A  complete  description  of  procedures  tor  preparing  the 
input  data  and  for  dimensioning  the  arrays  is  presented  in 
Appendix  I.  A  sample  input  is  shown  in  Appendix  III. 

A  few  points  should  be  noted  regarding  particular  input 
values.  Each  groundwater  element  is  identified  by  number. 

The  element  numbering  should  start  at  I  and  increase 
seguencially  by  1  generally  following  tne  routing  path  taken 
by  the  groundwater  as  it  flows  through  the  model.  A  new  set 
of  element  numbers,  starting  at  1,  should  be  provided  for 
each  subarea.  Figure  7  shows  some  examples  of  element 
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Figure  7.  Examples  of  numbering  storage  elements. 
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numbering . 

There  are  tour  types  of  groundwater  storage  elements, 
each  having  particular  recharge  and  discharge  properties. 

The  type  of  element  is  identified,  by  the  model,  solely  on 
the  basis  of  input  parameters  of  recharge  to  the  element. 
Tests  to  determine  if  recharge  to  an  element  may  originate 
from  the  unsaturated  zone  or  precipitation  (KEPREC) ,  if 
evaporation  is  allowed  (NEGBEC),  and  if  recharge  from  a 
surface  water  body  occurs  from  an  element  (NEGQGfc),  are 
expressed  as  *T*  if  it  is  allowed  or  1 E*  if  it  is  nor. 

Hence,  a  shallow  storage  element  of  Type  I  has  h EPF  EC  -  T, 
NEGBEC  =  E  and  NEGQGW  =  T  or  f.  Type  2  elements  have  RIPE EC 
=  T,  NEGR EC  -  T  and  NEGQGW  =  1  or  F.  beep  storage  elements 
are  identified  by  REPREC  -  f,  NEGEEC  =  F  and  NEGQGW  =  F. 

The  proportion  of  the  surface  area  of  the  watershed 
occupied  by  a  shallow  groundwater  storage  element  and 
(PCTAE)  the  proportion  of  groundwater  outflow  from  an 
element  allocated  to  other  various  storage  units  (PC'IELM) 
are  expressed  as  decimal  fractions.  1CELM  designates  where 
the  discharge  from  an  element  is  routed.  Values  from  1  to  89 
are  reserved  for  other  storage  elements,  90  to  98  indicate 
discharge  to  rivers  and  99  defines  output  to  the  lake.  The 
first  value  of  PCTELM  and  TCELh  refers  to  discharge  to  the 
lake  and  must  always  be  specified  even  if  no  discharge 
occurs  from  the  storage  element  to  the  lake.  In  this  case, 
the  first  value  of  PCTELM  is  0.  ana  the  first  value  of  TOELfl 


is  99 
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All  parameters  and  variables  used  by  this  model  require 
S.I.  units.  It  is  important  to  keep  the  units  of  the  input 
parameters  consistent  with  the  definitions  given  in 
Appendix  I. 

Output 

Output  consists  of  calculated  hydrologic  and 
hydrochemical  values  and  input  parameters  and  vanacles.  All 
input  parameters  are  printed,  providing  the  useL  with  a 
check.  The  following  values  are  printed  at  the  end  of  each 
time  interval: 

1.  the  record  counter, 

2.  a  subarea  rdentifier, 

3.  air  temperature,  in  °C, 

4.  precipitation,  in  cm, 

5.  potential  evaporation,  in  cm, 

6.  calculated  actual  evaporation,  in  cm  (does  not 
include  actual  evaporation  irons  the  groundwater 
zone) , 

7.  total  flow  from  the  impervious  portion  of  each 
subbasin,  in  m3, 

3.  moisture  content  of  the  upper  unsaturated  zone,  in 
cm , 

9.  moisture  content  of  the  lower  unsaturated  zone,  in 
cm, 

10.  recharge  to  the  groundwater  zone  rrom  the 
unsaturated  zone,  in  cm. 
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11.  infiltration  to  the  ursaturated  zone  or  groundwater 
zone  with  nooverlying  unsaturated  zone,  in  cm, 

12.  groundwater  discharge  to  the  lake,  in  cm, 

13.  impervious  area  flow  to  the  lake,  in  cm, 

14.  groundwater  discharge  to  the  rivers,  in  cm, 

15.  quantity  of  surface  runoff  to  the  lake,  in  cm, 

16.  total  outflow  from  the  lake  via  groundwater,  ns3, 

17.  the  lake  stage,  in  m  above  a  datum, 

18.  total  recharge  to  the  lake,  m  dp  3 , 

19.  total  surface  water  outflow  ,  in  m 3 , 

20.  chemical  quality  of  the  lake  water,  in  mg/1, 

2  1.  ice  thickness,  m  cm. 

An  Example  of  output  is  p re sen  ted  in  Appendix  IV. 

Calibration  of  the  Model 

Parameters  that  are  not  easily  measured  in  the  field, 
can  only  be  determined  by  a  calibration  procedure  in  which 
the  data  ace  varied  systema tically  from  trial  to  trial  to 
derive  a  set  of  optimum  values  which  will  best  simulate  the 
actual  watershed  records.  In  seine  models  this  parameter 
selection  procedure  is  carried  out  as  part  of  the  algorithm 
(Fields  and  fcatson,  1975).  However,  in  this  worx,  a  trial 
and  error  slelection  process  has  been  utilized.  The 
parameters  calibrated  are: 

PCTIMP:  the  percentage  of  the  basin  which  is  impervious 
QGWP  :  previous  groundwater  outflow  from  an  element 
HELM?  :  previous  recharge  to  an  groundwater  element 
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maximum  moisture  content  of  the  upper 
unsa turn ted  zone 

LSZMAX:  maximum  moisture  content  of  the  lower 
unsaturated  zone 

U5Z  Z  P  :  initial  moisture  content  of  the  upper 
unsaturated  zone 

L 5 Z Z P  :  initial  moisture  content  of  the  lower 
un saturated  zone 


DT  :  storage  delay  time  from  a  groundwater  e  feme  lit 
SGWM  :  maximum  ailowaole  recharge  to  the  groundwater 
zone  from  the  unsaturated  zone 
GwOLK  :  lake  discharge  via  groundwater 
PCTELM :  proportion  of  outflow  from  an  element 
fohile  it  appears  as  if  there  are  many  parameters  that  trust 
he  calibrated,  the  problem  is  simplified  because  most  values 
lie  within  a  narrow  range. 

PCf IMP  can  be  establ isned  from  a  topographic  map  hy 
measuring  the  area  of  surface  water  (excluding  the  lake)  and 
uroan  areas  with  runoff  flowing  directly  to  tne  lake. 

Because  PC1ELM  has  the  effect  of  increasing  surface  runoff 
to  the  lake  and  hence  increasing  the  level  of  the  lake,  some 
minor  adjustments  may  he  required.  Generally,  the  effect  of 
increasing  PCTIM  P  will  be  to  decrease  groundwater  recharge 
and  discharge  to  the  lake  and  increase  runoff  from 
impervious  areas. 

QGWP  and  RSLMP  are  most  easily  estimated  by  using  as  a 
starting  point,  a  month  following  a  dry  period  or  the  winter 
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where  evaporation  or  groundwater  outrlow  have  almost 
depleted  the  groundwater  storage  in  the  elements.  At.  this 
time  QGWP  and  HELM?  are  assumed  to  he  zero.  After  results 
from  several  years  of  records  are  available,  a  fetter 
estimate  may  he  obtained  by  examining  groundwater  outflow 
and  recharge  during  the  final  simulation  period  for  these 
moii  ths. 

USZM  and  LSZM  control  the  amount  of  moisture  retained 
in  the  unsaturated  zone,  which  indirectly  limits  the  amount 
of  water  that  is  recharging  the  lake  in  the  torn  of 
groundwater  discharge  or  surface  runoff.  These  are 
especially  influential  during  the  spring  when  the  snow  pack 
is  melted  and  runoff  is  increased.  Throughout  the  summer  in 
the  prairies,  it  has  been  shown  generally  than  the  moisture 
capacity  of  this  zone  should  not  allow  recharge  to  the 
groundwater  system  unless  an  unexp ectant ly  long-term  or 
intense  rainfall  occurs.  USZM  and  LSZfl,  along  with  OSZMP  and 
LSZMP ,  can  only  be  determined  through  a  trial  and  error 
procedure,  out  initial  values  should  be  between  10  and  15  cm 
for  LSZM  +  OSZH  in  Alberta  (Laycock,  1971).  In  general,  an 
increasingly  larger  total  unsaturateo  zone  requires  more 
water  to  satisfy  this  zone's  moisture  capacity  before 
moisture  can  percolate  to  the  groundwater  zone,  thus 
delaying  the  time  for  recharge  tc  the  groundwater  zone. 

DT,  the  storage  delay  time,  and  IT,  the  time  tor 
recharge  to  a  groundwater  element,  influence  the  time  and 
volume  of  groundwater  that  is  discharged  to  the  lake.  In 
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this  model,  the  time  step  is  cce  month,  so  BT  is  always  I. 
While  the  value  of  DT  is  usually  between  1  and  5,  it  still 
must  be  established  by  trial  and  error. 

The  storage  delay  time  and  time  for  recharge  to  an 
element  were  plotted  to  illustrate  the  influence  that  they 
have  cn  recharge  to,  storage  in,  and  discnarge  from  a 
groundwater  storage  element,  figure  8  shows  a  plot  of  DT:ET 


Figure  8.  Analysis  of  groundwater  recharge 
times  and  discharge  times. 

versus  the  percentage  of  the  recharge  to  an  element  tnat  is 
discharged  and  the  volume  of  water  remaining  in  storage.  As 
DT:KT  increases,  the  rate  and  total  volume  of  recharge  to 
the  lake  and  the  amount  of  groundwater  discharge  to  the  lake 
will  decrease. 

Three  further  parameters,  HGWM,  GwOLK  and  PCIEL^ 
require  calibrating  to  determine  optimum  values  for 
simulating  the  watershed  responses.  The  maximum  allowable 
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recharge  to  the  groundwater,  RGViM,  should  be  calibrated  by 
trial  and  error.  The  value  of  JKGWfi  controls  groundwater 
discharge  and  surface  runoff.  It  is  expected  to  lie  between 
0  and  5  cm.  GWOLK  controls  the  lake  level  in  the  calculation 
of  water  removed  from  the  lake  during  each  time  step.  The 
value  should  be  estimated  by  trial  and  error.  The  proportion 
of  outflow  from  an  element  that  is  discharged  (PCTLLM) ,  to 
other  elements,  rivers  or  the  lake,  should  first  be 
estimated  from  the  geologic  data.  Finial  adjustments  to 
PC T EL M  require  calibration. 

Calibration  should  initially  only  involve  the  water 
routing  components  of  the  model  to  approximate  lake 
elevation,  flow  via  rivers  and  discharge  or  lake  water  by  a 
surface  outlet.  Only  small  adjustments  should  he  required  to 
these  parameters  during  the  next  tew  runs  to  match  the  water 
chemistry  of  the  lake  and  rivers,  because  the  chemical  data 
from  the  watershed  can  be  reasonably  well  determined,  the 
extent  of  fit  afforded  the  first  few  simulations  may  provide 
a  further  independent  check  on  the  appropriateness  of  the 
hydrologic  parameters. 

It  is  best  to  work  with  long  term  simulations  for  two 
reasons.  A  longer  run  provides  more  varied  climatic  records 
to  which  the  model  must  use  in  simulating  the  lake  records, 
giving  a  better  probability  of  predicting  the  lake- watershed 
when  variables  are  unknown  or  changed.  Also,  a  longer  run 
will  reduce  the  influence  of  the  choice  of  the  irtial 
parameters,  such  as  QGw?,  HELM?,  lake  stage,  or  lake  water 
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quality  on  later  results.  Experience  with  several  data  sets 
show  that  the  initial  conditions  only  inflence  the  initial  6 
or  7  simulation  periods. 

Ve c if icat io n  of  the  Model 

Verification  of  a  model  is  an  integral  part  of  the 
development.  This  procedure  ensures  that  the  model  has  no 
undetected  programming  or  logic  errors.  In  order  to  gain 
confidence  in  this  model,  it  was  tested  with  many 
representative  combinations  of  input  parameters.  The  testing 
was  carried  out  to  ensure  that  any  possible  set  of  climatic 
data,  unsaturated  zone  moisture  conditions  and  groundwater 
element  parameters  could  be  correctly  handled  by  the  model. 

The  model  was  also  verified  by  comparing  calculated 
results  with  three  sets  of  published  records.  The  procedure 
which  routes  moisture  through  the  unsaturated  zone  was 
tested  with  Ayers'  (  1 9  7  J )  precipitation,  potential 
evaporation  and  unsaturated  zone  moisture  storage  data.  The 
values  computed  by  this  mcdel  duplicated  the  values 
calculated  by  Ayers  (1970).  Bouting  of  moisture  through  tne 
groundwater  zone  was  verified  by  comparing  the  model's 
calculations,  with  results  presented  by  Dooge  (I960)  for  a 
set  of  test  data.  The  values  computed  by  the  model  were  the 
same  as  the  calculations  that  appeared  not  to  be  in  error. 

In  both  Ayers'  (1970)  and  Dooge's  (I960)  papers,  some  errors 
were  found  in  the  data  and  calculations.  The  correspondence 
between  their  calculations  and  the  model's  indicates  that 
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the  errors  in  Ayers  (1970)  and  Dooge  (I960)  were  not 
significant  in  the  calculations  ana  probably  resulted  from 
printing.  Ihe  ice  forming  model  was  checked  with  data  from 
Fertuck  et  al.  (1970)  and  ice  measurements  obtained  at 
Baptiste  Lake.  Ihe  calculated  values  compare  favorably  to 
the  measured  ice  thickness.  Results  are  listed  m  lanie  2. 


Table  2. 

VERIFICATION  OF  THE  ICE  FORKING  MODEL 


Data 

from  Fertuck  et  al. 

(  19  7  0) 

ice 

thickness 

p  re  c  l  p . 

temp. 

( actua 1) 

(model ) 

Dec: 

50  cm 

64.7  cm 

1.4  7  cm 

-16.4 

°C 

Jan : 

80  cm 

89.6  cm 

4.17  cm 

-27.  8 

°C 

Feb : 

92  cm 

105.4  cm 

1  *  j 8  cm 

-15.5 

°c 

Data  from  Baptiste 

Lake 

ice 

t  hickness 

p  r  e  c  i  p . 

tem  p. 

(actual ) 

(model) 

Nov : 

-? 

c  m 

22.9  cm 

0.  56 

cm 

-2.06 

°C 

Dec : 

32 

cm 

33.1  cm 

4.24 

cm 

-10.6  1 

°C 

Jan : 

43 

cm 

42. €  cm 

1.86 

cm 

-14.40 

°c 

Feb : 

51 

c  ro 

45.0  cm 

0.  69 

cm 

-0.90 

°c 

Mar: 

6  2 

cm 

49.6  cm 

2.  53 

cm 

-3.90 

°c 

III.  SENSITIVITY  ANALYSIS 


The  influence  of  the  hydrologic,  climatic  and 
physiographic  parameters  of  a  watershed  can  Le  studied  and 
their  significance  to  the  water-balance  can  be  evaluated 
using  the  mathematical  model.  This  type  or  study,  known  as  a 
sensitivity  analysis,  involves  the  systematic  variation  of 
one  or  more  model  parameters  while  the  rest  are  held 
constant. 

In  this  chapter,  nine  model  parameters  are  stuaieu  to 
examine  their  influence  on  the  water  levels  and  chemistry  of 
a  prairie  lake.  First,  the  influence  of  the  physical  bounds 
of  a  lake,  including,  variations  of  its  surface  area  (1)  and 
total  volume  (ii) ,  are  examined.  Next,  the  formation  of  ice 
(lii)  is  discussed  to  determine  the  functional  relationsnip 
between  the  salinity  of  prairie  lakes  and  parameters  or  size 
and  initial  total  dissolved  solids  content.  The  signif icance 
of  groundwater  quality  (iv)  in  controlling  the  salinity  or 
prairie  lakes  is  studied  with  the  model,  by  systematically 
varying  the  groundwater  inflow  chemistry.  Kali’s  influence  on 
natural  environment  of  a  watershed  is  examined  to  examine 
hypothetical  situations.  For  example,  when  land  is  cleared 
(v)  or  drained  (vi)  for  agricultural  activities,  or  the 
watershed  is  prog ressi vely  urbanized  (vii) .  Finally,  the 
influence  of  climatic  variacles  (precipitation  (vm)  and 
potential  evaporation  (ix)  )  is  discussed  with  emphasis  on 
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how  accurate  these  data  should  he  in  order  to  adequately 
model  a  lake-watershed . 

The  watershed-lake  system  usea  for  the  sensitivity 

analysis  was  designed  to  present  a  set  of  conditions  that 

would  re  typical  for  a  variety  of  lakes  on  the  prairies. 

Five  groundwater  storage  elements  (four  shallow  and  one 

deep)  are  defined  in  a  waters hed  with  a  land  area  of  300  km2 

and  a  lake  area  of  approximately  9  km2.  The  geological 

representation  of  the  five  groundwater  storage  elements  are: 

element  I  =  a  silt  unit  with  an  overlying  unsaturated  zone 

element  2  =  a  silt  unit  with  no  overlying  unsaturated  zone 

element  3  =  a  till  unit  witn  an  overlying  unsaturated  zone 

element  4  =  a  till  unit  with  no  overlying  unsaturated  zone 

element  5  =  a  deep  till  unit 

The  schematic  representation  of  the  watershed  and  values  or 
the  parameters  are  illustrated  on  Figure  9.  An  example  input 
card  deck  is  listed  in  Appendix  III.  All  simulations  were 
run  over  an  eight  year  time  period  using  a  time  step  of  one 
month . 

Lake  Surface  Area t  Lake  Voluae  and  Ice  Formation 

Lakes  and  reservoirs  in  the  Canadian  prairies  are 
important  water  resources.  However,  problems  such  as 
salinity,  contamination  or  declining  water  levels  have 
significantly  impaired  the  development  of  these  lakes.  In 
Alberta,  the  provincial  government  has  been  active  m 
research  to  try  to  reverse  this  trend  by  stabilizing  the 
lake  levels  and  improving  drainage  to  and  from  the  lakes 
(Laycock,  1968,  1971;  Gull  Lake  Task  Force,  1970;  Beid, 
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Schematic  representation  of  the  watershed 
model  used  in  the  sensitivity  analyses. 


Crowthers  and  Partners,  1973;  Stanley  and  Associates,  1974; 
EPEC  Consulting,  1S76;  Planning  Division,  1979) 

This  section  will  provide  insight  into  an  understanding 
of  how  lake  or  reservoir  physiography  alterations  influence 
water  quality  and  lake  levels.  Specifically,  the  following 
three  questions  will  te  considered; 

1.  How  do  water  quality  and  levels  vary  if  the  area  of  the 
lake  is  expanded  or  contracted  without  a  related 
increase  or  decrease  of  the  total  volume  of  lake  water? 

2.  What  is  the  impact  of  increasing  or  decreasing  the 
volume  of  water  in  the  lake  witnout  a  corresponding 
enlargement  or  reduction  in  lake  surface  area? 

3.  Does  the  formation  of  ice  on  lakes  significantly  effect 
the  salinity  of  lakes  during  tne  winter? 

Again,  it  is  assumed  that  no  salts  are  precipitated 
from  the  laxe.  In  a  real  lake,  the  pattern  of  salinity 
increase  will  ce  influenced  by  saturation  levels  with 
respect  to  several  minerals.  It  is  also  important  to  note 
that  the  saturation  level  is  dependent  on  the  temperature 
and  pH  of  the  lake  water.  Both  of  these  parameters  vary  over 
a  year  in  a  prairie  lake  and  hence  a  complete  analysis  of 
lake  salinity  changes  requires  that  a  model  incorporate  the 
saturation  processes.  The  lake  waters  considered  in  this 
study  are  typically  undecsatura ted  with  respect  to  many 
common  minerals.  Thus,  the  assumption  is  valid  as  a  first 
approximation. 

The  change  in  the  total  dissolved  solids  content  of  a 
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lake  with  different  surface  areas  is  shown  in  figure  10.  In 
each  case,  the  total  volume  of  the  lake  is  constant.  The 
most  obvious  trend  is  that  as  the  lake  area  to  lake  volume 
ratio  increases,  the  salinity  of  the  lake  is  increased  and 
the  variation  in  total  dissolved  solids  on  an  annual  basis 
is  enhanced,  kith  a  g r eater  lake  surface  area,  potential 
evaporation  and  precipitation  have  a  greater  inrluence, 
causing  an  overall  increase  in  salinity. 

Lakes  of  a  small  surface  area  undergo  less  annual 
elevation  fluctuations  than  dc  those  with  a  larger  area,  and 
the  overall  lake  surface  elevation  is  greater  (figure  10).  As 
the  surface  area  of  the  lake  decreases,  the  area  available 
for  evaporation  also  decreases,  and  hence  less  water  is  lost 
through  evaporation.  This  relationship  is  especially  evident 
during  the  winter  when  a  latqe  proportion  cf  lake  water 
forms  ice  in  a  shallow  lake. 

The  infuence  of  changing  lake  volumes  on  salinity  were 
simulated  by  adding  a  large  volume  of  water,  approximately 
equal  to  multiples  of  the  initial  lake  volume,  to  the  lake. 
The  water  is  assumed  to  be  added  ly  lowering  the  lake 
bottom,  thus  the  lake  surface  area  and  lake  level 
fluctuations  are  constant.  The  purpose  of  the  analysis  is  to 
study  the  effect  of  enlarging  the  volume  of  a  lake  or 
resevoir  while  its  surface  area  remains  the  same. 

The  increase  in  salinity  as  the  lake  volume  diminishes 
(figure  II)  is  due  to  groundwater  recharge  to  the  lake. 
Groundwater  inflow  to  the  lake  becomes  the  most  significant 
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Figure  10.  Analysis  of  various  lake  surface  areas  in  influencing  lake  levels  and  salinity. 
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Figure  11.  Analysis  of  various  lake  volumes  in  influencing  lake  salinity. 
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process  in  this  watecsaed,  especially  in  the  spring,  A  large 
volume  of  groundwater  mixing  with  a  smaller  volume  of  lake 
water  will  tend  to  alter  the  lake  salinity  towards  that  of 
the  groundwater. 

The  ratio  of  surface  area  to  total  lake  volume  is  very 
significant  in  controlling  both  the  salinity  of  prairie 
lakes  and  reservoirs  and  their  surface  elevations.  Thus  if 
one  were  physically  anle  to  alter  the  lake  by  concen trating 
the  same  volume  of  water  into  a  lake  with  a  smaller  surface 
area,  or  increasing  the  lake  volume  without  expanding  its 
surface  area,  the  salinity  of  the  lake  would  decrease  and 
the  lake  Jevels  would  increase.  This  is  tantamount  to 
deeping  a  lake.  In  reality,  increasing  the  la^e  surface 
elevation  will  make  the  surface  area  larger  and  the 
previously  discussed  controlling  factors  will  influence  the 
lake.  The  reason  for  the  increase  in  lake  salinity  is  that 
as  the  surface  area  to  lake  volume  ratio  increases, 
precipiation  and  potential  evaporation  exhibit  an 
increasingly  greater  influence.  Concentrating  water  into  a 
lake  with  a  smaller  surface  area  has  a  net  effect  of 
increasing  the  land  area  of  the  watershed  and  will  cause  the 
same  results.  These  results  correspond  to  Laycock's  (1968, 
1971)  lake  stabilization  studies,  in  which  he  proposed  that 
evaporation  losses  would  decrease  if  the  surface  area  of  the 
lake  availacle  for  evaporation  were  reduced  and  the  depth 
increased.  Hence,  the  lake  level  would  be  higher  and  the 
salinity  lower.  In  the  same  reports,  it  was  stated  that  by 
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expanding  a  watershed,  through  lake  or  stream  flew  diversion 
across  a  watershed  boundary,  higher  lake  levels  would 
result. 

Ice  formation  is  also  very  important  in  influencing  the 
salinity  of  lakes  during  the  winter,  especially  if  the  lake 
or  reservoir  is  saline ,  or  if  the  p  ro portion  of  lake  water 
that  is  taken  up  during  ice  formation  is  large  relative  to 
total  storage  of  lake  water.  A  greater  proportion  of  lake 
water  will  be  frozen  in  shallower  lakes  than  deep  lares 
(Figure  10  ana  11),  and  hence,  all  dissolved  lake  salts  will 
be  concentrated  in  a  smaller  volume  of  water. 

Groundwater  Chemistry 

There  are  two  reasons  for  conducting  a  sensitivity 
analysis  using  groundwater  chemistry  as  a  variable.  Firstly, 
tc  provide  insight  into  how  lake  salmity  may  be  controlled 
by  groundwater  chemistry.  Secondly,  an  analysis  can  help 
determine  how  accurately  this  component  should  te  defined  by 
field  measurements  in  order  to  adequately  model  the 
watershed.  Typically,  the  quality  of  groundwater  is  quite 
variable,  especially  in  the  shallow  groundwater  regime,  and 
is  hard  to  know  entirely. 

The  salinities  of  the  various  groundwater  storage 
elements  used  in  the  six  trials  are  listed  on  Table  3.  These 
salinities  are  typical  of  observed  groundwater  conditions  in 
the  prairies. 
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Table  3. 


CHEMISTRY  OF 

GROUNDS  ATEB  STORAGE 

ELEMENTS 

(all 

sal  init ies 

in  mq/L) 

trial  #1 

trial 

til 

trial 

#3 

element 

salinity 

element  salinity 

element  salinity 

1 

200 

1 

400 

1 

400 

2 

200 

400 

2 

800 

3 

25  C 

3 

500 

3 

5  00 

4 

2  50 

4 

50  0 

4 

500 

5 

300 

5 

600 

5 

600 

trial  #4 

trial 

#5 

trial 

tit 

element 

sa 1 i nit  y 

element  salinity 

element  salinity 

1 

400 

1 

500 

1 

8  00 

2 

400 

2 

500 

2 

800 

3 

1  OCO 

3 

625 

3 

1000 

4 

500 

4 

62  5 

4 

1000 

5 

600 

5 

750 

5 

1200 

Results  obtained  with  this  hypothetical  lake- wa tershed 
system  (Figure  12),  show  that  the  lake  is  not  particularly 
sensitive  to  even  larqe  changes  in  the  groundwater  chemistry 
of  part  of  the  basin,  tut  that  changes  in  the  groundwater 
chemistry  of  the  entire  basin  may  significantly  alter  the 
salinity  of  the  lake.  For  example,  a  100%  increase  in  the 
salinity  of  a  small  grcundwater  storage  element,  such  as  in 
run  #3,  or  a  larqe  element,  run  #4,  resulted  in  only  a 
minimal  change  to  the  lake.  Run  tit  was  conducted  with  the 
groundwater  salinity  100%  higher  in  all  storage  elements 
than  the  initial  run  (#2)  and  the  resultant  lake  salinity  is 
approximately  75 %  higher.  An  overall  increase  of  25 %  of  the 
groundwater  salinity  (run  #5)  raised  the  salinity  of  tne 
lake  by  approximately  20%. 


The  implications  of  this  analysis 
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Figure  12.  Analysis  of  various  groundwater  salinities  in  influencing  lake  salinity. 
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would  be  that  m  a  prairie  watershed  system,  similar  to  that 
modelled  heie,  influences  which  cause  an  increase  in  the 
groundwater  salinity  over  a  small  part  of  the  watershed 
would  not  signif icantiy  alter  the  chemistry  ox  the  lake. 

Land  Use 

Man  has  greatly  influenced  all  aspects  of  the 
environment.  Urban  development,  agriculture,  forestry  and 
land  reclamation  can,  within  a  time-span  of  a  few  years  or 
even  a  few  months,  afxect  all  systems.  This  section  will 
demonstrate  the  possible  impact  ot  two  common  activities 
within  a  watersned;  { i )  increasing  the  proportion  of  the 
watershed  that  is  covered  by  an  unsaturated  zone,  and  (n) 
increasing  the  area  or  the  watershed  tiiat  is  occupied  by  an 
impervious  area.  Laycock  (1971,  1973)  suggested  that  a 

change  in  the  land  from  forest  cover  to  cropland  or  pasture 
would  produce  more  surface  runoff  and  less  groundwater  flow 
to  a  lake.  Thus  lake  levels  wcula  increase  and  the  quality 
of  the  lane  would  be  enhanced,  assuming  that  erosion  not 
increase. 

The  effect  of  changing  the  land  cover  for 
lake- watershed  management  programs  is  worth  investigating 
via  sensitivity  analyses.  Unfortunately,  this  requires  data 
on  vegetation  types,  rooting  depths  and  moisture  use  by 
various  plant  species,  which  is  beyond  the  scope  of  this 
model.  However,  the  modelling  of  land  use  changes  can  be 
attempted  by  the  general  grouping  of  land  covers  into 
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forest,  wet-lands,  cropland  and  pasture.  These  various 
categories  affect  the  areai  extent  of  the  unsaturated  zone 
and  its  moisture  capacity.  Therefore,  two  analyses  shall  ae 
conducted  by:  (a)  simulating  the  draining  of  wet-lands  by 
increasing  the  portion  of  the  watershed  that  is  occupied  by 
an  unsaturated  zone,  and,  (ii)  studying  the  change  from 
forest  cover  to  pasture  or  cropland  by  reducing  the  moisture 
capacity  of  the  unsaturated  zone  (Laycock,  J9t7,  1966,  1971; 

Gray  et  ai. ,  1970) . 

Decreasing  the  unsaturated  zone  moisture  capacity  in 
the  model  resulted  in  greater  recharge  to  the  groundwater 
storage  elements,  and  hence,  greater  groundwater  discharge 
to  the  lake  because  less  precipitation  was  retained  in  the 
unsaturated  zone.  Surface  runoff  increases  because  of  the 
tendency  for  percolation  to  exceed  the  maximum  groundwater 
infiltration  rate  during  periods  of  intense  precipitation. 

Of  these  two  factors,  groundwater  discharge  is  the  more 
prominent  and  acts  to  raise  the  lake  levels  and  increase  the 
salinity  (Figure  13).  A  noticeable  change  in  luxe  levels  did 
not  result  because  water  had  not  significantly  exceeded  the 
unsaturated  zone  moisture  capacity. 

The  increased  surface  runoff  and  increased  lake  levels, 
cor respondi ng  to  a  decreasing  unsaturated  zone  moisture 
capacity,  is  in  agreement  with  Laycock's  (1971,  1973) 

findings.  However,  the  increased  groundwater  discharge  and 
greater  lake  salinity  are  opposite  to  his  findings.  The 
reason  for  increased  groundwater  discharge  found  in  this 
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Figure  13-  Analysis  of  various  unsaturated  zone  moisture  capacities  in  influencing  lake 
1  eve  1 s  and  sal i n i t y . 


watershed  is  that  lore  water  is  able  to  recharge  tiie 
groundwater  zone;  hence,  more  is  discharged  to  the  lake  as 
groundwater,  raising  tne  salinity  of  tne  lake.  Should  a 
decrease  in  the  unsaturated  zone  moisture  capacity  be 
accompanied  by  a  larger  groundwater  recharge  rate,  the 
groundwater  discharge  and  lake  salinity  would  increase  as 
suggested  by  Laycock  (1971,  1973). 

The  next  part  of  this  sensitivity  analysis  will  examine 
the  influence  that  increasing  the  portion  of  the  rasin  which 
is  occupied  by  an  unsaturated  zone  has  on  lake  levels  and 
salinity.  The  analagous  real  situation  would  be  the  draining 
of  wet-lands,  which  are  lands  occupied  by  slcuqhs,  swamps, 
muskeg  or  high  watertanles.  It  should  be  noted  that  when  a 
wet-area  is  drained,  the  water  wrll  be  diverted  directly  to 
a  lake  or  a  river.  The  increased  surface  runoff  will  produce 
higher  lake  levels  and  increase  water  flow  through  the  lake, 
resulting  in  improved  lake  water  quality  (Laycock,  196d, 
1971).  However,  in  this  study  only  the  unsaturated  zone  area 
is  enlarged;  the  drainage  from  former  wet-lands  does  not 
increase.  It  is  assumed  that  the  land  has  been  drained  for  a 
sufficient  length  of  time  for  the  watershed  to  reach  an 
equilibrium  with  the  increased  area  underlain  by  an 
unsaturated  zone. 

Increasing  the  proportion  of  the  watershed  that  is 
covered  by  an  unsaturated  zone  reduces  the  salinity  of  the 
lake  and  lowers  the  surface  elevations  (Figure  14).  It  is 
proposed  that  the  total  recharge  to  the  lake  is  decreased  as 
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igure  lA.  Analysis  of  various  areal  extents  of  an  unsaturated  zone  in  influencing  lake 
1  eve  1 s  and  salinity. 
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the  extent  of  the  unsaturated  zone  is  enlarged  because  less 
water  (especially  groundwater)  jls  discharging  to  t re  lake. 
Lake  water  chemistry  responds  in  the  expected  manner  with  a 
decrease  in  the  total  dissolved  solids.  In  this  example, 
groundwater  is  the  major  contributor  of  dissolved  solids  to 
the  lake  and  a  decrease  of  groundwater  inflow  to  the  lake 
will  significantly  reduce  lake  salinity.  Also,  fluctuations 
of  the  lake  water  chemistry  decrease  as  the  unsaturated  zone 
area  increases  (Figure  14).  The  greater  variability  m 
chemistry  among  simulations  during  the  winter  is  due  to  the 
effect  that  the  formation  of  ice  nas  on  different  lake 
salinities.  Thus,  the  net  effect  of  draining  wet- areas,  as 
indicated  by  tnis  study  (decreased  lake  levels  and 
salinity) ,  would  appear  to  be  due  to  reduced  groundwater 
recharge.  It  can  be  concluded  that  as  the  area  ol  toe 
watershed  that  is  covered  by  an  unsaturated  zone  increases, 
less  water  will  reach  the  groundwater  system  or  discharge  as 
surface  runoff  because  the  water  re gui re me n t s  ol  the 
unsaturated  zone  must  first  be  fulfilled  before  any 
groundwater  recharge  occurs.  Another  implication  of  this 
analysis  is  that  the  lake  is  quite  sensitive  to  the 
proportion  of  the  basin  that  is  occupied  by  an  unsaturated 
zone  and  this  parameter  should  therefore  be  defined 
carefully.  In  this  watershed,  a  10a  increase  m  the 
unsaturated  zone  area  resulted  in  an  approximate  15  mg/L 
drop  in  lake  salinity.  Similarly,  lake  levels  ate  much 
higher  in  the  spring  if  the  unsaturated  zone  is  not  very 
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large,  due  tc  surface  runoff  and  groundwater  discharge. 

The  roam  effects  of  urbanization  on  a  watershed  are 
primarily  an  increase  m  the  amount  of  impervious  area  by 
the  construction  of  roads  and  Buildings.  This  causes  a 
decrease  in  the  time  required  for  runoff  and  contaminants  to 
reach  a  lake  or  river  through  storm  sewers  and  lined 
drainage  ditches  (Crawford  and  Linsley,  19t>6).  Water  falling 
on  these  impervious  areas  will  ne  directly  discharged  to  a 
lake  or  a  river.  As  previously  stated,  isolated  impervious 
areas  are  not  considered  to  be  part  of  the  impervious  area 
parameter  because  surface  runoff  ircm  here  must  flow  over 
pervious  land  before  reaching  surface  water  bodies. 

The  i  tract  that  increasing  the  portion  of  the  watershed 
that  is  impervious  has  on  the  lame  stage  and  lake  salinity 
is  illustrated  on  Figures  15.  As  greater  porporticns  of  the 
watershed  are  covered  by  impervious  areas,  more 
precipitation  is  routed  directly  to  the  lake,  with  a 
resulting  decrease  in  lake  salinity.  Even  a  small  increase 
of  the  impervious  area,  involving  as  little  as  l>  of  the 
watershed  in  this  example,  caused  a  decrease  of  the  lake 
salinity  by  10  mg/L.  The  reason  tor  this  is  that  surface 
runoff  becomes  increasingly  important  relative  to 
groundwater.  The  lake  level  is  higher  and  the  fluctuations 
are  greater  when  the  impervious  area  increases. 

The  results  of  tn is  analysis  correspond  to  Crawford  and 
Linsley' s  (1966)  experience  with  several  watersheds, 
suggesting  that  the  influence  of  impervious  areas  found  in 
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Figure  15-  Analysis  of  various  impervious  areas  in  influencing  lake  levels  and  salinity. 
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this  study  is  common.  Precipi t at ion  is  no  longer  routed  to 
the  lake  via  groundwater  storage  elements,  rut  becomes 
surface  runoff.  Thus  the  flow  length  has  been  reduced  and 
the  time  delay  for  water  to  reach  the  lake  decreased, 
enabling  the  surface  runoff  from  a  given  storm  to  reach  the 
lake  at  faster  rate  and  in  greater  amounts. 

The  sensitivity  analysis  indicates  that  the  impervious 
cover  is  an  important  watershed  parameter  controlling  lake 
salinity.  Accurate  values  should  re  oitamcd  Ly  either  field 
study  or  a  detailed  calibration  procedure. 

Precipitation  and  Po te n tial  Evaporation 

For  many  watersheds  cn  the  prairies,  there  are  often  no 
meterological  recording  stations  located  in  the  specific 
areas  being  studied.  Because  data  must  be  interpolated  or 
extrapolated  from  the  nearest  stations,  they  are  accompanied 
by  uncertainties.  Parameters  that  are  not  always  measurable, 
such  as  potential  evaporation  from  lakes,  reservoirs, 
forested  land  and  agricultural  areas,  can  ne  estimated  by  a 
number  of  formulae  and  methods.  The  accuracy  of  these 
estimates  is  often  douttful  because  of  the  many  influencing 
factors.  This  part  of  the  sensivity  analysis  will  examine 
the  level  of  accuracy  required  for  precipitation  and 
potential  evaporation  data  in  order  to  adequately  model  lake 
level  and  salinity  fluctuations. 

The  effect  on  lake  levels  of  over-estimating  or 
under-estimating  precipitation  and  potential  evaporation  by 


up  to  25/t  in  various  sizes  of  watershed s,  are  shown  or 
Figures  16n  arid  d,  17b  and  d,  18d  and  d.  As  expected,  as 
potential  evaporation  increases,  the  elevation  of  the  lake 
surface  declines  and  increasing  precipitation  will  cause  the 
lake  level  to  rise. 

The  changes  of  the  lake  level  are  also  dependent  on  tne 
size  of  the  watershed.  When  precipitation  varies  by  as  much 
as  25%  on  the  large  watersheds  used  in  this  example  (greater 
than  203  km2),  the  level  may  only  change  by  a  small  amount; 
less  than  .25  m,  as  shown  on  Figures  17b  aid  18b.  But  when 
the  small  (50  km2),  watershed  was  tested,  a  25k  variation 
resulted  in  a  lake  level  change  or  almost  2  m  (Figure  16b). 
In  cases  where  potential  evaporation  was  varied,  a  similar 
trend  could  be  observed.  A  deviation  of  approximately  .05m 
occurs  in  tne  large  watersneds  (Figures  1 7 d  and  18d),  but 
the  smaller  watershed  shows  that  a  25k  change  resulted  in  a 
lake  level  variation  of  almost  1  m  (Figure  16d).  It  should 
also  be  noted  that  the  deviation  from  the  normal  lake  level 
in  a  small  watershed  is  greater  when  precipitation  or 
potential  evaporation  are  under-estimated,  than  if  they  are 
over-estimated. 

The  effect  on  lake  chemistry  of  changing  these 
parameters  is  more  difficult  to  explain.  Another  parameter, 
watershed  size  plays  an  important  role  in  determining  lake 
salinity.  Figures  16a  and  16c  show  the  result  of  increasing 
or  decreasing  precipitation  and  potential  evaporation  by  up 
to  25k  in  the  small  watershed  used  in  the  previous  analysis. 
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Figure  16.  Analysis  of  varying  precipitation  and  potential  evaporation,  separately,  in  influencing 
lake  levels  and  salinity  in  a  50  km2  watershed. 
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igure  16.  Analysis  of  varying  precipitation  and  potential  evaporation,  separately,  i 
influencing  lake  levels  and  salinity  in  a  50  km2  watershed. 
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Figure  17*  Analysis  of  varying  precipitation  and  potential  evaporation,  separately,  in 
influencing  lake  levels  and  salinity  in  a  200  km?  watershed. 
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Figure  17-  Analysis  of  varying  precipitation  and  potential  evaporation,  separately,  in 
influencing  lake  levels  and  salinity  in  a  200  km2  watershed. 
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Figure  18.  Analysis  of  varying  precipitation  and  potential  evaporation,  separately,  in 
influencing  lake  levels  and  salinity  in  a  300  km2  watershed. 
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Figure  18.  Analysis  of  varying  precipitation  and  potential  evaporation,  separately,  i 
influencing  lake  levels  and  salinity  in  a  300  km^  watershed. 
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It  is  net  surprising  to  otserve  that  as  precipitation 
increases,  lake  salinity  is  lowered,  and  that  as  potential 
evaporation  increases,  lake  salinity  is  raised,  indicating 
that  precipitation  and  potential  evaporation  from  the  lake 
surface  become  the  most  important  processes  controlling  the 
chemistry  of  the  lake.  However,  as  the  area  of  the  test 
watershed  increases  relative  to  the  lake  area,  groundwater 
discharge  becomes  a  more  significant  factor  controlling  lake 
chemistry.  Because  the  greundwater  of  the  test  watersheds 
have  a  relatively  high  salinity,  a  greater  contribution  of 
groundwater  to  the  lake,  {relative  to  other  sources  of  water 
such  as  surface  runoff  and  precipitation  falling  directly  on 
the  lake),  will  cause  the  salinity  of  the  lake  to  approach 
that  of  the  groundwater.  Increasing  potential  evaporation  or 
decreasing  precipitation  over  a  large  watershed  results  in  a 
lower  lake  salinity  (Figure  18a  and  16c),  because  there  is 
less  groundwater  being  discharged  to  the  lake.  A  similar 
effect  an  the  large  watershed  can  be  noted  when 
precipitation  is  increased  by  25a  and  potential  evaporation 
is  reduced  ry  25%;  groundwater  recharge  is  increased  ana 
therefore,  groundwater  discharge  is  also  greater, 
contributing  more  total  dissolved  solids  to  the  lake  and 
increasing  its  salinity. 

In  the  previous  sensitivity  analyses,  either 
precipitation  or  potential  evaporation  were  varied.  However, 
the  more  realistic  situation  is  for  them  tc  vary  together. 

If  precipitation  increases,  there  are  fewer  days  in  which 
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evaporation  may  occur.  Thus  potential  evaporation  would 
decrease  and  vice  versa.  Figure  19  illustrates  the  results 
fcr  model  trials  for  a  15$  increase  and  decrease  in 
precipitation  with  a  respective  15;a  decrease  and  increase  in 
potential  evaporation.  Here,  the  lake  surface  varies  up  to 
,15m  and  the  lake  water  quality  varies  by  only  25  mg/L. 

After  conducting  this  sensitivity  analysis,  the 
following  conclusions  can  ce  stated  for  these  and  similar 
wa  ter sheds : 

1.  An  error  in  estimating  precipitation  and  potential 

evaporation  may  not  significantly  alter  the  calculated 
lake  salinity  and  lake  surface  elevations  or  a  large 
lake- watershed,  out  over-estimating  or  under-estimating 
precipitation  and  potential  evaporation  may  cause 
significant  deviations  in  a  small  watershed  (less  than 
100  km2).  In  particular,  under-estimating  precipitation 
or  over-estimating  potential  evaporations  causes  a  much 
larger  lake  level  change  than  under-  and  over-estimating 
potential  evaporation  and  precipitation,  respectively. 
Therefore,  in  watersheds  similar  to  those  tested  here, 
it  may  be  concluded  firstly,  that  more  accuracy  is 
required  for  a  small  lake- water  sheds  than  for  a  large 
one,  and,  secondly,  that  if  the  actual  values  for  a 
small  watershed  are  unknown,  better  lake  level  and 
salinity  simulations  may  be  obtained  when  potential 
evaporation  is  under-estimated  and  precipitation  is 
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igure  19*  Analysis  of  varying  precipitation  and  potential  evaporat ion  ,  together  , 
influencing  lake  levels  and  salinity. 
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The  results  of  the  simulations  car  he  gerer alizaile  in 
saying  that  precipitation  to,  and  potential  evaporation 
from  a  lake  in  a  small  watershed  plays  a  more  prominent 
role  in  controlling  lake  salinity  than  does  groundwater 
discharge.  This  conclusion  assumes  that  the  lake  is  the 
major  groundwater  discharge  point  for  the  watershed  and 
that  as  the  area  of  the  watershed  varied,  in  these 
tests,  the  area  of  the  lake  remained  essentially 
constant.  It  follows  tnat  in  a  larger  watershed  where 
groundwater  is  a  major  contributor  to  the  lake, 
groundwater  discharge  is  the  main  process  controlling 
lake  salinity.  In  this  case,  decreasing  potential 
evaporation  or  increasing  precipitation  will  increase 
groundwater  recharge,  hence,  increased  groundwater 
discharge  to  tne  lake.  This  groundwater  inflow  will 
increase  tne  salinity  of  the  lake. 
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IV.  APPLICATION  TO  THE  BAPTISTE  LAKE  HATEESHED 


This  chapter  will  demonstrate  a  practical  application 
of  the  model  to  the  Baptiste  Lake  watershed.  Important 
features  of  the  basin  such  as  physiography,  qeoloqy, 
hydrogeology  and  hydrology  will  re  reviewed  in  the  following 
sections.  The  purpose  of  this  description  is  to  give  the 
reader  a  fetter  insight  into  the  cnoice  or  parameters  for 
the  simulation.  Following  this  discussion  is  an  analysis  or 
the  simulation.  Finally,  a  few  specific  applications  of  the 
model  to  the  Baptiste  Lake  watershed  are  presented. 

General  D isc r ipti  on  of  the  jfatershed 

Baptiste  Lake  is  located  in  the  Eastern  Alberta  Plains 
physiographic  region,  approximately  130  km  north  of 
Edmonton,  Alberta  (54°  45'  N,  113°  33*  *)  .  The  watershed 
(Figure  20)  covers  an  area  of  approximately  306.9  km2,  of 
which  53%  is  forest,  25%  is  standing  water  consisting  of 
Baptiste  lake,  swamps  and  rivers,  16%  is  agricultural  ana 
less  than  !%  is  urban;  mainly  cottage  development  along  tne 
shore  of  the  lake.  The  laud  tc  the  north  and  east  of  the 
lake  is  primarily  rarm  land,  while  that  to  the  west  and 
south  is  primarily  swamp  and  forest. 

Land  elevations  range  from  578  m  above  mean  sea  level 
at  the  surface  of  Baptiste  Lake,  to  over  fc90  m  a.m.s.i,  west 
of  the  lake.  The  area  of  the  lake  is  approximately  9.2  km2 
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with  a  maximum  depth  of  27.5  m  and  a  volume  of  6.5  *  107  m3 
(±  10%)  at  an  elevation  of  578.35  m  a.m.s.l.  Twelve 

tributary  streams  flow  into  the  lake.  One  outlet ,  Baptiste 
Creek,  drains  the  lake  and  flows  east  to  the  Athabasca 
River. 

Field  kor k 

Two  programs  involving  deep  test-hole  drilling  were 
undertaken  tc  investigate  tne  regional  geology  and 
hydrogeology  and  to  determine  whether  the  buried  preqlacial 
bedrock  valley  reported  by  Eorneuf  (1973)  is  present  beneath 
Baptiste  Lake.  Thirteen  test-holes  were  drilled  with  a 
rotary  rig  through  the  glacial  drift  and  into  the  upper  rew 
meters  of  bedrock  (Figure  21).  All  boreholes  were  logged 
geophysically.  A  summary  cf  the  results  of  this  program  is 
presented  in  Appendix  VIII.  Headers  interested  in  a  more 
complete  description  of  the  pr eg  rams  can  refer  to  Crowe 
{ 1 S73)  . 

Three  observation  wells  were  completed  in  the  sand  and 
gravel  deposits  of  the  buried  valley.  These  welis  were  used 
tc  monitor  shallow  groundwater  levels  and  to  obtain  samples 
for  chemical  analysis.  Pumping  tests  were  conducted  with  tne 
wells  to  determine  hydraulic  properties  of  the  major  aquifer 
in  both  the  north  and  south  buried  valleys.  The  hydraulic 
head  was  recorded  (Appendix  VIII),  and  groundwater  samples 
obtained  from  August  1377  to  August  1S76. 

During  early  May  1978,  25  shallow  test-holes  were 
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Figure  21.  Location  of  the  test  holes 
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augered.  Water  table  wells  were  completed  at  23  sites 
(Figure  21).  These  veils  were  used  to  monitor  groundwater 
level  fluctuations  and  chemistry.  The  results  of  the  shallow 
drilling  program  are  listed  in  Appendix  VIII,  The  water 
table  observation  wells  were  completed  with  a  large  slotted 
interval  to  ensure  that  the  water  tarle  was  intersected. 

Once  the  wells  nad  stabilized,  water  levels  were  monitored 
during  the  summer  and  fall  of  1978  (Appendix  VIII).  During 
the  period  of  study,  two  sets  of  water  samples  were 
collected  for  analysis  by  the  Pollution  Control  Division, 
Alberta  Environment  Laboratory,  or  at  the  University  of 
Alberta. 


Geology 

The  uppermost  bedrock  unit  underlying  Baptiste  Lake  is 
the  LaBiche  Formation  (Femak,  1944  ;  Green,  1972),  which  is 
Upper  Creataceous  in  age.  It  consists  of  dark  greenish-grey 
marine  shale.  The  hills  west  of  Baptiste  lake  consist  cf  a 
younger  unit,  the  Upper  Cretaceous  Wapiti  Formation  (Green, 
1972),  a  unit  composed  of  argillaceous  sandstone,  bentonitic 
mudstone  and  coal. 

The  bedrock  topography  of  the  area  has  been 
re-interpreted  from  the  studies  by  Borneuf  (1972)  and 
Carlson  (1976)  to  include  data  from  the  present  test 
drilling  program  and  from  water  well  logs  (Figure  22), 
Baptiste  Lake  is  situated  in  a  deep  bedrock  depression  which 
is  part  of  a  preglacial  buried  bedrock  valley  (Figure  22). 
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Figure  22. 


Bedrock  topography,  Baptiste  Lake  area 
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This  valley  trends  from  west  to  east  at  an  elevation  oi 
approximately  525  m  above  mean  sea  level  (15  m  below  the 
lake  bottom) .  Bedrock  elevations  in  the  surrounding  area 
typically  exceed  575  m  above  sea  level  west  of  the  lake  and 
56C  ir  east. 

Preglacial  sand  and  gravel  is  found  resting  on  the 
bedrock  surface  at  tne  base  of  the  buried  valleys.  These 
deposits  consist  of  fine  sand  to  medium  gravel,  and  contain 
a  considerable  amount  of  clay. 

Quaternary  deposits  were  originally  mapped  on  a  very 
large  scale  by  Kjearsgaard  (  1  972).  Additional  data  from  the 
shallow  drilling  program  and  from  air  photo  interpretation 
facilitated  a  more  detailed  interpretation  or  the  surficiai 
geology  (Figure  23).  The  till,  lake  deposits,  alluvium  and 
dune  deposits  found  in  the  area  are  often  mixed  or  reworked. 

Till  is  the  dominant  material.  It  overlies  the  bedrock, 
generally  with  thicknesses  of  up  to  30  ®.  however,  overlying 
the  preglacial  deposits  in  the  bedrock  valleys,  the  till 
attains  thicknesses  of  up  to  100  m.  The  dark  brown-grey  to 
dark  grey  till  lias  a  clay-loam  texture  with  occasional 
interbeds  of  sand  or  gravel.  The  parent  material  of  the  till 
is  the  LaBiche  formation  (Kjearsgaard,  1972).  The  till  is 
described  ty  Kjearsgaard  (1972)  as  an  ablation  or 
recessional  type  with  associated  iandforms  such  as  ice 
contact  rings  on  the  west  side  of  the  lake  and  knob  and 
kettle  topography  near  the  southeast  end  of  the  lake.  Along 
the  west  shore  of  Baptiste  lake,  the  till  has  been  reworked 
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Figure  23.  Surficial  deposits,  Baptiste  Lake  area. 
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by  glacial  ice  and  melt-water.  It  is  silty  and  contains  some 
sand  lenses. 

Lacustrine  deposits  are  founa  generally  as  a  thin  unit 
{less  than  0,5  m)  inantlinq  till.  In  depressions,  they  can  be 
up  to  3  m  thick. 

Along  the  northwest  portion  of  the  study  area  which 
overlies  the  tuned  valley,  the  unconsolidated  sediments 
consist  of  layers  of  till,  alluvium  and  silt.  At  the  surface 
they  have  often  been  reworked  and  mixed,  and  are 
cummulatively  described  in  this  thesis  as 
alluvium-lacustrine  de posits  (F igure  23). 

Medium  to  coarse  grained,  rust-orange  aeolian  sands  are 
found  southwest  of  the  lake  and  with  a  minor  occurrence  on 
the  east  side  of  the  south  basin,  Ihe  maximum  thickness 
recorded  during  this  study  was  30  m  at  1-20-66-24-&4M. 

Hydro -geology 

Ihe  patterns  of  groundwater  flow  were  investigated 
using  data  obtained  from  the  deep  and  shallow  drilling 
programs,  from  the  piezometers  and  from  existing  well  logs. 
However,  most  test-holes  drilled  during  this  study  and 
domestic  wells  are  located  close  to  the  lake  and  the  major 
portion  of  the  western  half  of  the  watershed  has  not  been 
drilled. 

Four  hy drost rat igraphic  units  have  been  defined:  (1) 
sand  and  gravel  deposits  situated  alonq  the  base  of  the 
buried  preglacial  valleys,  (2)  the  major  till  unit,  covering 
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most  of  the  basin,  (3)  a  silty  aliuvium-lacust rire  unit 
overlying  till  and,  (4)  aeolian  sands  which  occupy  the 
smallest  portion  of  the  basin.  These  units  correspond  to  tne 
surficial  deposits  previously  discussed.  The  Labiche 
Formation  is  of  little  hydrogeological  significance.  No 
existing  domestic  wells  are  completed  in  this  bedrock  unit 
and  test  drilling  did  not  define  any  potential  aguifers.  The 
Wapiti  Formation,  located  west  of  Baptiste  Lake,  contains 
more  significant  yields  of  groundwater  (Borneuf,  1972). 

Groundwater  flow  in  tne  Baptiste  Lake  area  can  be 
studied  by  examining  the  hydraulic  conductivity  of  the 
hy drostr atigraph ic  units.  Values  of  hydraulic  conductivity 
of  the  buried  channel  deposits  were  determined  from  aguifer 
tests.  Hydraulic  conductivities  for  the  other  units  were 
obtained  by  estimation.  The  till,  which  probably  contributes 
most  of  the  groundwater  to  the  lake,  was  estimated  to  have  a 
hydraulic  conductivity  in  the  order  of  I 0~ 6  cm/sec.  The 
hydraulic  conductivity  of  the  alluvium-lacustrine  deposits 
is  an  average  for  the  total  thickness  or  the  unit  and  is 
estimated  to  be  10~4  cm/sec.  The  aeolian  sands  possess  the 
highest  hydraulic  coad ucti vit y ,  being  in  the  order  of 
1C-1  cm/sec.  The  hydraulic  conductivity  of  the  sand  and 
gravel  in  both  the  north  and  south  bedrock  valleys  were 
found  to  be  approximately  1.2*10-4  cm/sec  and 
1.6*10~4  cm/sec  respectively.  Thus  water  movement  through 
the  till,  including  infiltration,  is  generally  guite  slow, 
whereas  flow  through  the  aeolian  sands  will  be  relatively 
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rapid. 

The  water  table  was  found  at  depths  ranging  from  0  ru  to 
15  m  below  ground  surface,  and  fluctuated  according  to 
recharge  from  precipitation.  The  depth  to  water  was  usually 
greatest  in  areas  of  extensive  till,  because  the  uplands 
west  of  the  lake  consist  mainly  of  swamps  and  muskeg,  the 
water  table  is  at  the  ground  surface.  Along  tne  trend  of  the 
preglacial  valley,  the  watertarle  is  usually  at  or  near  tne 
ground  surface. 

It  is  suspected  that  the  main  groundwater  recharge  area 
for  the  Baptiste  Lake  watershed  is  the  upland  area  to  the 
west  and  southwest  of  the  lake.  Water  is  either  ponded  in 
shallow  surface  depressions,  which  drain  to  the  lake,  or  it 
slowly  percolates  into  the  till.  Water  in  the  shallow 
hyd rostra ti g raphic  units  primarily  discharge  to  streams  or 
directly  to  Eaptiste  Lake.  Some  water  probably  percolates  to 
the  deeper  groundwater  systems,  which  primarily  discharge  to 
Baptiste  Lake.  Groundwater  from  the  deeper  till  units  cay 
possibly  flew  into  the  buried  channel  deposit.  The  chemistry 
of  groundwater  in  the  turied  valley  sediments  indicates  that 
it  probably  originates  in  the  Wapiti  formation  because  of 
the  similarity  of  the  chemistry  of  groundwaters  from  these 
two  units. 

Rater  levels  at  piezometers  IB,  2B  and  3B  remained 
relatively  constant  during  the  period  of  study,  varying  only 
by  a  maximum  of  1.5  m  at  piezometer  2B  and  2  m  at  piezometer 
IB.  The  hydraulic  head  of  the  wells  completed  in  the  buried 
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valley  deposits  on  the  west  side  of  the  lake  ate  above  the 
ground  surface  in  all  three  observation  wells  and  at  three 
domestic  wells.  However,  hydraulic  heads  in  domestic  wells 
completed  in  these  sediments  on  the  east  side  or  the  lake, 
at  NW -36-66-24-W4 PI ,  are  con sider all y  lower.  This  head  loss 
between  piezometers  IB,  2E,  3E,  the  lake  and  the  wells  at 
N W-36-6 6- 2 4- Whtf ,  is  interpreted  as  groundwater  discharge 
occurinq  from  the  buried  valley  to  the  lake  and  outflow  from 
the  lake  via  the  buried  valley.  On  the  basis  of  the  shallow 
groundwater  chemistry,  it  is  thought  that  some  qroundwater 
moves  upward  from  the  turied  valley  through  the  overlying 
till  at  the  southwest  and  northwest  ends  of  Baptiste  Lake. 

Groundwater  Chemistry 

Chemical  analyses  were  conducted  on  groundwater  samples 
obtained  from  wells  and  from  piezometers  completed  in  the. 
sucficial  deposits  (Appendix  VII).  The  chemistry  of  the 
shallow  groundwater  regime  exhibits  extreme  variability  with 
no  discernable  patterns.  This  situation  is  typical  of  many 
shallow  groundwater  systems  in  Alberta  (Schwartz,  1979). 
However,  when  the  analyses  were  grouped  according  to  the 
four  hy drost ratigraphic  urits  in  which  the  well  or 
piezometer  was  completed,  the  mean  ionic  concentrations, 
summarized  in  Table  4,  show  seme  patterns. 

Groundwaters  in  the  turied  valley  deposits  and  the 
overlying  surficial  deposits  have  very  different  chemical 
character is t ics.  Groundwater  from  the  ruried  valley  is 
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characterized  ny  high  concen t ra tions  of  soaiut  and  low 
calcium  and  magnesium,  waters  from  the  surficial  materials 
typically  have  higher  calcium  and  magnesium  concentrations 
and  much  lower  sodium  concentrations.  Bicarbonate  is  the 
major  anion  present  in  both  deep  and  shallow  groundwater. 

1  he  chemistry  of  the  buried  valley  deposits  may  be  the 
result  of  discharge  of  groundwater  from  the  fcapiti  Formation 
into  sand  and  gravel  situated  in  the  buried  valley. 

Table  4. 

MEAN  IONIC  VALUES  Of  GRCUNLU ATEb  IN  THE 

LAKE  lATEBSHED 

(all  values  in  mg/L) 


Unit 

Ca 

Mg 

Na 

K 

HC  03 

5C4 

Cl 

ID  5 

Till 

98 

33 

63 

5.  1 

630 

74 

2 

1060 

Alluvium-Lac. 

10  3 

27 

16 

3.  0 

433 

5  1 

1 

735 

Aeolian 

1  1  9 

10 

12 

5.  C 

16  7 

1  8 

<  1 

300 

Buried  Valley:  N 

3  7 

14 

357 

3.  1 

6  8  7 

223 

98 

1120 

S 

<2 

2 

305 

1.  6 

6  5 

6b 

24 

757 

Var ia hi  lity 

in  chemistry 

is 

al  so 

evident 

in  the 

groundwater  within  the  buried  valleys  flowing  towards 
Baptiste  Lake.  Higher  concentrations  of  all  major  ions  are 
observed  at  the  north  valley  as  compared  to  the  south. 
Chloride  and  sulfate  are  especially  more  prevalent  in  the 
north  valley.  Total  dissolved  solids  is  approximately  1120 
mg/L  in  the  north  channel  as  opposed  to  with  760  mg/L  in  the 
south  channel. 

Groundwater  in  the  glacial  deposits  are  generally 
character 1st ic  of  shallow  groundwaters  in  the  prairies; 
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being  of  a  calciu m-magnesiu m-t icarbonate  type.  Analyses  or 
waters  from  the  aeolian  sands  have  the  lowest  concen trations 
of  dissolved  chemical  species  of  any  unit.  This  is 
especially  evident  at  water  table  well  #10  where  the 
groundwater  analyses  indicated  a  total  dissolved  solids 
concentration  of  57  mg/L.  it  is  thought  that  these  waters 
are  derived  from  precipitation  whicn  has  rapidly  infiltrated 
into  the  sands  and  has  not  had  a  sufficient  residence  time 
to  dissolve  much  salt.  Groundwater  chemistry  from  the  till 
and  alluvium-lacustrine  deposits  are  guite  similar.  However, 
the  till  is  generally  a  few  mg/L  higher  in  all  ionic  species 
except  calcium.  In  the  till,  the  average  sodium  and  calcium 
concentra t ions  are  almost  the  same. 

Study  of  1  a p t i s t e  Lake 

The  principle  purpose  of  any  model  is  in  its 
application  to  understanding  the  processes  involved  in  a 
real  system,  and  ultimately  tc  predictmq  the  responses  of  a 
system  to  various  kinds  of  stresses.  In  this  section,  the 
lake- watershed  model  is  applied  to  the  Baptiste  Lake 
watershed.  In  particular,  this  section  will  consider  the  use 
of  the  model  in  studying  various  processes  in  the 
lake-wa te rshed  system,  such  as  the  nature  of  groundwater 
inflow  to  the  lake,  and  in  predicting  the  impact  of  future 
cottage  development  along  the  shoreline.  Finally,  data 
reguirements  for  lake- watershed  studies  or  this  type  will  be 
discussed  in  a  general  way  to  assess  the  possi bill  ties  for 
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future  model  studies. 

Eaptiste  La*.e  Simulation  Parameters 

To  achieve  a  satisfactory  simulation  of  any 
lake- watershed,  the  model  must  re  atle  to  approximate  the 
hydrologic  processes  in  the  basin.  The  hydrological 
variables  which  the  model  simulates  are  lake  staqe,  salinity 
of  the  lake,  outflow  from  the  lake  via  groundwater  and 
surface  water,  recharge  to  the  groundwater  zone,  unsaturated 
zone  moisture,  thickness  of  ice  during  the  winter,  surface 
runofr  and  groundwater  discnarqe  to  the  lake.  Where 
possible,  the  parameters  required  ror  simulation  are  taKen 
from  real  data.  However,  as  indicated  previously,  some  could 
only  be  estimated  oy  trial  and  error.  These  parameters  are 
discussed  in  detail  in  the  following  pages.  Appendix  III 
contains  a  listing  of  the  Baptiste  Lake  model  input  cara 
deck. 

The  Baptiste  Lake  watershed  was  modelled  as  a  single 
basin  or  subarea.  Neither  the  accuracy  of  the  field  data  nor 
the  structure  of  the  model  justified  an  increased  complexity 
of  watershed  parameters,  for  example,  to  account  for  the 
twelve  streams  flowing  into  the  lake. 

The  Baptiste  Lake  watershed  was  divided  into  eight 
groundwater  storage  elements,  as  illustrated  by  figure  24: 
five  shallow  and  three  deep.  The  elements  correspond  to  the 
hydrostratigraphic  units  discussed  previously.  They  are: 
elements  I  and  2,  units  which  represent  the 
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LEGEND: 

3  ELEMENT  NUMBER  ©  RIVER 

is  PROPORTION  OF  DISCHARGE  (°)  OUT  OF  BASIN 

FROM  ELEMENT  w 

0  LAKE 


Figure  24.  Schematic  representation  of  the 
Baptiste  Lake  watershed  model. 
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alluvium- lacustrine  deposits  with  and  without  an  overlying 
unsaturated  zone;  elements  3  and  4,  which  correspond  to  the 
shallow  till  unit,  again  with  and  without  an  overlying 
unsaturated  zone;  element  5,  an  element  representing  the 
aeolian  sand;  element  6,  the  deep  till  unit;  and  two 
elements,  7  and  8  corresponding  to  the  sand  and  gravel 
deposit  along  the  north  and  south  buried  valleys, 
respectively.  The  percentage  of  the  watershed  surface  that 
is  occupied  by  the  elements  is  listed  on  Table  5. 

Table  5. 

GHCUNDJATSE  STORAGE  ELEMENT  PARAMETERS  ROE  THE 
SIMULATION  OF  THE  3AP1JSTL  LAKE  WAT  EES  HEP 

previous  g.v. 


element  % 

area 

salinity 

DT 

RI 

outflow 

recharg 

1 

silt 

u .  z . 

15 

400 

mg  /L 

1. 

]. 

0. 

cm 

0. 

cm 

2 

silt 

2 

4  00 

mg  /L 

1. 

1  . 

n 

U  « 

cm 

0. 

cm 

3 

till 

u.  z . 

4  6 

500 

mg/L 

1. 

1  . 

0. 

c  m 

0. 

cm 

4 

till 

3  5 

500 

mg  /L 

1. 

1. 

0. 

cm 

0. 

cm 

5 

sand 

u»  z. 

2 

250 

mq/L 

1. 

1. 

0. 

cm 

0. 

cm 

6 

deep 

till 

0 

6  00 

mg/L 

7. 

1  . 

.  005 

cm 

0. 

cm 

7 

N.  channel 

0 

1  1  25 

mg/L 

4, 

1. 

.CGI 

cm 

.001 

cm 

8 

S.  channel 

0 

7  75 

mg/L 

4. 

1. 

.  004 

cm 

.  004 

cm 

The  storage  delay  time  coefficit nts  for  the  groundwater 
storage  elements,  previous  groundwater  outflow  from  the 
storage  elements,  and  previous  recnarge  to  the  elements  were 
obtained  by  trial  and  error.  Model  runs  generally  commenced 
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ir.  the  spring  so  that  the  recharge  to  the  shallow  elements 
and  outflow  during  the  previous  month,  which  was  the  last 
month  of  winter,  might  he  assumed  to  he  zero.  This 
assumption  was  supported  in  all  subsequent  runs  which  showed 
the  outflow  and  recharge  to  be  almost  zero  by  the  end  or  the 
last  month  of  winter.  Initial  estimates  of  the  amount  of 
previous  outflow  from  the  deep  elements  were  estimated  and 
then  reficed  by  comparing  them  to  discharge  calculated 
during  the  latter  years.  Groundwater  storage  delay  time 
coefficients  are  parameters  which  can  only  be  estimated  by 
trial  and  error.  Lower  values  were  initially  used  for  the 
shallow  elements  than  the  deep  elements  because  it  was 
thought  that  the  residence  time  for  groundwater  in  the 
shallow  elements  would  be  less  than  the  deep  elements.  The 
values  of  delay  time,  previous  groundwater  recharge  and 
previous  groundwater  discharge  to  give  the  ’best  fit’ 
simulation  are  listed  on  Table  5. 

Other  watershed  parameters,  sucn  as  groundwater 
discharge  from  an  element,  unsaturated  zone  moisture 
capacity  and  initial  moisture  content,  the  percentage  or  tne 
basin  which  is  impervious,  and  the  groundwater  recharge  rate 
were  first  estimated  from  available  data  then  refined  during 
the  later  runs.  The  amount  of  groundwater  discharge  from 
each  element  that  is  proportioned  to  other  elements  or  to 
the  surface  waters  was  first  estimated  and  finally  arrived 
at  by  trial  and  error  adjustments  to  the  value.  Tne  tinal 
proportions  of  discharge  are  illustrated  on  Figure  2h. 
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Initial  estimates  of  the  upper  ana  lower  unsaturated  zone 
moisture  capacities,  their  present  moisture  content,  the 
maximum  allowable  recharqe  to  the  qroundwater  zone  and  the 
percentage  of  the  basin  that  is  impervious  were  subsequently 
refined,  producing  a  better  fit  to  firstly,  the  fluctations 
of  the  lake  level  and  secondly,  laxe  salinity*  These  values 
are  listed  on  Table  6. 

During  the  calioraticn  process  and  sensitivity  analysis 
(Chapter  III),  it  was  found  that  the  lake  level  fluctuations 
and  lake  salinity  were  more  sensitive  to  some  parameters 
than  to  others.  Small  chanqes  of  the  percentage  of  the  basin 
which  is  impervious,  unsaturated  zone  moisture  capacity  and 
the  maximum  qroundwater  recharqe  rate  resulted  in 
particularly  sensitive  lake  chanqes*  For  example,  a  0.5% 
change  of  tne  impervious  cover,  or  a  0. 5  cm  variation  or  the 
unsaturated  zone  moisture  capacity  and  the  maximum; 
groundwater  recharge  rate  significantly  altered  the  lake 
system.  Parameters  that  were  not  very  sensitive  were  the 
storage  delay  time  coefticient  for  a  qroundwater  storage 
element,  previous  groundwater  outilow  from;  an  element, 
previous  groundwater  recharge  to  an  element,  the  proportion 
of  groundwater  outflow  from  an  element  and  initial  moisture 
storage  in  the  unsaturated  zone,  here,  relatively  large 
changes  could  be  made  without  resulting  m  significant 
changes  to  the  lake  levels  or  chemistry. 

The  salinity  of  the  groundwater  storage  elements  were 
estimated  from  chemical  analyses  studied  in  the  groundwater 
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chemistry  section,  and  are  listed  on  Table  5.  Characteristic 
salinities  of  surface  runoff  and  precipitation  were 
determined  from  analyses  collected  by  the  hater  Quality 
Control  Branch  of  Alnerta  Environment.  These  values  are  also 
presented  on  Table  6. 


Table  6. 

BAPTISTE  LAKE  WATEHSHEfl  ^^dLATICN  PAPA M BILKS 

number  of  subareas:  1 

area  of  the  basin:  299.70  k m2 

percentage  of  impervious  area:  13k 

upper  unsaturated  zone  moisture  capacity:  5.0  cm 
lower  unsaturated  zone  moisture  capacity:  10.0  cm 
previous  upper  unsaturated  zone  moisture:  0.0  cm 
previous  lower  unsaturated  zone  moisture:  0.0  cm 
maximum  groundwater  recharge  rate:  3.0  cm 

salinity  of  precipitation:  10  mq/L 
salinity  of  impervious  flow;  10  mq/L 
salinity  of  int erf low/overiand  flow:  100  mg /L 


Climatic  Variables 

Total  monthly  precipitation  and  average  monthly 
tern perature  were  obtained  from  the  Athabasca  2  meteorolgy 
station  (Canada  Atmospheric  Environment  Services,  1970  - 
1977).  The  station  is  located  approximately  15  km  east  of 
Baptiste  Lake.  Appendix  X  lists  the  precipitation  and 
temperature  measurements. 

There  are  many  problems  assoc iatea  with  trying  to 
measure  evaporation  from  a  watershed,  for  this  reason  many 
empirical  methods,  which  are  based  on  measurable  climatic 
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parameters,  have  been  proposed.  A  simple  but  very  practical 
method  was  developed  by  Thornthwaite  (1948)  and  Ihornthwaite 
and  Mather  (1955) .  The  formula,  based  on  mean  monthly 
temperature  and  length  of  day,  (discussed  1 y  Gray  et  al.. 


19  70)  is: 


Ep  -  1.62 


10  T„ 

4> 


(39) 


where:  Ep  =  monthly  potential  evaporation,  in  cm 
Tm  -  mean  monthly  temperature,  in  °C 
a  =  an  exponent  defined  as: 

a  =  6. 75 x  10  7- 63-  7  71  x  10‘5- 02  +  1.79x  IQ2-  6  +  0.492*  *  *  *  ■ 


6  -  heat  index,  defined  as: 


^  =  daylight  and  latitude  correction  factor 
The  values  of  potential  evaporation  calculated  by  this 
method  are  adjusted  for  changes  in  daylight  hours  and  the 
latitude  of  the  recording  station  by  employing  a  correction 
factor  provided  by  Gray  et  al.  (1970). 

This  method  of  determining  potential  evaporation  has 
been  selected  for  use  at  Eaptiste  Lake  because  of  its 
applicability  to  a  watershed  where  the  detailed  climatic 
data  that  are  required  by  other  evaporation  formula  are  not 
available.  More  importantly,  Laycock  (1967,  1966,  1971, 

1973)  has  shown  the  applicability  of  Thorn thwaite 's 
technique  for  numerous  water- calance  studies  m  the  Canadian 


103 


prairies.  Tests  of  the  various  evaporation  models  have  led 
to  the  conclusion  that  this  formula  produces  the  most 
realistic  water  balance  (Laycock,  1971), 

The  values  of  mean  monthly  temperature  required  in  the 
Thorn thwaite  formula  were  obtained  from  the  Athabasca  2 
meteorolgical  station.  The  calculated  potential  evaporation 
values  are  listed  in  Appendix  X. 

Lake  Functions 

Lake  functions  for  Baptiste  Lake  were  obtained  from  the 
Hydrology  3ranch  of  Alterta  Environment.  These  include  a 
capacity-stage  curve  and  an  area-stage  curve,  drawn  from  a 
hydrographic  survey  of  Baptiste  Lake.  The  curves  and  the 
lake  survey  are  included  in  Appendix  V.  The  functions  that 
define  the  two  curves  were  determined  to  fit  the  data  along 
the  upper  portions  of  the  curves  as  closely  as  possible 
because  this  includes  the  range  of  natural  lake  level 
fluctuations.  The  functions  are: 

Area-Stage  Curve 


(42) 


Capacity-Stage  Curve 


S=  557.82  + 


(43) 


where:  S  =  elevation  of  the  lake  surface. 


in  m  a.m.s.l 
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A  =  area  of  the  lake  surface,  in  km2 
V  =  volume  of  water  m  the  lake,  in  m.3 

A  stage-discharge  curve  for  the  outlet  creek  was  not 
available.  Surface  water  outflow  from  the  lake  was 
calculated  according  to  a  method  oiscussea  by  kigham  (1970): 

Q=dL.eJl/KL  . (44) 

°°  *L 

where: Q0  =  surface  water  outflow,  in  m 3 

V0  -  volume  of  the  lake  above  the  surface 
water  outlet  level,  in  m3 
Tl  =  time  of  uniform  recharge  to  the  lake, 
in  months 

Kl =  storage  delay  time  of  the  lake, 
in  months 

The  lake  surface  elevations  were  recorded  by  the  fcater 
Survey  of  Canada.  All  initial  parameters  describing  the  lake 
are  listed  in  Table  7. 


Table  7. 

BAPTISTE  LAO  SIMULATION  PARAMETERS 

surface  outflow  elevation:  577. 76  m 
initial  lake  elevation:  576.60  m 
initial  lake  salinity:  200  mg/L 
lake  recharge  time:  0.5  months 
lake  storage  delay  time:  1.0  month 
lake  discharge  via  groundwater:  4,0  cm 
initial  ice  thickness:  100  cm 
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The  Simulation  of  Baptiste  Lake 

A  satisfactory  simulation  of  lake  surface  levels  and 
lake  salinity  fluctuations  at  Baptiste  Lake  was  achieved 
with  the  model.  A  comparison  cf  the  actual  field 
measurements  of  lake  level  and  salinity  with  the  model 
calculations  are  illustrated  on  Fiqure  25  and  figure  26 
respective ly . 

Although  the  actual  and  calculated  values  are  not 
always  the  same,  the  patterns  of  fluctuation  follow  the  same 
trends.  The  discrepancies  between  the  two  sets  of  results  is 
thought  to  be  caused  by:  (l)  model  error,  (ii)  data  error, 
and,  (lii)  inaccuracy  of  hydrologic  functions. 

Model  error  is  probably  the  most  important  and  is 
attributed  to  the  simplicity  cf  the  model.  For  example,  the 
geology  of  the  watershed  has  been  partitioned  into  only 
eight  hyarcstratigraphic  units.  Variability  in  seme  surface 
features,  such  as  vegetation,  have  not  been  considered,  and 
it  is  assumed  that  spring  runoff  commences  at  the  same  time 
each  year.  Also,  the  lake  is  not  considered  to  be 
stratified,  precipitation  of  salts  in  the  lake  is  assumed 
not  to  occur,  and  groundwater  salinity  assigned  to  the 
storage  elements  is  an  averaged  general  value. 

Data  errors  result  from  not  having  a  me teoroig ical 
station  withir  the  watershed.  Precipitation,  potential 
evaporation  and  temperature  are  measured  15  km  away  at 
Athabasca.  In  addition,  there  are  no  appropriate  ways  of 
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determining  whether  the  averaged  groundwater  chemistry 
values  are  a  reasonable  characterization  of  the  complex 
natural  situation.  However,  the  sensitivity  analyses 
indicate  that  these  data  errors  should  not  be  a  significant 
source  of  error  in  the  simulation. 

The  last  group  of  errors,  the  inaccuracy  in  the 
hydrologic  functions,  also  cause  a  significant  deviation  of 
the  model  calculations  from  actual  values.  The  surface  water 
outflow  probably  causes  the  sinqle  largest  source  of  error. 

A  stage- d ischarge  curve  was  net  available  for  the  outlet 
creek.  However,  the  function  used  to  estimate  surface  water 
outflow  was  thought  to  provide  a  reasonable  estimate  because 
it  has  been  shown  by  Wigham  (1970)  to  be  an  acceptable 
formula  for  calculating  discharge  from  a  lake.  Other 
problems  will  arise  in  the  model  study  from  applying  an 
inaccurate  function  to  the  stage- area  and  stage- capacity 
curve.  While  it  is  difficult  tc  assess  the  accuracy  of  the 
lake  survey,  it  can  only  be  considered  approximate. 

Differences  between  real  and  simulated  lake  levels 
could  occur  also  as  a  result  of  changes  to  the  outlet  creek, 
resulting  in  increased  or  decreased  outflow  from  the  lake. 
Fcr  example,  the  outlet  creek  could  have  been  blocked  by 
beaver  dams  or  slumping  of  its  banks  (Water  Eesources 
Mangement  Division,  1979)  causing  higher  than  predicted  lake 
levels  during  the  latter  part  of  1977.  A  bridge  replaced  the 
existing  culverts  during  the  fall  of  1970  and  the  stream 
embarkment  was  removed  in  the  summer  of  1971  (Water 


-Ml}  U'i  v*>i.iirSp6-ni  *  *  ,«•*>**  ic  •> ' {t : 


109 


Resources  Management  Division,  1979).  Thus,  possibly 
explaining  why  the  model  predicted  higher  lake  levels  than 
were  actually  recorded  during  1972. 

The  simulation  results  point  to  the  importance  of 
groundwater  in  the  hydrologic  balance  ot  Baptiste  Lake. 

Table  8  lists  the  percentages  of  the  various  inflows  to,  and 
outflows  from  Baptiste  Lake,  for  the  period  1972  to  1977, 
calculated  with  the  model.  During  this  time,  the  groundwater 
contributed  an  average  of  13X  of  ail  recharge  to  the  lake. 

1  a  1 1  e  8 . 

PROPORTIONS  C±  THE  HID RO LOGIC  COM  ROh  £K IS  AT  BAPTIST  £  TAKE 


g.  w. 

%  inflow: 

s.  w . 

prec. 

% 

Cj  •  Vi  4 

outflow 

S.  V. 

-• 

evap 

19  72 

22 

68 

10 

13 

72 

1  5 

19  73 

1  3 

68 

19 

1  4 

70 

1  6 

197a 

18 

74 

8 

8 

83 

9 

19  75 

9 

74 

1  7 

14 

6  9 

1  7 

1976 

7 

76 

17 

1  8 

6  0 

22 

19  7  7 

9 

75 

16 

13 

72 

15 

The  reader  is  cautioned  that  Table  8  and  the  given 
value  of  13:4  for  tne  average  groundwater  recharge  to  the 
lake  are  only  an  estimates.  These  values  are  based  on  over 
50  trial  runs;  not  just  one  simulation.  And,  while  these 
numbers  arc  as  accurate  as  the  model  and  the  data  permits, 
the  values  are  calculated  only  to  provide  the  reader  with 
insight  into  the  significance  of  groundwater  recharge  to 
Baptiste  Lake. 

Tc  further  verify  this  result,  the  watershed  was 
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modelled  assuming  no  groundwater  discharge  to  the  lake  and 
no  discharge  from  the  lake  via  groundwater.  Ail  water  was 
assumed  to  enter  the  lake  as  eitaer  surface  runoff  or  as 
precipitation  falling  directly  on  the  surface.  The  resulting 
lake  level  f luctuations,  Figure  27,  are  similar  to  the 
actual  values.  However,  a  closer  fit  can  be  obtained  when 
groundwater  component  of  the  watershed  budget  is  included . 
Lake  salinity,  without  groundwater  input,  is  less  than  half 
that  as  when  groundwater  discharge  is  included  in  the 
simulation  (Figure  23).  Also,  the  salinity  fluctuations  mat 
exist  at  Eaptiste  Lake  are  not  shown,  lhese  results  present 
a  striking  illustration  of  the  importance  cf  the  groundwater 
contribution  to  Baptiste  Lake. 

Further  Study  at  Baptiste  Lake 

In  the  previous  two  sections,  a  simple  hydrologic  model 
simulated  lake  level  and  salinity  fluctuations  at  Eaptiste 
Lake.  Now,  a  preliminary  extension  is  undertaken  to  evaluate 
one  possible  impact  of  watershed  development  at  Baptiste 
Lake. 

This  study  is  concerned  with  the  possible  impact  that  a 
new  cottage  development,  situated  on  the  shore  of  Baptiste 
Lake,  will  have  on  the  lake.  The  cottaqe  development  was 
simulated  by  including  a  shallow  groundwater  element  which 
covers  approximately  5*  (15  km2),  of  the  surface  of  the 
watershed,  having  hydrogeological  characteristics  of  the 
shallow  till  storage  element  with  no  overlying  unsaturated 
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zone,  (element  #4).  A  total  dissolved  solids  content  of  2000 
mg/L,  approximately  10  times  that  of  the  lake,  is  designed 
to  reflect  an  increased  chemical  parameter  over  background 
concentrations.  While  salinity  may  not  be  the  test  chemical 
parameter  to  simulate  cottage  development,  it  is  used  nere 
for  two  reasons.  Firstly,  salinity  is  the  only  chemical 
parameter  routed  through  this  model  and,  secondly,  the 
purpose  of  this  simulation  is  only  to  illustrate  in  a 
general  way  the  effect  that  any  contaminant  may  have  on  lake 
chemistry.  All  groundwater  discharging  from  the  element 
flows  directly  to  the  lake.  Other  element  parameters  remain 
the  same  as  defined  for  element  #4. 

The  resultant  change  in  the  laxe  water  quality. 

Figure  28,  is  that  the  salinity  of  the  lake  increased  ty 
approximately  50  mg/L,  thus  indicating  a  significant  rise  in 
lake  salinity  over  background  conditions.  This  simulation 
assumed  the  worst  hydrologic  conditions,  such  as  no 
overlying  unsaturated  zone  and  all  water  is  discharged 
directly  tc  the  lake.  With  no  unsaturated  zone,  the  water 
table  is  at  the  ground  surface  and  all  contaminants  will 
immediately  enter  the  groundwater  system. 

It  can  probably  be  assumed  that  a  development  would  be 
situated  in  an  area  where  the  water  table  is  well  below  the 
ground  surface,  presenting  the  possibility  ot  contaminant 
attenuation  in  the  unsaturated  zone.  Another  analysis  was 
conducted  with  the  development  located  on  the  same  till 
unit,  but  included  an  overlying  unsaturated  zone.  This 
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resulted  in  restricted  infiltration  to  the  groundwater  zone 
and  hence,  a  negligible  influence  on  lake  water  quality. 

However,  dealing  with  total  dissolved  solids  as  a 
parameter  to  study  lake  contairmatioc  from  cottage 
developments,  may  be  misleading.  It  should  be  remembered 
that  this  analyses  shows  tnat  a  single  parameter,  in  hign 
concentrations  in  the  effluent  above  background  lake  values 
could  have  a  significant  influence  on  the  chemistry  of  the 
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V.  CONCLUSIONS 


The  model  developed  in  this  thesis  provides  a  useful 
means  of  simulating  the  hydrological  and  hydrochemical 
processes  of  a  lake- watershed  system.  The  model  is  simple 
and  is  designed  to  he  used  when  only  an  estimation  of  the 
various  hydrologic  components  of  the  watershed  is  required. 
Only  a  limited  amount  of  basic  physical,  climatic, 
hydrological  and  hydrochemical  data  is  required  by  the  model 
for  hydrological  simulation  of  a  lake- watershed  system. 

It  is  possible  to  use  this  kind  of  model  to  learn  more 
about  the  processes  and  parameters  controlling  the  physical 
and  chemical  hydrology  of  lakes  in  a  prairie  setting. 
Sensitivity  analyses  of  actual  or  hypothetical  cases  can 
lead  to  a  better  understanding  of  how  the  system  functions. 
Although  caution  is  warranted,  generalizat le  results  for 
sets  of  lake  systems  should  be  possible  from  the  analyses. 
Sensitivity  analyses  of  hypothetical  lake- wa te rsh ed  systems 
showed: 

1.  If  it  were  physically  possible  to  alter  a  la*e  by 

concen  tratirig  the  same  volume  of  water  into  a  lake  wita 
a  smaller  surface  area,  the  salinity  of  the  lake  should 
decrease  and  lake  levels  rise.  The  reason  for  this 
occurrence  is  that  evaporation  exhibits  less  of  a 
influence  on  a  lake  with  a  small  surface  area  and  a 
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greater  influence  ever  lakes  with  a  large  surface  area. 
Tne  analysis  of  the  hypothetical  lake- watershed  showed 
that  if  the  lake  volume  was  increased,  without  a 
corresponding  enlatqenient  of  its  surface  area,  the  same 
results  as  above  are  obtained. 

The  study  of  the  formation  or  ice,  an  important  process 
on  prairie  lakes,  indicated  that  the  more  water  that  is 
taken  up  from  tne  lake  during  tne  rormation  of  ice,  the 
greater  the  salinity  cf  the  lake  (the  salinity 
fluctuations  in  shallow  lakes  are  greater  than  in  deep 
lakes,  both  containing  equal  volumes  of  water). 

Small  variations  ic  groundwater  chemistry  or  large 
variations  over  a  part  of  tne  watershed  did  not 
significantly  alter  the  salinity  of  the  hypothetical 
prairie  lake.  Large  increases  or  decreases  in  the 
quality  of  groundwater  over  the  entire  watershed  should, 
however,  alter  the  salinity  of  the  lake. 

The  unsaturated  zone,  primarily  controls  the  amount  of 
water  percolating  to  the  groundwater  zone  and  hence, 
indirectly  influences  the  amount  of  groundwater  that  is 
discharged  to  a  lake.  Decreasing  the  unsaturated  zor.K 
moisture  capacity  in  the  test  watershed  caused  increased 
surface  runoff  and  groundwater  discharge,  which  resulted 
in  higher  lake  levels  and  salinity. 

As  the  proportion  of  a  watershed  that  is  covered  ty  an 
unsaturated  zone  decreases,  the  amount  of  groundwater 
recharge  increases.  Hence,  the  volume  of  groundwater 


that  discharges  to  the  lake  is  increased,  and  the 
quality  of  tne  lake  water  tends  towards  that  of  the 
groundwater, 

7.  Increasing  the  area  of  this  hypothetical  watershed  that 
is  impervious,  for  example  by  urbanization,  causes 
precipitation  to  flow  directly  to  the  lake,  decreasing 
the  lake  salinity  which  increases  fluctuations  of  lake 
levels,  and  generally  increases  the  lake  stage. 

8.  Variations  m  precipitation  and  potential  evaporation 
did  not  cause  large  fluctuations  in  lake  levels  and 
salinity  in  the  large  watershed  used  in  the  sensitivity 
analysis,  but  did  cause  significant  deviations  in  the 
small  lake- watershed .  The  division  between  large  and 
small  watersheds  is  estimated  to  be  near  100  km2.  In 
particular,  under-estimating  precipitation  or 
over-estimating  potential  evaporation  caused  a  much 
larger  variation  than  if  fceth  parameters  were  either 
over-estimated  or  under-estimated.  The  major 
generalization  which  nay  be  drawn  from  this  study  is 
that  precipitation  and  potential  evaporation  exhibit  a 
greater  influence  on  lake  levels  and  salinity  in  a  small 
watershed  than  groundwater  discharge.  The  opposite  is 
true  in  a  large  lake-watershed. 

Modelling  tne  Baptiste  Lake  watershed  has  been  a  userui 
aid  fer  studying  hydrologic  processes  in  this  basin.  The 
model  has  only  been  applied  to  this  one  watershed  but  its 
implications  are  that  it  can  successfully  be  used  for 
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studying  the  water  balance  of  many  prairie  lake- watershed . 

In  the  case  or  Baptiste  Lake,  the  model  was  able  to 
simulate,  with  a  reasonable  degree  o£  accuracy,  a  two  year 
record  of  water  levels  and  chemistry.  .Examination  of  the 
relative  importance  of  various  components  of  the  hydrologic 
cycle  of  the  lake  showed  that  groundwater  inflow  probably 
constitutes  a  significant  proportion  cf  all  recharge  to  the 
lake,  which  was  calculated  to  be  an  average  or  13$.  A  model 
of  a  lake-watershed  system  can  be  used  to  study  various 
processes  and  the  impact  of  changes  to  the  watershed. 

A  further  capability  of  this  model  is  to  extend  a 
simulated  lake  record,  such  as  was  done  at  Baptiste  Lake,  to 
predict  future  changes  in  the  watershed.  This  use  of  the 
model  is  particularly  useful  for  application  to  watershed 
development,  changes  in  climatic  patterns  or  examining  the 
past  natural  history  of  the  lake. 

This  model  study  has  only  examined  a  few  of  the 
numerous  natural  and  man-made  processes  that  govern  lake 
water  cbenistry  and  level  fluctuations.  The  results  obtained 
in  this  thesis  indicate  that  further  study  of  watersheds 
with  this  model  rs  justified  and  development  or  the  model  rs 
warranted.  Such  developments  could  include:  the  expansion  of 
the  model  to  study  watershed  parameters  in  more  detail,  the 
development  of  functions  to  stratify  lakes  ana  to  allow  for 
precipitation  of  salts  from  lake  water,  and  accounting  for 
chemical  reactrons  involving  individual  nutrients  as  they 
are  routed  through  watersheds  to  lakes. 


EEFEBENCES 


1  19 


Ayers,  Ii.  D.  ;  197  3;  Basin  Yield;  Section  X  of  Handbook  on  tne 

Principles  of  Hydrology;  Editor-in-Chief:  D. a.  Gray; 
fcater  Information  Centre  Inc.;  pp,  10.1-10.20. 

Bcrneuf,  D. ;  1973;  Hydrogeology  of  the  Tawatinaw  Area, 

Alberta,  NTS  831;  Alberta  Research  Council  Report  72-11. 

Carlson,  V.  A.  ;  1  9  76;  bedrock  Topography  of  the  Tawatinaw  dap 
Area,  NTS  831,  Alberta  ;  Alberta  Research  Council  dap. 

Canada  Atmospheric  Environment  Services;  1970-1977;  Monthly 
Records,  Meteorological  Observations  in  Canada; 
fisheries  and  Environment  Canaaa. 

Crawford,  N. H. ,  and  B.  K,  Linsley;  1966;  Digital  Simulation 
in  hydrology;  Stanford  watershed  Model  IV;  Dept,  of 
Civil  Eng.,  Stanford  Univ.,  Tech.  Rep.  No.  39. 

Crawford,  N.H.,  and  A. S.  Donigian;  1973;  Pesticide  Transport 
and  Runoff  Model  for  Agricultural  Lands;  Grfice  of 
Reasearch  and  Development,  U.S.  Envir,  Protection 
Agency,  Rep.  No.  E EA-66 0/2-74-0 13. 

Crowe,  A.;  1978;  Baptiste  Lake  Groundwater  Investigation; 

Preliminary  Report;  Alberta  Environment,  Environmental 
Protection  Services,  Earth  Sciences  Division  Report. 

Dooge,  J.C.I.;  I960;  The  Routing  of  Groundwater  Recharge 

through  Typical  Elements  of  Linear  Storage;  Int.  Assoc, 
of  Sci.  Hydrology  General  Assemtly  of  Helsinki  I.U.G.G.; 
Publ.  No.  52,  pp.  286-300. 


120 


EPEC  Consulting  Ltd.;  1976;  Cookinq  Lake  Area  Study,  Vol,  I: 
Formulation  of  Management  Alternatives;  Feport  prepared 
for  the  Cooking  Lake  Area  Stuay  Management  Committee, 
Planning  Division,  Alberta  Environment. 

Feniak,  M . ;  1944;  Athabasca  -  Barrhead  Map  Area,  Alberta; 

C.S.C.  Taper  44-6;  Can.  Dept  or  Mines  ana  Resources, 
Mines  and  Geology  Branch. 

Fertuck,  L.J.,  J.  W.  Spyker  and  W . H . W .  Hus banc;  1970;  A 

Numerical  Method  for  Estiaating  Ice  Growth  as  a  Function 
of  Air  Temperature,  Wind  Speed  and  Snow  Cover;  Proc. 
Thrd.  Western  Can.  Heat  Transfer  Com.,  Edmonton, 
Alcerta,  May  28-29,  1970. 

Fields,  D. E.  and  S. B.  Watson;  1975;  OPTEM  -  A  Hydrologic 

Transport  Model  witn  Parameter  Optimization;  Oak  Biage 
National  Laboratory,  Eep.  OE NL— NSF-EATC- 1 4. 

Fritz,  P.  and  H.R.  Krouse;  1973;  Wabamun  Lake  Past  and 

Present:  An  Isotope  Study  of  tne  Water  Budget;  Proc.  of 
the  Symp.  on  the  Lakes  of  Western  Canada,  1973, 

Edmonton,  Canada;  Pub.  No.  2,  Univ.  of  Alberta  Water 
Resources  Centre.,  pp.  244-25S 

Gray,  D.M.,  G. A.  McKay  and  J. M.  Wigham;  1970;  Energy, 
Evaporation  and  Evapo transpiatior.  ;  Section  III  or 
Handbook  cn  the  Principles  of  hydrology; 

Editor-in-Chief:  D.M.  Gray;  Water  Information  Centre 
Inc. ,  pp .  3 . 54-3. 6  1 . 


Green,  fi.  ;  1  9  72  ;  Geological  Map  of  Alberta;  Alberta  Reaseach 

Council  Map  No.  35. 

Gull  Lake  lask  Force;  1970;  Gull  Lake  Study,  Vol  1;  Report 
prepared  for  the  Gull  Lake  Co-ordinating  Committee, 

Water  Resources  Division,  Alberta  Dept,  of  Agriculture. 

Huff,  D.O.;  1968;  Simulation  of  the  hydrologic  Lransport  of 

Radioactive  Aeresols;  Ph,  C.  thesis,  Stanford  Univ. 

Ingersoll,  L. E. ,  G.J.  Zobel  and  A. C.  Inqersoll;  1948;  Heat 
Conduction,  with  .Engineering  and  Geological 
Applications;  McGrav  -  Hill  publishers. 

Kjearsgaard,  A. A.;  1972;  Soil  Survey  of  the  Tawatinaw  Map 

Sheet  (831) ;  Alberta  Soils  Survey  No,  29. 

Laycock,  A.H.;  1967;  Water  Defeciency  and  Surplus  Patterns 

in  the  Prairie  Provinces;  Prairie  Provinces  Water  Eoara, 
Prairie  Farm  Behabilita t ion  Administration  Hep.  No,  13. 

- ;  1968;  Water  Balance  of  the  Gull  Lake  Basin,  A  System 

Analysis  for  the  Period  1921-1968;  Report  for  the  Red 
Deer  Regional  Planning  Commission. 

- ;  1973;  Lake  Level  Fluctuations  and  Climatic  Variation 

in  Central  Alberta;  Proc.  of  the  Syrcp.  on  the  Lakes  of 
Western  Canada,  1973,  Edmonton,  Canada;  Pub.  No.  2, 

Oniv.  of  Alberta  Water  Resources  Centre.,  pp.  38-96. 

- ;  1971;  Water  Balance  in  the  Cooking  Lake  Region; 

Presentation  to  the  Public  Hearings  of  the  Environmental 
Conservation  Authority  on  Cooking  Lake,  Edmonton. 


'  .  j  ’ .  !  •  «  »  4  * 


122 


Linsey,  R.K.,  M. A.  Kohler  and  J .  L . H .  Paulhus;  1975; 

Hydrology  for  Engineers;  2nd  Edition,  McGraw  -  Hiii 
publishers. 

Nursall,  J.R.,  J.B.Nuttall  and  ?.  Fritz;  1971;  The  Effect  of 
Thermal  Effluent  on  Lake  Wabamun;  XVI 11  Limnological 
Congress,  Leningrad,  1971. 

Patterson,  M.R.  et  al. ;  1974;  A  User’s  Manual  ror  the 

FORTRAN  IV  Version  of  the  Wisconsin  hydrologic  Transport 
Model;  Oak  Ridge  National  Laboratory,  Rep. 
OENL-NSF-EATC-7. 

Penman,  H.  I.  ;  1948;  Natural  Evaporation  from  Open  Water, 

Eare  Soil  and  Grass;  Eroc.  Royal  Society  of  London, 
Series  193a,  pp.  120-145. 

Planning  Division,  Alberta  Environment;  1979;  Sylvan  Lake 
E^egulat.icn  Study;  Report  prepared  for  the  Red  Deer 
Regional  Planning  Commission. 

Pounder,  E.R.;  1965;  Physics  of  Ice;  Perga ton  Press. 

Prairie  Provinces  Water  Board;  1952;  Evaporation  froc  Lakes 
and  Reservoirs  on  the  Canadian  Prairies:  A  Study  Based 
on  30  Years  of  Meteorological  Records,  1921-1950;  Report 
#5. 

Saska tche wan-Nelson  Basin  Board;  1972;  Water  Supply  for  the 
Saskatchewan-Nelson  Basin,  Appendix  4:  Hydrology, 


Voi.  1. 


■ 


123 


Schwartz,  F.il.  and  D.N,  Gallup;  1978;  Some  Factors 

Coutollinq  the  Major  Ion  Chemistry  Gf  Small  Lakes: 
Examples  From  the  Prairie  Parkland  of  Canada; 
Hydrobiolgoia,  Vol.  56,  No.  1,  pp.  65-81. 

Schwartz,  f.  W.  ;  1979  ;  Interim  Report  or.  a  Hydrogeological 

Investigation  of  the  Muskeg  River  Basin,  Alberta;  Report 
prepared  for  the  Alberta  Cil  Sands  Environmental 
Research  Program;  Project  hG  I.  i,  ACSEKF  Report  68. 

Stanley  and  Associates;  1974;  Cooking  Lake  Study,  Vol.  2: 
"water  Inventory  and  Demands;  Report  prepared  for  the 
Cooking  Lake  Area  Study  Management  Committee,  Planning 
Division,  Alberta  Environment. 

Thornthwai te ,  C. K. ;  1948;  An  Approach  Towards  a  Rational 

Classification  of  Climate;  Geographical  Review,  Vcl. 

38,,  No.  1.,  pp.  55-94, 

Thorn th waite ,  C. H.  and  J. I.  Mather;  1955;  The  water  Balance; 
Laboratory  of  Climatology,  Publ,  in  Climatology,  Vol.  8, 
No,  1. 

Trew,  D.O.,  E.J.  Beliveau  and  E.I.  i'onge;  1979  ;  Baptiste 

Lake  Study  Summary  Report;  Water  Quality  Control  Branch, 
Pollution  Control  Division,  Alberta  Environment. 

Ward,  R.C.;  1975;  Principles  cf  Hydrology;  McGraw  -  Hill 

publishers. 

Water  Quality  Control  Eranch;  1977;  Baptiste  Lake  Study 
Proposal;  Pollution  Control  Division,  Alberta 


invironment. 


124 


Water  Hesoucces  Management  Division;  1979;  Baptiste  Lake 
66-24-W4,  correspondence  file:  file  #13011;  Alberta 
Envircnioen  t. 

Wigham,  J.  M.  ;  19  70;  PeaK  Flow  -  Rainfall  Events;  Section 

VIII  of  Handbook  on  the  Principles  of  Hydrology; 
Editor-in-Chief:  D.M,  dray;  Water  Information  Centre 
Inc.  ,  pp.  6.  4  7-8.  64. 


- 


12  5 


APPENDIX  I 
MODEL  USEE’S  GUIDE 


126 


USEE’S  GUIDE 


bifiensio  n  1  ng  Arrays 


There  are  several  subscripted  variables  used  in  the 
program  and  these  must  be  dimensioned  to  the  correct  array 
size,  as  indicated  in  the  following: 

HEAL  LKSTG, LKVOL , ICE, ICEF, SUB  {  10) ,T1ILE(20)  , 

r  (nr)  ,  P  E  (nr)  ,  AE(nr)  ,T£i1P  (nr)  ,  1  NE  ILT  (ns)  ,CIHP  (ns)  , 

CSRO  (ns) ,INTfi  (ns)  ,PCTI M?  (ns)  , AES DB  (ns) ,RGfc (ns)  , 

EGWI1AX  (ns)  ,  QGWLKT  (ns)  ,  GSZMP  (ns)  ,  CGfcLK  (ns)  ,0SZK  (ns), 
USZMAX  (ns)  ,  L5ZM?  (ns)  ,LSZ»  (ns) , LSZtTAX (ns) ,CR¥ (ns)  , 
yINTR  (ns) ,QBVT (ns) , QIM? (ns)  ,QG «RV  (ns,ne) ,HELM (ns, ne)  , 
BEL  M  P  (ns,ne)  ,PCTAB(ns,ne) , QGw (ns,ne) ,  QGfeP  (ns,ne) , 

S (ns , ne)  ,2? [ns , ne) , CGw  (ns , ne) ,  £Gk  LK (ns , ne)  , DT (ns, ne)  , 
ET  (ns, ne) ,C0  (ns,ne)  ,C)  (ns,ne)  , C2  (ns ,ne) ,C3 (cs, ne)  , 

CL  (ns,ne)  ,C  5  (ns  ,ne)  , ECTELN(ns,ne,6) 

IN TEG EE *2  NOPP (ns , ne) , TOE LS  (ns , ne, 6 ) , NGW E ( ns) 

LOGICAL* 1  NEGEEC (ns,ne) ,REPREC  (ns,ne) , NEGQGW (ns,ne)  , 

ELENAM  (7) 

where:  ns  =  number  of  subareas 
nr  =  number  of  records 

ne  =  number  of  groundwater  storage  elements 
in  each  subarea 


User  Supplied  Functions 


Three  functions  are  reguired  to  be  included  at  the  end 
of  the  program  by  the  model  user.  These  are: 


A  E  E  A  L  K  ( 1  k  s  t  g ) 


ELEVLK  (  lkvcl) 


OUTELO  (lkstg) 


VOLLK  (lkstg) 


A  function  used  to  calculate  the  surface 
area  of  a  lake  in  km2,  given  a  lake 
stage  in  m  above  a  datum. 

A  function  which  calculates  the  stage  of 
a  lake  in  m  above  a  datum,  given  the 
volume  of  water  in  the  lake  in  m3» 

A  function  which  calculates  the  surface 
water  outflow  from  a  lake  in  m3, 
given  the  lake  stage  in  m  above  a  datum. 
This  function  finds  the  volume  of  water 
in  the  lake,  in  m3,  given  the  lake 
stage  in  m  above  a  datum. 


The  calling  parameters  for  the  subroutines  are: 

lkstg  =  elevation  of  the  lake,  in  m  above  a  datum 
lkvcl  =  the  volume  of  water  in  the  lake,  in  m3 


- 


. 

Input  Variables 

initial  program  parameters: 


TITLE 

NR  EC 
NSUE 

:title  of  data  set 
: number  of  records 
:number  of  subareas 

subarea  parameters: 


AHSUB 

PCTIMP 

:area  of  a  subarea  (km) 

: fraction  of  subarea  which  is  impervious 
(C  <  PCI  IMP  <  1) 

NGWE 

:numfcer  of  groundwater  storage  elements  ir. 
a  subaraea 

FGWMAX 

:maximum  allowable  recnarge  from  the  unsaturated 
zone  to  the  groundwater  zone,  in  a  subarea  (cm) 

initial  lake  parameters: 


G  w  0  L 

LKSTG 

:surface  water  cutlet  elevation  (m  above  a  datum) 
:initiai  elevation  cf  the  lake  surface  (m  above 
a  datum) 

ICE 

SNOW 

DTELh 

EIFLK 

TEMPP 

:  initial  ice  thickness  on  the  lake  surface  (cm.) 

: initial  depth  pf  sncw  cover  on  the  lake  (cm) 
:storage  delay  time  for  lake  (months) 

:duration  of  recharge  to  lake  (months) 

: average  monthly  air  temperature  of  the  previous 
month  (°  C ) 

unsaturated  zone  parameters: 


USZ  MAX 

:maximum  alicwatie  moisture  storage  in  the  upper 
unsaturated  zone  (cm) 

USZM? 

: present  moisture  content  of  the  upper 
u n saturated  zone  (cm) 

LSZMAX  : maximum  allowable  moisture  storage  in  the  lower 


LSZMP 

unsaturated  zone  (cm) 

: present  moisture  content  of  the  lower 
un saturated  zone  (cm) 

groundwater  element  parameters: 


ELENAM 

:7  character  identifier  for  a  groundwater 
storage  element 

DT 

FT 

KELMP 

: storage  delay  time  (months) 

:time  or  recnarge  to  an  element  (months) 

:recharge  to  an  element  during  the  previous 

QGK? 

month  (cm) 

:groundwater  discharge  from  an  element  during 
the  previous  month  (cm) 

REPREC  :check  to  see  if  recharge  from  precipitation  or 
from  the  unsaturated  zone  is  allowed  (1  or  F) 
NEGREC  :test  to  determine  if  evaporation  may  occur  from 


NEGQGW 

an  element  (T  or  F) 

:check  to  determine  whether  recharge  can  originate 
from  a  surface  water  body  (1  or  F) 

. 
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PCT  A  K 

NOPP 
PCT  ELM 

l'CELM 


ifraction  of  outrlow  from  aii  element  that  is 
routed  to  other  elements  or  the  surface 
<0  <  PCT AR  <  1) 

: number  of  divisions  of  an  element  (0  <  NOPP  <  1) 

:portion  of  discharge  from  an  element  that  is 
routed  according  to  NOPP  (0  5  PCTiLM  <  1) 

:element  number  or  surface  water  body  to  which 
groundwater  outflow  is  routed 


chemistry 

CP 

CLK 

ciap 

CSRO 

CGW 


pa  raiaeters : 

: water  guality  of  precipitation  (mq/L) 

: initial  salinity  oi  the  lake  (mg/L) 

:salinity  of  flew  from  impervious  areas  (mg/L) 
:salmity  of  surface  runoff  (mq/L) 

;salinity  of  a  groundwater  storage  element  (mg/L) 


climatic  variables: 

TEMP  : average  monthly  air  temperature  (°C) 

PE  : total  montnly  potential  evaporation  (cm) 

P  :tctal  monthly  precipitation  (cm) 


Input  Card  Instructions 


The  variables  input  into  the  program  must  conform 
to  the  following  cards,  order,  format  and  columns: 


CARD  # 

VARB. # 

VARIABLE 

FORMAT 

COLUMNS 

1 

1 

TITLE 

20A4 

1 

-  60 

2 

1 

NREC 

15 

1 

—  c; 

2 

NS  UB 

15 

6 

-  1C 

3 

CP 

F  1  0  a  0 

1  1 

-  2  0 

4 

TEMPI? 

i  5.  0 

2  1 

-  25 

5 

ICE 

15.  0 

26 

-  30 

6 

SNOW 

F 5.  0 

3  1 

-  35 

3 

1 

LKSTG 

I  10.  0 

1 

-  10 

2 

SWC-L 

i  1  0 . 0 

1  1 

-  20 

3 

CLK 

i  10.  0 

2  1 

-  30 

4 

DTE  LK 

i  10.  0 

3  1 

-  40 

5 

R1FLK 

i  10.  0 

4  1 

-  50 

6 

GWOLK 

f  1  0.  0 

5  1 

-  60 

■ 
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CARD  ft 

V  ALB.  # 

V  ARI ABLE 

FOE MAT 

COLUMNS 

41 

1 

ABSUB 

110.0 

1 

-  10 

2 

RGKMAX 

F  1  0.  0 

1  ! 

-  20 

3 

C.IMP 

i  10.  0 

2  1 

-  30 

4 

CSHO 

1  10.  0 

3  1 

-  40 

5 

PCTIMP 

F  5 .  0 

4  1 

-  45 

6 

USZMAX 

F  5.  0 

4  6 

-  50 

7 

LS2MAX 

1  5.  0 

5  1 

-  55 

8 

USZMP 

F5.  0 

56 

-  60 

9 

LSZHP 

F5 . 0 

6  1 

-  65 

10 

NG  RE 

15 

6  6 

-  70 

52 

1 

ELEN  AM 

7  A  1 

1 

-  7 

2 

REPBEC 

L  1 

8 

3 

NEGREC 

L  1 

9 

4 

NEGQGW 

L  1 

10 

5 

DT 

r  5.  0 

1  1 

-  15 

6 

RT 

F 5.  0 

16 

-  20 

7 

CGKP 

F5. 0 

2  1 

-  25 

8 

SELMP 

F 5.  0 

2  6 

-  30 

9 

PCT  AE 

F 5.  0 

3  1 

-  35 

10 

C  Gw 

F  10.  0 

3  6 

-  4  5 

1  1 

NO  PP 

15 

46 

-  50 

6  3 

! 

PCI EL  M 

6  F  5  .  0 

1  - 

5,  6  - 

10 

2 

TOELM 

615 

31  - 

35,  36  - 

40 

7* 

1 

P 

ILF  5. 0 

1  - 

5,  6  - 

10 

2 

PE 

12F5.  0 

1  - 

5 ,  6  - 

10 

3 

TEMP 

12F5.  0 

1  - 

5,  6  - 

10 

notes:  1 :  1  card  f o  t 

each 

subarea , 

followed 

by  ; 

2:  1  card  for 

each 

element , 

followed 

by; 

3:  all  6  PC  T  EL  M  and  TOELM  values  on  1  card 
*:  12  records  pec  card,  followed  Ly  icre  of 

the  same  record  if  required 


. 
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APPENDIX  III 
SAMPLE  INPUT  CARL  DECK 


LO 


O  •  O^-nr^fOrHhOO  OHDrOrHrHCOCO'O 

oo  CNJHDrocoo-OLOooooooooo . 

O  . r'-LOOJrHCOLOtOLO 

•  c\jroorOrHC\jcjoJOOoooooo  i  i  i  * — i  i  i  i  i 

O'  I 

o  r-n  rHOno  O' ior>o  OvONO'NNKiO'o 

•  O'Mi^sfCCO'XlOOOOOOOOO . 

o  . .  •^NOOvfO'OL'i 

rH  OinoJlOHHOfJOOOOOOOO  I  HHKH  I  I  I 

I  I  I  I 

in  -Kl  fO  n  Kl  CO  O'f'^rHrHO'OJHDO  O  CJ  LO  O' rH  HD  100*0 

♦LO  COMflvDvDsrcOOOOOOOOOO . 


o 


K1  •rr 1  .CO  .CO  »C0  • 

tO  HOO'OU'OO'OO'O 

M  «0  O  O  O  O 

to  ^  <r  ld  in  ^ 

—i  •  •  in  in  in  in 

<  HO 

Z  o 

<  rH 

>-  ♦  OO'OOOO'OO'  * 

HO  <10  CJO'O-CMOCOO 

MCJ  rH  o  o*  ro 

I>rH  ♦  ♦  •  • 


. .  .C'JOOJrHrOOJO'LO 

UO OJ O Lfl OJrH  rHOJ O O O O O O O O OJ CJ» — ICO t — 1 1 — 1 1 — I  I 

I  I  I  I  I  I  1 

LTicOcocOr-i  00"r0  orOnDr^coorojDCO 

<Tvj-ocncHmc\ihoooooooo . 

. .  •  •  •  •OCOHv^-sPCOOl^ 

mrOCO^H^^HOOOOOOOOCOHHH  I  rH.-HrH 

till  ill 

h  i~i  o  r  o  cj  so  rH  ■  o  <r  ih  o  <r  so  O' o  in  >  t  o  ^  h  cj  rH  ch 

O'O'O  OrHOnDLnO'OJCJLOO'rHOJlOO . 

. CHvOvf'jTOlTlCJvO 

HCJC\JvOHOOOOHOCOK)COrOO  I  }  I  f — I  till 

I 

O'fOO  O'LOOLOLOOLOCOHDOLOrOrOOJO'COCOrHlOvOLOOJ 

O'HD  rH  rH  T  ^  O  HD  rH  LO  O  HD  O' COrO  OHO  O' . 

. CvJrOrHsfHDrOIOLO 

fO  rH  rH  O' rH  rH  LO  rH  O' O' CGI^  O' HDOO 


HOOOO 
CD  sO  CO  rH 


LU 

to  ♦ 

rH  O 

-J  I 

LU 

□ 

O  ‘O  •  • 
HOCOfOrH 


KioonDOO‘Ojr'".cOojcOoojHDCJCJCJcOLOOTOO'0'LO 

LOHDrOC  J  rH  LOf^*  O' O’ fO  LOOJ  HD  C  J  O' rH . 

. O' CO  LOO' CO  OCJ  O' 

rOrHrOh-vOOjcjLOOHOoroojejcoc  j  hh 


LO 

o 

o 


CO  LO  CO  O' CO  1- O  sf  H  r  ^  rH  O  OJ «”  J  o '  rH  rH  vO  r  O  HD  0“  LO  H  vj* 

O'HUlOOl  '•LOO’O'OJLOCJLO  'OOrH . 

. LO  OHDO'CJCJHDOJ 

sj“  xO  CO  H  LO  rH  O' LO  LTir^"-©  O' -O’ COH  CO  rHrHrHrHrHrHrHrH 
rH  HH  r~  I 

rOO'vOCOOJHO'vH'vfrOO'COOKlOO'rHHOOJO'COO'sr 

H-OO'HxOlOHrOO'OLOr-lcOCOrOxO . 

. souvi-vovj’r^ov* 

0  J  O  LO  LO  H  ro  CO  O-' o  o  rH  OJ  <r  O '  H  CO  r-l  rH  rH  rH  H  rH  rH  rH 

OJ  LO  o  o  o  a  j  •  0  O  CO  H  C  U  CO  O  sO  O' r- 1  HD  rH  OJ  sO  CO  OJ  CO  aD 

ro{^o-<ro  O'COrHKir^iO'vocoor^o . 

. voo-o-roo-rncvio* 

LOOJrHO'CJO'HDHDO'CJrOrOrOO  COf^rHrHrHrHrHrHrHrH 

r!Hrl  r— I  rH  i — I  r— I 

<rHOO'Ll_OLLOL.LOLi.LOLLOO  COC  jrOrH'OLOlOO'O'OJrHxOOxf  LO  O  HD  CO  HD  rH  HD  fO  LO 

3  LO  •  U_  LTi  h- Lf  i  U_ '0  OU- O' L'' HD  K)  fO  LO  rOCJf‘')rHrO  CO  rH  CO  OCJ . 

I  hOh  *H“  •  f —  •  h~  *L_  •  •  •  •  »  •  •  •  ♦  •  •  •  •  ♦  •  •  *0  OJi — IrHCOO'rHrH 

LU  »0'NJ  Kl  — J  rH  rH  OJ  CO  OJ  OJ  f  O  rHf^- CO  OJOrO  O' OOJrHrHrHrH  rHrH 

^  r^OJCO  to  -J  rH  rH  rH 

<HDT-  O  O  LO  LOMOrOO'rHOJCOrHrHrOrOOJCJHDrOrHrOOO'O'rHrHOr^LOnD 
-JO  LO  H-K»HKl-JK)-JKlHHDCJsf  rOrHC/'O'O-f^O'COrHO'O’OJCJIO . . 

«J  J  »-J  •— J  «  .KlfOrHrOO'OHDHD 


nr-> 

o 

lu  r- 
x  r- 
to  LO 

DC 


•o  *o  *o  *o  • 

LOi  OJ  LO  c  J  rH  C  J  rH  C  J  r^. 


MMMMCL 

to  in  h  h  D 


O  OJ  OJ  K)  O' OJ  O  r-l  OJ  LO  ro  rH  O' K)  OJ  LO 


CO 


o  • 
oo 
o 


o 

o 


LO  »fO 
•LO 
O 


jo  n  Ki  n 


O'H'nNK'IHKlO  sOKJrHrHCOOO'O 

W^^lCONOl/lOOOOOOOOO . 

. * . r^LocjrHcoinroLT) 

(\JK)OtOHC\JC\JC\JOOOOOOOO  I  I  I  I — (  1111 

I 

r-trHr^KKj-ioo'O  sDr'-O'h-r^rwo 

CM^vjOCOCfjOOOOOOOOO . 

. O^OOsJOOLT) 

OLflCJKIHHOCJOOOOOOOO  I  rHrHrHrH  I  I  I 

till 

CJ  O'NHHvtOJvOO  C\JLT)0'rHv£)K'Kj-0 

O'  coMno^<rcrjooooooooc . 

. •COOCOr-UOCO<Mi,1 

1/lNOLTlCJHHCJOOOOOOOOC'JC'JHCJHHH  I 

I  I  I  I  t  I  I 

mcoconrH^v^-ro  Kivor^coorosoco 


Ki  ‘CO  »c0  ‘CO  ‘CO  *co  *co  • 

Hoa'oa  ocr'ocr  oa'o  if)  lo  <r<rocc)0' mcjnoooooooo 


< 

Q 


♦  o  o  o  o  lo  o  cj  r-  . ocor^sJ-soroor^ 

sd'<rLOLOCJsOfHrs-  K)rOOJ<rH^HOOOOOOOOCJHHH  I 

r-l  I  I  I  I  I  I  I 

•  •  sO  o  o  so  r-»  r^rHorocjvOrHsor-.  LTi^O'rHr^-ojrHCT' 

HO  O'  aCf'C/'OHO^LDO'OOOOOOOC . 

o  . .  •  •O'NOsf  sJTO  LOCO  vO 

rH  HCOCJvOHOOOOOOOOOOO  1  1  I  i — f  !  I  I  I 

I 

OO'OO'OO'OO-OO  •!>  *o>  •  o  ro  O' CMfKMmn  o  cj  o  hn  in  rooo^  co  co  rHirjvO  in  cj 

LOCHCJO'sOO  LOO  0J0sO0-O0'O0'sDrHrHr^OvDrHL0r^rHr0<fh-<\J0va . 

•CjrOrHsfsOKirOLn 


O 

CO  • 
Di  o 
<C0 

LU  rH 

>- 


h-  oooo 

IOOH 

(JjrHCJ 


KlrHrHs^rHrHfOrHrHCJrHOJrOOJrHCO 

rH  rOOOv£)Ox'tOJ(N-OrHrOOOO'r<fCJCOLO<frOO'<rLO 

o  o  LOojrocjHr-nr^O'COi/iO'OsOO'roo  •  . . 

o  O  . .  .O' COLO  O' coo  CO  O' 

•  •  rOrnhOr^vOCOCOLOLOLOfONOLOvOr^sO  rHrH 


LU 


I 


♦  LO 

o  o 
o 


LU 

< 


LU 


r- 

ir> 


•O  *0  ‘O  ‘O  ‘LO  ‘h*.  ♦ 
rHCOrHrHrHCOrHC  JrHf  "*•  f'''-  O 


rH  <f  CO  LOCO  OCQv.no  vOO  O' <fvO  COLT)  CJnCHOD  rOsO'^'LOr^v^* 

o  o  0'HLrtoor-Ln<rc0s0c0LncrmvOH . 

o  O  . . . LO  sO  sO'd'OJ  CO  sO  CO 

>-  ...  vt^rOM/lHCM/lOHHOa'aHO'HHHHHHHH 

Q  r-l  i — I  f — 1 1 — i  i — 1  r — I 

D  -Cj . .  •  •  JOO'.OCOCMrHO'..tinorOr'-r''l,-HCO'OrHNOCMC/-COO'N}- 

I—  OCOJOrH  rH  rH  r-l  rH  i-HKIrH  rH  r^HDI  'rH-C}' JOf^  Kl.OfOLfir'rHOHDlT) . 

to  rHf— »  •  HDLrK3-H5<rrv-vDsT 

CNJ  O  LO  Ifl  K)  CO  O- CH  C\J  rH  <\J  CM  CJ  rH  rH  rH  rH  r— I  rH  rH  rH  rH 
r— 1 i — I  r — I  r — I  r-^  r~i  r~l  r—i  r—\  r—i  r~i  r—^ : — { 

WlClOOO(OH)OlJ  1^-010  cr  O  O' O' HD  rH  CO  HD  CO  CM  CO  CO 

KH^-O-sJ  O'h-COrHHjOv-ct'lv-CO.OCOlO . 

. sO-d’<)'Kl<)-KiCM'J- 

irK'JrHa'CJO'O.OOJOHOHOC-HHrHHHHrHHH 

i — 1 1 — I ' — I  r — 1 1 — I  > — I  p — 1 1 — 1 1 — |  i — | 

I—  rH  O  rv.  U.  Lfl  LL.  Lfl  LL-  lfl  LL.  lfl  U-  LO  LL.  fO  Li_  LOLL-  Lfl  O'  CO  CM  K1  rH  sf  lfl  10  CM  U1  JO  IT)  rH  O'  rH  CO  O  vO  CO  -O  rH  sO  JO  |j~| 

to  o  •u-r'h-cOLi.r'>-r^.u_oLj.oLi.|v-u.r^r--HDr^rofOLOKicjoi'.ocMocM’'-rH . 

HH  CO  <7"  I —  "t—  •(—  -h-  -h-  -IX.  -LL.  -U_  •  * OC’J'-li-l  COO'rHrH 

H-  -ONI  M  M  — I  •  rH  rH  CM  CO  CM  CM  JO  rH  ©O' O' CJ  f^F'' CO  O' rHrHrHrH  rHrH 

O.  HD  CM  to  to  to  -j  Z  Z  rH 

W  Ul  lfl  l/l  Ot-HO<tLO<LnrOCT'rHCMCOrHrHJOOOCOCOLOf^rHr^HnsJ-rHrHOIv-LflHD 

(DO'lfl  (—  01-0— IO— IODWI-  sOXrHXrHCMsJ-JOrHIO'.J-HM^.OHDO.I-CMO'.a-O' . 

— I  * — I  * — •  * — •  *  •(_>  . JOJOrH/OsJ-OMDsO 

M  M  M  M  <C  Q.  OCMCMJOHfCMOrHCMCMrHCMJOsOHfsf 

(O  O  Z  to 


tO  tO  1—  K 


, 


APPENDIX  IV 


SAMPLE  OUTPUT 


LIST IS 


BAPTISTE  LAKE  STUDY:  EIGHT  YEA3S  OF  DATA 


1  hi 


OsD'O 

O'HDO 

OvOO‘1 

LU 

O 

O'  O' 

O'  O' 

O'  O' 

n 

OOO 

ooo 

ooo 

LU 

ooo 

ooo 

coo 

— J 

LA  L£*1 0 

LA  LA© 

LA  LAO 

LU 

Of^-CJ 

OCCH 

okcj 

N* 

OOO 

ooo 

ooo 

n. 

Kl 

Kl 

Kl 

Cl 

O 


LU 

o 

o 

O 

to 

• 

• 
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CLX 

< 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

x  • 

O' 

vf 

A 

A 

A 

A 

O' 

A 

A 

o 

O' 

um 

A  A _ 

A 

A 

CJ 

ro 

<r 

o 

o 

o 

o 

o 

A 

M 

L0 

M 

M 

A 

o 

o 

o 

o 

o 

D-AO 

o 

A 

1 — 1 

o 

o 

o 

o 

o 

o 

o 

A 

r — 1 

O 

CJ 

o 

O 

o 

o 

o 

o 

o 

o 

AH  1 

A 

A 

A 

A 

ro 

A 

CJ 

A 

A 

A 

O 

o 

O' 

A 

M 

A 

A 

vf 

<f 

cm 

ug  J 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

A 

A 

A 

A 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

O 

O 

O 

o 

o 

o 

o 

o 

o 

o 

_ i 

CJ 

ro 

A 

M 

A 

o 

A 

A 

M 

A 

L0 

O' 

MX 

o 

o 

ro 

O' 

A 

o 

o 

o 

o 

o 

A 

A 

CJ 

CJ 

<T 

A 

A 

o 

o 

o 

o 

o 

MU 

DD 

A 

A 

ro 

o 

o 

o 

o 

o 

o 

A 

o 

A 

A 

O 

A 

o 

o 

o 

o 

o 

Mill!  I  I  I  I 


ox 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

c 

o 

o 

o 

O 

LU  AU 
0:M 

oo 

h- 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

AX 

L0 

CJ 

CJ 

1 — 1 

A 

A 

A 

A 

A 

A 

CJ 

A 

A 

A 

A 

o 

A 

A 

A 

A 

CO 

A 

AlL'X 

orr.u 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

ro 

O' 

A 

O' 

O' 

o 

ur 

A 

L0 

A 

A 

A 

IQ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

A 

r — 1 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

I  — 1 


MX 

A 

A 

A 

A 

r — 1 

A 

o 

o 

o 

O' 

OJI 

O 

o 

o 

o 

A 

A 

A 

C0 

A 

A 

o 

o 

O 

o 

A 

o 

O 

O 

o 

o 

o 

o 

A  0-  U 

a. 

o 

o 

o 

o 

O 

O 

o 

O 

o 

o 

A 

A 

o 

A 

o 

o 

o 

O 

o 

o 

o 

o 

x  • 

o 

M 

M 

M 

A 

M 

A 

A 

A 

A 

A 

CJ 

oi-r 

o 

A 

CJ 

to 

O' 

o 

O 

o 

o 

o 

M 

A 

A 

L0 

rH 

A 

A 

O 

o 

o 

o 

o 

AUU 

« 

• 

i — =t 

A 

A 

CJ 

A 

o 

o 

o 

o 

o 

sT 

O' 

o 

A 

O' 

A 

A 

o 

o 

o 

o 

o 

<(  A  A 


o  • 

O' 

cO 

CJ 

•3" 

A 

A 

O' 

A 

A 

O' 

CL  1 —  7Z 

A 

A 

A 

A 

O' 

o 

o 

O 

O 

o 

A 

A 

A 

A 

O 

O' 

o 

o 

O 

o 

o 

<oo 

>CL 

O' 

CJ 

A 

M 

A 

o 

o 

o 

O 

o 

O' 

A 

A 

O' 

A 

CJ 

o 

o 

o 

o 

o 

UJ 

A 

r — 1 

A 

A 

A 

O' 

O 

CJ 

A 

A 

'd- 

A 

O' 

A 

ro 

A 

O 

vf 

A 

M 

o 

M 

A 

o 

o 

o 

u 

A 

M 

A 

M 

A 

A 

CJ 

A 

A 

A 

M 

CJ 

A 

A 

O' 

A 

O' 

o 

o 

o 

ujD 

XU 

A 

A 

O' 

CJ 

A 

o 

-3- 

A 

o 

CJ 

A 

A 

A 

O' 

A 

L0 

A 

o 

A 

CJ 

CJ 

CJ 

X 

A 

A 

# 

o 

o 

O 

O 

o 

o 

o 

O 

o 

o 

O 

O 

O 

O 

o 

O 

O 

o 

O 

o 

o 

o 

CL 

A 

O' 

M 

A 

A 

A 

<T 

O' 

O' 

A 

A 

A 

<r 

<T 

A 

CJ 

O' 

A 

o 

o 

o 

XU 

UJ 

CJ 

A 

A 

CJ 

A 

CJ 

O 

o 

A 

A 

A 

'd' 

CJ 

O' 

A 

A 

M 

A 

A 

A 

h- 

A 

i — 1 

A 

A 

1 

A 

rH 

• 

1 

A 

A 

t — i 

A 

1 

A 

1 

1 

1 

LUO 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

ro 

Hv)  LOAMAO'OAcjA<fLOAMAO'OACJA<rAA 

I NkNf\fsc0(0c0c0C0(0C0C0C0c0OC'(7'0'0'0'^ 
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APPENDIX  V 


BAPTISTE  LAKE  HYDE  CLGG1C  D  AY  A 


Baptiste  Lake 


Surface  Elevations 


(all  values  are  recorded  on  the  last  day  of  the  month) 


(units  are 

meters  . 

aheve  mean  sea  level) 

19  72 

1973 

1974 

1975 

1976 

1977 

1978 

Jan 

— 

— 

— 

— 

— 

— 

— 

Fee 

- 

- 

- 

- 

- 

- 

- 

Mar 

- 

— 

- 

- 

- 

— 

- 

Apr 

•? 

■? 

578. 867 

578.  02  6 

577.962 

5  78.  059 

Ma  y 

578. 066 

577. 809 

578.  18  i 

577. 822 

577.  83  1 

> 

578. 062 

J  un 

577. 989 

578.  169 

5  77.  9J4 

57  7.943 

577.843 

578. 0 1 7 

5  78.  06  7 

Ju  1 

577. 683 

578.  035 

578. 367 

578.  29  7 

577.  8  1  9 

578. 257 

5  78.  107 

A  u  g 

5  7  7.  7  62 

57  7.  968 

5  78.  148 

5  78.  34 3 

577. 858 

5  78. C5  3 

578. 139 

Sep 

57  7.  72  7 

5  77.  93  7 

5  78. C23 

5  78.  224 

5  7  7.  86  1 

5  78.  09  9 

5  78.  765 

Oc  t 

5  77.  7  C  9 

577. 9  16 

577.  9  19 

5  7  7.  9  77 

577. 843 

5  76.  08  1 

57o,  27  1 

Nov 

- 

- 

- 

- 

— 

— 

- 

Dec 

- 

- 

- 

- 

- 

- 

- 

-  :lake  was  frozen 
?  :lake  level  not  recorded 


* 


Recorded  icy  the  Water  Survey  of  Canada 
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Ice  Thicknesses  on  the  Iso  Basins  of  Baptiste  Lake 


Baptiste 

S  o  u  t  h 

Eaptiste 

North 

Nov 

2  0, 

1  9  7  c 

frozen  over 

Nov 

10, 

1976 

frozen  < 

o  ver 

Dec 

15, 

1976 

20-25  cm 

Dec 

1, 

1976 

23 

cm 

Jan 

5, 

1977 

36-46  cm 

Dec 

15, 

1976 

20-26 

cm 

Feb 

16, 

1977 

49  cm 

Jan 

c; 

-  9 

1977 

3  1-33 

CO 

Mar 

5, 

1977 

5  1  cm 

Jan 

15, 

1977 

3  1 

CO 

Mar 

23, 

197  7 

6  1  cm 

Feb 

2, 

1977 

38-46 

cm 

Apr 

25, 

1977 

ice  out 

Fe  fc 

16, 

1977 

5  1-54 

cm 

Mar 

9, 

1977 

51 

cm 

Mar 

23, 

1977 

56-6  1 

cm 

Apr 

O  c 

^  i 

197  7 

ice  - 

out 

*ice 

thickness  supplied  by  D. 

Beliveau, 

i«ater 

Quail ty 

Control  Etanch,  Alterta  Environment. 


- 

m 


wmri  iisinv*  on<mv#q 


a yu  wnnootw 


«>n  nn 


•  «l  Hit  JO  ONIJMTWO  S13V143M 


575 
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' I  S' lU'D  UJ  N0I1VA3H3  3XV1 


FIGURE  30.  BAPTISTE  LAKE  AREA  -  CAPACITY  CURVES 
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APPENDIX  VI 


BAPTISTE  LAKE  H YDEOCH EMICAL  LATA 
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South  Basin  of  Baptiste  Lake,  sampled  from  top  of  lake,  1976-1977- 
analyses  supplied  by  Water  Quality  Control  Branch,  Alberta  Environment 


154 


155 


South  Basin  of  Baptiste  Lake,  sampled  from  middle  of  lake,  1976-1977- 
analyses  supplied  by  the  Water  Quality  Control  Branch,  Alberta  Environment 


!i 


156 


i 


157 


South  Basin  of  Baptiste  Lake,  sampled  from  bottom  of  lake,  1976-1977. 
analyses  supplied  by  the  Water  Quality  Control  Branch,  Alberta  Environment 


I 


>1 
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North  Basin  of  Baptiste  Lake,  sampled  from  top  of  lake,  1976-1977- 
analyses  supplied  by  the  Water  Quality  Control  Branch,  Alberta  Environment 


160 


1 61 


j. 


North  Basin  of  Baptiste  Lake,  sampled  from  middle  of  lake,  1976-1977- 
analyses  supplied  by  the  Water  Quality  Control  Branch,  Alberta  Environment 


162 


Ji 
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North  Basin  of  Baptiste  Lake,  sampled  from  bottom  of  lake,  1976-1977. 
analyses  supplied  by  the  Water  Quality  Control  Branch,  Alberta  Environment 
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APPENDIX  VII 

GROUNDWATER  CHEMISTRY  DAI  A 
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Groundwater  Chemistry:  Peizometer  IB 
LSD  4-  15-6  7-24- K4M 


(all  analyses  in  mg/1) 


Date 

Ca 

ag 

Na 

K 

Cl 

S04 

HC03 

CO  3 

i 

f  e 

S  l 

1977 

Aug 

10: 

3  1 

1 1 

452 

2.3 

260 

1  1  7 

754 

<5 

0.64 

1.45 

10.  2 

Aug 

16: 

3  1 

1 1 

452 

2.6 

269 

1  18 

747 

6 

0.  66 

10.47 

8.5 

Aug 

23: 

7 

10 

4  10 

3.5 

2  C  7 

150 

638 

<5 

0.  28 

8,70 

6.  1 

Au  g 

30: 

7 

1 1 

354 

-5.3 

134 

1  79 

480 

40 

0 . 4  ^ 

0.  16 

6  .  2 

Sep 

6: 

4  7 

1  4 

375 

3.6 

128 

2  3  1 

7  1c 

6 

0.  55 

6,36 

8.  7 

Se  p 

14: 

4  5 

14 

3  75 

3.6 

123 

222 

6  St 

7 

0. 57 

6.  34 

8.  7 

Se  p 

20: 

4  6 

15 

365 

3.  1 

1  12 

236 

6  85 

- 

0.52 

<.05 

8.8 

Sep 

27: 

48 

1  6 

375 

3 . 3 

10  1 

224 

733 

— 

0.56 

1  0.  4  7 

9. 6 

Oct 

4: 

16 

S 

2  18 

1 . 8 

62 

1  40 

407 

— 

0.  52 

0.  3  1 

4.  6 

Oct 

19: 

9 

14 

365 

3.4 

100 

22o 

605 

<5 

0.50 

0.  18 

7.  5 

Oct 

26 : 

25 

14 

350 

3.  1 

93 

2  14 

635 

— 

0.  40 

6. 30 

8.  5 

Nov 

1: 

50 

16 

365 

3.3 

92 

2  33 

735 

— 

0.52 

7.95 

9.  6 

Nov 

8: 

3  5 

15 

3  6  5 

3.4 

99 

2  30 

7  16 

<5 

0.53 

7.  52 

9.  6 

Nov 

15: 

5  1 

16 

365 

3.0 

98 

24  1 

732 

<5 

0.  4  7 

7.53 

1  0.  5 

Nov 

22: 

48 

16 

360 

4.  ! 

103 

2  35 

7  19 

6 

0.  53 

6.  97 

9.  8 

Nov 

30: 

5  0 

15 

370 

3.2 

95 

229 

697 

1  5 

0.65 

7.  67 

9.  6 

Dec 

7: 

4  8 

15 

356 

2.9 

95 

23  1 

732 

— 

0.52 

8.  04 

9. 4 

De  c 

2  1: 

27 

15 

390 

3.5 

94 

224 

6  69 

— 

0.39 

9.2  1 

5.  8 

1978 

Jan 

3: 

47 

15 

365 

3.2 

93 

233 

723 

<5 

0.53 

6.  86 

8.  6 

Jan 

10: 

48 

15 

356 

2.9 

95 

234 

7  13 

9 

0.56 

7.  70 

9.6 

Jan 

2  5: 

34 

15 

356 

2.8 

9  1 

227 

737 

— 

0.  52 

9.05 

9.  2 

Jun 

27: 

4  4 

15 

359 

3.  1 

94 

230 

744 

- 

0.  50 

9.9  0 

9.  6 

Ju  1 

25: 

4  9 

14 

350 

2.6 

92 

244 

733 

<5 

0.46 

11.24 

9.  1 

Aug 

14: 

44 

16 

354 

2.8 

90 

246 

736 

<5 

0.  50 

11.32 

9,  5 

♦analysed  cy  the  Pollution  Control  Lab,  Alberta  Environment 


in  in  o  vo  in  co  vo 


- 

' 

Groundwater  Chemistry:  Piezometer  IE 
(all  analyses  in  mtj/L,  except  conductivity 


which 

is  in 

microsi erne ns/cm) 

Date 

19  7  7 

IDS 

Cond 

TIC 

IOC 

NC3+N02 

Ik  N 

C-PC4 

T-P04 

Aug 

10 

1269 

2  150 

150 

7 

0.  0  10 

2.  36 

0.  167 

0.7  1 

Aug 

16 

1 2  5  £ 

2200 

149 

13 

0.  0  10 

2.  63 

0.  054 

0.  <43 

Aug 

23 

1  107 

1925 

129 

30 

0.  0  1b 

2.69 

0.  047 

1.  16 

Aug 

30 

965 

1604 

126 

24 

0.  003 

1.8  1 

0.  03  1 

0.  4  1 

Sep 

6 

1  158 

1850 

129 

12 

0.033 

2.  73 

0.  040 

0.5  8 

Se  p 

1  4 

1  133 

1825 

129 

<2 

0.  0  14 

2.95 

0.  047 

C .  o  3 

Sep 

20 

1115 

1  850 

134 

12 

0.  0  16 

3.  13 

0 .  0  4  6 

1.26 

Sep 

27 

1  128 

1850 

128 

1  7 

0.  005 

2.60 

u  •  0  2  7 

0.68 

Oct 

4 

646 

1  106 

87 

4 

^  Ago 
\Jm  U/  U 

2.6  9 

0.015 

0.  75 

Oct 

19 

1  009 

1666 

137 

3 

0.  0  16 

1.16 

0,067 

0.3  2 

Oct 

26 

10  1  1 

16  56 

130 

18 

0.  034 

1.80 

0.052 

0.  50 

Nov 

1 

112  1 

1850 

144 

7 

0.  008 

2.  70 

0,  2  55 

0.  69 

Nov 

0 

1  103 

1850 

1  3  6 

8 

0.0  11 

64.  80 

0.  016 

0.  72 

Nov 

15 

1  137 

1870 

129 

26 

0.  0  17 

2.  5fc 

0.  0  16 

0.  72 

Nov 

22 

1  126 

1  7  33 

1  4  7 

15 

0.  025 

2.69 

0.036 

0.  6 3 

Nov 

30 

1  120 

1  743 

150 

9 

0.022 

13.33 

0.016 

0.  74 

Dec 

7: 

1  108 

1773 

129 

36 

0.  0  26 

2.  77 

0.0  16 

0.  88 

Dec 

2  1 

108  2 

169  1 

138 

12 

0.  02b 

2.  42 

0.  067 

0.59 

Jan 

3 

1  104 

1753 

154 

<2 

0.  0  1  1 

2.50 

0.  0  19 

0.56 

Jan 

1  0 

1112 

1  750 

140 

16 

0.  02  1 

2.  66 

0.  027 

0 .  o  3 

Jan 

25 

1089 

1  754 

129 

37 

0.  0  13 

2.  44 

0.  02  7 

0.  70 

Ju  n 

27 

:  1111 

1  763 

137 

- 

0.  0  17 

2.29 

0.  04  1 

- 

Jul 

25 

:  1115 

1  79  1 

150 

- 

0.  068 

2.  1  0 

0.  057 

- 

Aug 

1  4 

:  1118 

1  760 

156 

9 

0.  002 

2.68 

0.  044 

0.  72 

♦analysed  by  the  Pollution  Control  Lab,  Alberta  Environment 


168 


Groundwater  Chemistry:  Piezometer  2B 
LSD  1-2  1-66-24-K4H 


(ail  analyses  in  mg/L) 


Date 

1977 

C  a 

Mg 

Na 

K 

Cl 

S04 

HC03 

CO  3 

fe 

S  l 

Aug 

10: 

<2 

2 

298 

1.3 

22 

8  1 

707 

1  4 

0.40 

8.  19 

7.  3 

Aug 

16: 

Aug 

23: 

Aug 

30: 

<2 

1 

3  12 

1.  7 

25 

7  5 

667 

27 

0.42 

0.  33 

3.  2 

Sep 

6 : 

2 

2 

3  08 

1.7 

23 

75 

678 

2  5 

0.  42 

3.38 

2.  5 

Sep 

14: 

4 

2 

3  18 

1 . 8 

24 

73 

6  6  7 

3  1 

0.  43 

3.43 

2.  4 

Sep 

20: 

2 

2 

3  08 

1.  7 

22 

79 

654 

1  4 

0.37 

<.05 

2.  1 

Sep 

27: 

3 

O 

3  1  7 

1.6 

22 

74 

685 

2  I 

0.47 

4.37 

2.3 

Oct 

4: 

<2 

1 

1  94 

0.9 

13 

40 

4  38 

13 

0.53 

0.40 

1  .  4 

Oct 

18: 

Oct 

26: 

<2 

2 

300 

1  .  7 

22 

64 

632 

3  S 

0.40 

1 .02 

2.  2 

Nov 

1: 

3 

2 

288 

1  .  6 

22 

72 

644 

43 

0.44 

1.  59 

1.  4 

Nov 

8: 

3 

2 

302 

1.7 

23 

83 

655 

37 

0.43 

1.59 

1.  0 

Nov 

1  5: 

2 

1 

3  17 

1 .  7 

30 

84 

64  1 

42 

0.4  1 

1.63 

1.3 

Nov 

2  2: 

<2 

2 

300 

1.6 

24 

64 

640 

42 

0.  43 

1 .56 

1 . 6 

Nov 

30: 

<2 

2 

3  15 

1.6 

24 

6  1 

638 

45 

0. 58 

1.  53 

1.  2 

Dec 

7: 

<2 

2 

3  08 

1.6 

23 

66 

647 

40 

0. 4  1 

1.61 

1.  3 

Dec 

2  1: 

<2 

270 

1.4 

2  1 

53 

579 

35 

0.38 

1  .47 

0.  8 

1  9  78 

Jan 

3: 

<2 

2 

302 

1.  7 

23 

60 

63  1 

49 

0.4  1 

1.42 

1.  0 

Jan 

10: 

<2 

2 

308 

1 . 6 

25 

57 

6  13 

58 

0.  46 

1 .50 

1.0 

Jan 

25: 

<2 

1 

298 

1  .  6 

23 

58 

634 

47 

0.4  1 

1.69 

1. 4 

Feb 

13: 

<2 

2 

300 

1.4 

25 

62 

624 

53 

0.  4  1 

1.33 

1.  0 

Ma  r 

8: 

<2 

2 

298 

1 . 6 

24 

57 

5  13 

57 

0.42 

1  .  14 

1.  0 

Jun 

27: 

2 

<  1 

3  15 

1  .  8 

27 

58 

698 

18 

0 .  4  4. 

2.  1  1 

4.  1 

Ju  1 

2  5: 

<2 

2 

3  20 

1.6 

28 

65 

687 

24 

0.42 

1.78 

3.  6 

♦  analysed  by  the  Pollution  Control  Lab,  Alberta  ^Environment 


' 

Groundwater  Chemistry:  Piezometer  2 B 


(ail  analyses  in  wg/L,  except  conductivity 
which  is  in  microsiemens/cm) 


Date 

IDS 

Cord 

TIC 

ICC 

NC3+N02 

1‘KN 

0-P04 

T-  ?0  4 

]  9  7  7 

Aug 

10 

768 

1295 

148 

4 

0.0  06 

1.46 

0.  008 

0  .  0  6 

Aug 

16 

Aug 

23 

- 

- 

132 

<2 

0.  0  17 

0.  67 

0.  03  1 

<.05 

Aug 

30 

772 

1282 

133 

8 

0.  003 

1. 26 

0.00  7 

<.05 

Se  p 

6 

770 

1280 

— 

- 

0,  006 

1.6  7 

0.029 

<.05 

Sep 

1  4 

78  1 

1  100 

125 

1  1 

0.0  12 

1.65 

0.  04  1 

<.05 

Sep 

20 

752 

1280 

134 

12 

0.  0  10 

1.  66 

0.0  15 

<.  05 

Sep 

27 

778 

1280 

1  3  1 

10 

0  •  004 

1.80 

0.  0  10 

<.05 

Oct 

4 

479 

83  1 

122 

<2 

0.  0  15 

1.2  1 

0.  0  13 

<.05 

Oct 

19 

Oc  t 

26 

739 

1254 

123 

2  7 

0.0  17 

1.  18 

0.  0  14 

<.05 

Nov 

1 

749 

1310 

122 

9 

0.  0  10 

1.  76 

0.0  11 

<.  05 

No  v 

8 

773 

1350 

138 

4 

0.  0  12 

1.  72 

0.  007 

<.05 

Nov 

15 

794 

13  10 

133 

1  4 

0. 009 

1.6  1 

0.011 

<.  j5 

Nov 

22 

750 

1288 

142 

10 

0.  022 

1.8  1 

0.  022 

<.  05 

Nov 

3  0 

764 

1268 

136 

<2 

0.  0  19 

2.  85 

0.  0  18 

<.  05 

De  c 

7 

760 

12  76 

1  2  6 

22 

0.017 

1.  76 

0.  0  29 

<.  05 

Dec 

2  1 

:  668 

1  132 

134 

12 

0.  024 

1.  73 

0.  022 

<.05 

19  7  6 

Jan 

3 

75  0 

1263 

1  42 

14 

0.  0  17 

1.  76 

0.  0  18 

<.05 

Jan 

1  0 

753 

126  1 

144 

6 

0.  03  1 

1.  75 

0.  0  16 

<.05 

Jan 

25 

742 

126  6 

12  1 

33 

0.  0  12 

1.6  1 

0.  035 

<.  05 

Feb 

13 

:  75  1 

1272 

1  4  C 

10 

0.  0  13 

1.  66 

0.  038 

<.  05 

!1ar 

8 

:  743 

1275 

134 

<2 

o 

• 

o 

1.  76 

0.  0  47 

0.0  7 

Jun 

27 

:  767 

1280 

125 

<2 

0.011 

1 .68 

0.  0  19 

0.  0  3 

Ju  1 

25 

:  780 

1285 

I  5  C 

— 

0.  057 

1.  62 

0.  0  14 

— 

^analysed  fcy  the  Pollution  Control  Lab,  Alberta  Environment 


Groundwater  Chemistry:  Piezometer  3B 
LSD  6-2  I-66-24-W4JS 

(all  analyses  in  mg /l,  except  conductivity 
which  is  in  fflicrosiemens/ciii) 


Date  Ca 

19  77 

Na 

K 

Cl 

304 

tiCG  3 

CC  3 

1 

S  e 

si 

Oct  8 :  8.4 

2 

28  I 

1 . 6 

26 

28 

694 

<5 

0.55 

0.88 

11.5 

Da  te 

19  77 

IDS 

Con  d 

TIC 

T  CC 

N03+NC3 

Ik  N  0-P04  1-P04 

Oct  d : 

6  8  9 

12  10 

- 

- 

<0.  1 

-  -  - 

^analysed  by  the  Pollution  Control  Lab,  Alberta  Envirnment 


. 

Groundwater  Chemistry  from  the  Water  laLle  Veils 

(ail  analyses  in  mg/L) 


Well 

No. 

Location 

(k4M) 

Sampling 

Date 

Ca 

Mg 

Na 

K 

HCG3 

SC4 

Cl 

1 

1-26-66-24 

9/8/78 

33  7 

1  15 

30 

5.9 

453 

263 

500 

2 

3-26-66-24 

8/8/78 

1  1  3 

4  1 

23 

4.  6 

564 

55 

<  1 

3 

4 

16-27-66-24 

6-22-66-24 

damaged 

8/8/78 

120 

42 

14 

3.  7 

584 

6  1 

2 

C 

5-22-66-24 

8/8/78 

38 

1  1 

189 

3.  5 

622 

6  1 

2 

6 

7 

7-2  1-66-24 
10-2  1-66-24 

lost 

8/6/78 

1  76 

4  7 

29 

1 . 6 

500 

3  16 

C 

a 

12-21-66-24 

8/8/ 7  8 

228 

1  7 

1  7 

2.  2 

292 

3  1 

<  1 

q 

2-28-66-24 

8/8/78 

12  1 

4  I 

27 

2.  2 

56  6 

137 

<  1 

10 

4-28-66-24 

8/8/7  6 

1  0 

2 

7 

7.  7 

42 

<  10 

<  1 

1  1 

4-33-66-24 

9/8/76 

10  7 

72 

12 

6.8 

790 

c  c; 

3 

12 

9-29-66-24 

8/8/78 

96 

25 

28 

3.  2 

422 

61 

2 

13 

8-33-66-24 

8/8/76 

98 

17 

8 

2.  0 

36  0 

33 

<  1 

14 

1-  3-67-24 

8/8/78 

52 

1  7 

168 

5.3 

582 

133 

<  1 

1  5 

16-  3-67-24 

8/8/7  6 

158 

52 

67 

5.8 

77  1 

2  4  4 

2 

16 

1 5- 1 0-6  7-24 

8/8/78 

103 

26 

5 

2.  2 

385 

55 

2 

1  7 

8-  15-67-24 

9/8/78 

56 

18 

130 

3.  7 

52  6 

65 

2 

18 

19 

10-  1  1-6  7-24 
u-13-67-24 

destr oyed 
9/8/78  138 

4  9 

36 

4.  1 

590 

55 

4  5 

20 

4-  1  2-67-24 

8/8/78 

136 

46 

34 

6.  7 

636 

1  13 

3 

2  1 

1 6-34-66-24 

9/8/78 

123 

47 

7 

0.6 

645 

28 

2 

22 

2-  2-66-24 

8/8/76 

84 

25 

3 

4.6 

632 

6  1 

<  1 

23 

1  1-35-66-24 

9/8/7  8 

94 

34 

18  1 

6.6 

574 

158 

<  1 

♦analysed  cy  the  Pollution  Control  Lac,  Alberta  Environment 


■ 

■ 

Groundwater  Cheiristry  fro®  Water  Table  Wells 

(all  analyses  in  rcg/L,  except  conductivity 
which  is  in  nucrosiemens/cn) 


Well  Location 


No. 

(  W  4  M ) 

Si 

Fe 

1 

N02+N03 

TDS 

Con  d 

1 

1-26-66-24 

19.  0 

0.  55 

0.  14 

0.10 

1494 

26  4  0 

2 

3-26-66-24 

16.  4 

C.  24 

0.  15 

<.05 

5  16 

9  1  9 

3 

u 

16-27-66-24 

6-22-66-24 

dama  ced 

15.  2 

1.  15 

0.  13 

0. 06 

532 

942 

C 

5-22-66-24 

11.6 

0.  26 

0,  36 

0.06 

6  1  1 

10  18 

6 

7 

7-2  1-66-24 
10-2  1-66-24 

destroyed 

11.7  1.40 

0.  1  7 

<.  0  5 

82  1 

120  6 

8 

12-21-66-24 

1  7.  0 

0.  75 

0.  13 

2.  1  6 

2  16 

5  14 

9 

2-28-66-24 

12.  9 

1. 29 

0.  2  1 

0.  06 

609 

102  1 

10 

4-28-66-24 

9.  5 

4.  00 

0,05 

0.84 

57 

98 

1  1 

4—  33- 6  6  —  24 

18.  3 

0.  74 

G.  12 

0.  15 

645 

1  167 

12 

9-29-66-24 

12.4 

1.  2  1 

0.  13 

0,  1  2 

422 

738 

]  3 

8-33-66-24 

18.2 

1.  10 

0.  19 

0.06 

334 

603 

14 

1-  3-67-24 

13.  4 

C.  56 

0.  29 

0.  10 

663 

1088 

15 

16-  3-67-24 

14.5 

0.  43 

0.  13 

0.  08 

909 

144  1 

16 

15-10-6  7-24 

12,  1 

0. 49 

0.  2  1 

2.  64 

395 

68  7 

17 

8-  15-67-24 

14.  4 

G,  49 

0.  20 

0.  10 

534 

875 

18 

19 

10-  1  1-67-24 

4-  1  3-67-24 

destroyed 

17.5  0.37 

0.  1  1 

0.  08 

6  It 

1073 

20 

4-  12-6  7-24 

14.5 

0.  17 

0.  08 

0.  1^ 

652 

1080 

2  1 

1  6- 34- 66-24 

22.  1 

0.  12 

0.  2  1 

0.  12 

529 

843 

22 

2-  2-66-24 

5.  1 

0.  17 

0.  12 

0.  06 

550 

1002 

23 

1  1  -  35-t>  6-24 

1  7.  0 

0.  93 

0.  18 

0.  12 

79  1 

1226 

♦analysed  Ly  the  Solution  Control  Lai, 


Alberta  Environ cent 
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Groundwater  Chemistry:  Water  iable  Weils 

(ail  analyses  in  mq/L,  except  conductivity 
which  is  in  microhmos) 

Well  Location 


No. 

{ W  4  M ) 

Na 

HCO  3 

604 

Cl 

C  o  it  d . 

1 

1-26-66-24 

29 

725 

260 

a  io 

35  10 

2 

3-26-66-24 

22 

54  1 

54 

<2.  5 

84  6 

3 

1 6—  2  7-6  6-2  4 

9 

4  1  7 

7 

<2.5 

620 

4 

6-22-66-24 

16 

552 

62 

<2.  5 

848 

c 

5-22-66-24 

2  14 

597 

6  5 

<2. 5 

960 

6 

7-2  1-66-24 

43 

767 

22 

<2.  5 

1090 

7 

10-2  1-66-24 

33 

493 

320 

5.  0 

1095 

8 

12-2  1-66-24 

7 

255 

12 

<2. 5 

432 

9 

2-28-66-24 

32 

- 

167 

<2.5 

995 

10 

4-26-66-24 

6 

— 

20 

3.  0 

98 

1  i 

4-33-66-24 

9 

_ 

52 

<2,  5 

1050 

12 

13 

9-29-66-24 

8-33-66-24 

dest ro  yed 

7 

3  1 

<2.  5 

575 

14 

1-  3-67-24 

2  1 

- 

59 

<2.  5 

890 

15 

16-  3-67-24 

73 

— 

250 

<2. 5 

136  5 

16 

15-  1  0-67-24 

4 

— 

50 

<2.5 

64  5 

1  7 

8- 1 5-6  7-24 

95 

- 

34 

<2 .  5 

670 

18 

19 

20 

10-  1  1-67-24 
4-  13-67-24 
4-  12-67-24 

lost 

destroyed 

174 

1  50 

<2  .  5 

1000 

2  1 

16-34-66-24 

6 

— 

20 

<2.5 

805 

22 

2-  2-66-24 

2 

6  1 

<2.5 

1020 

♦analysed  at  the 

Geoloqy 

Dept  , 

Oniv.  or 

Alberta 

on  14/10/7 

. 
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APPENDIX  VIII 
PHYSICAL  GECUNDKATEE  DATA 
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0  fc  s  e  r 

Well 

1 

2 
— N 

J 

4 

5 

6 

7 

8 
9 

JO 

1  1 
12 
1  3 

14 

1  5 

16 
1  7 
18 

1  9 

20 

2  1 

2  2 
2  5 

24 

25 


Summary  of  the  Shallow  lest-Hole  Dr il liny  Program 


• 

lo  tal 

Depth  to 

Slotted 

Aquifer 

#  Location 

Depth 

bater 

Interva 1 

Lithology 

1-26-66-24- W  4  M 

7,  6m 

6.  4m 

4.6-  7.6  m 

sa  ndy-tili 

3-26-66-24- W4M 

6.  1  m 

5.  2m 

1.5-  6.1m 

san  d 

1  6-2  7-66-24- W  4  M 

4 .  6m 

2 . 6m 

1.5-  4.6m 

sand 

6-2  2-6  6-2  4- W  4M 

4.  6m 

1.  5m 

1.6-  4.  om 

cla  yey-t ill 

5  -  2  2-  6  6  -  2  4-  W  4  M 

7.  6m 

6.4m 

5.2-  7.6m 

san  dy-t ill 

7-2  1-66-24-H4M 

3.  Om 

0,  9m 

0.6-  3.0m 

sandy- till 

10-2  1— 66— 24-W4H 

3 .  0  m 

0.  6m 

0.3-  3.0m 

sand 

12-2  1-66  —  24-W451 

6,1m 

4.6  m 

1.  5-  6.  lm 

sand 

2- 2  8—  6  6-24- W  4M 

4,  6m 

2.  7m 

0.9-  3.7m 

sand 

4-28-66-24- »4M 

9.  Im 

6.  7m 

4.6-  9.  lm 

sand 

4-33-66— 24-w4M 

12.  2m 

1  1.  6m 

7.  6-12. 2m 

cla  yey-t ill 

9-29-66 -2 4-W4M 

3,  Om 

1. 8c. 

0.  6-  3, 0m 

sand 

8-33-6  6-24-  W  4  M 

12.2m 

11.0m 

6.  1- 12. 2m 

clayey-till 

1-  3-67-24-W4M 

10.  7m 

9,  5m 

6 .  1  -  1  0 .  7  m 

clayey-till 

16-  3-6 7-24- W4M 

9.  1m 

7.  9m 

3.0-  9,1m 

clayey-till 

15-  10-67-24- W4M 

4,6  m 

3.0m 

1.5-  4.  6m 

sand 

8-  15-67-24-W4M 

12.2m 

1  1 .  Om 

6 .  1  —  12.  2  ie 

clayey-till 

10-  1  1-67-24- W4M 

10.  7 hi 

10.  !m 

6 .  1  —  10.7m 

clayey-till 

4-  1  3-67-24- W4K 

6.1m 

4 . 6m 

2.4-  6.1m 

clayey- till 

4-  I  2-67-24-W4M 

1 5 . 2m 

1  4,  3m 

12.  5-  15. 2m 

clay ey-t ill 

16-34-66-24- H4M 

9.  1m 

7.3m 

5.5-  9.1m 

clayey-till 

2-  2-66-24- SUM 

3.  0m 

0,  9m 

0.6-  3.0m 

sandy-till 

1  1  -35-66-24- H4M 

12.  2m 

10.  7m 

7.  6-  1  2 . 2  m 

sandy -till 

16-  1  7-66— 24-W4M 

15.2m 

- 

no  piezometer 

14-28-66-24- W  4  M 

9.  1m 

8.2  m 

no  piezometer 

. 

■ 

Summary  of  the  Deep  Test-Hole  Drilling  Program 


Test 

Hole 

No. 

Dril ling 

Location  Date 

i 

Depth 

bedr ock 
Depth 
(Elev.  2) 

Eemar ks : 

76-  1 

4-  1  5-6  7-24- fo'  4M 

1 

54.  9m 

43 . 0m 
(54  5  m) 

flowing;  14 

L/min 

76  —  2 

1 -  20-6  6-2  4- fc  UM 

1 

102.  Im 

98.  2m 
(485m) 

flowing;  45 

L/min 

76-3 

13-36-66-24- E  4  M 

S 

56.  7m 

4  1.2ic 
( 545m) 

7  6-4 

1 2-26-66-24-W4M 

1 

36.6m 

20,  Im 
( 560m) 

76-5 

12-36-66-24- H4M 

1 

7 3.  2m 

62 . 5m 
(525m) 

76-8 

I3-33-66-24-W4M 

1 

30.  5m 

20 .  4m 
(580  m) 

76-  1  1 

14-  12—6  7—  2 4  —  W 4  M 

1 

42.  7m 

7.6  m 
(590m) 

77-  1 

6-2  1-66-24  - W  4?! 

2 

2  1.9m 

17.  Im 
(5  70m) 

77-2 

6-2  1-6  6-24- H4M 

£ 

5  1 . 8m 

47.2m 

(540m) 

flowing;  40  L/min 
static ;  2.1m  ags3 
piezometer  3B 

77-3 

4-  1  5-6  7-24- W4M 

2 

29.  0m 

(-) 

flowing;  15  L/min 
static;  7.0m  ags3 
piezometer  IB 

77-4 

7-  1  0-6  7-24- H4M 

2 

19.8m 

13.  Im 

(565m) 

77-5 

8-  1  0-6  7-24- '*41 

2 

2  7.4m 

1 4.  6m 
(570m) 

7  7-6 

1-20-66— 24- W4M 

2 

82.  3m 

59.4m 

(530m) 

flowing;  10  l/min 
static ;  9.8m  ags3 
piezometer  2B 

i;  1  -  Dec  6  to  Dec 

2  =  Jul  17  to  Aug 
2:  elevations  are  m 

10,  1976 

3,  1977 

m •  s  •  X • 

3:  ags  =  above  ground  surface 
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Depth  to  Water  in  the  Water  lable  Wells 


Well 

No, 

Date: 
Jun  8 

1976 

Jun  28 

Auq  8 

Auq  1 2 

OCt  1  4 

1 

2.79m 

2.  90m 

— 

3.30m 

2.46m 

2 

0 .  3  3  m 

C.  86m 

1 , 3  8  m 

- 

0.  38m 

3 

2.  6  9m 

2.7  3m 

2.81m 

- 

1. 94m 

4 

0.88m 

1.  1  8m 

1.29m 

- 

0 .  7  5  m 

5 

0.6  5m 

1.  09m 

1.40m 

— 

0.61m 

6 

0 .  79m 

C.  95m 

— 

— 

0.  35m 

7 

0.  53m 

0.  8  la 

0.  95m 

0.53m 

8 

5.  0  7m 

5 .  16m 

5.06m 

- 

4.97m 

9 

1 .  5  2  in 

1 .  6  8  m 

1.73m 

- 

l .  35m 

10 

5.  8  7m 

5.  63m 

5  .  b  7  m 

— 

5.  7 2m 

1  1 

1  1 .  04® 

1  1.  13m 

— 

1  1 . 09m 

1  0.  9  9m 

12 

2.24m 

2.  3  0m 

2.08m 

- 

- 

13 

3,  76m 

4.  09m 

4.05m 

- 

2.74  m 

14 

4.10m 

4.11m 

4.00m 

- 

2.51m 

15 

5.  4  5m 

5.  56m 

5.  79m 

— 

4.  6 5m 

16 

1  .  70m 

1.  9  1m 

2.  06  m 

— 

1.64m 

17 

2.90m 

2.  72m 

- 

2.  8  5m 

1.80m 

13 

8.  0 5m 

- 

- 

- 

— 

1  8 

1.2  1m 

1.  35m 

1 .  52  m 

- 

0.  69m 

20 

8 . 4  5m 

8.  4  1m 

8.44m 

— 

8.57m 

2  1 

2.  7  7m 

2.  86m 

— 

2.6  0m 

2.01m 

22 

0.7  8m 

1.  09m 

1.26m 

- 

1.0  1m 

23 

1  1. 0  7m 

1  1.  13m 

- 

- 

- 

. 

• 

1  78 


Water  Levels  in  the  Piezometers 


Date: 

1  A 

Depth1 

18 

Height2 

2  A 

Depth 

2B 

Height 

3  A 

Depth 

3b 

Height 

197  7 

Aug  10 

1 .  22m 

_ 

8.  5  lm 

_ 

1.  9  5m 

Aug 

16 

1 .  2  2  in 

- 

8.  20m 

— 

1 .  88m 

- 

A  ug 

2  3 

1.  2  5m 

- 

8.  82m 

- 

1.  9  3m 

- 

Aug 

30 

i .  35m 

— 

9,  COm 

— 

2.  10m 

— 

Sep 

6 

1 .  3  7  m 

— 

9.  04m 

— 

2.  05m 

— 

Sep 

14 

1.  32m 

- 

8.  63m 

— 

1.  89m 

- 

Sep 

20 

1.  35m 

- 

8.  95m 

— 

1 .  8  5m 

- 

Sep 

2  7 

1 .  3  3m 

— 

8.  9 2m 

— 

1 .  84  n 

— 

Oct 

4 

1.31m 

— 

8.  87m 

— 

1.  63  m 

— 

Oct 

18 

1 .  3  1m 

3.  5  m 

8.  7 Sm 

- 

1 . 64  c. 

1  1.3m 

Get 

26 

1.  33m 

2.  lm 

8.  76m 

7.1m 

1.  85m 

1  1.3m 

Nov 

1 

1.  33m 

2.  5m 

8.  79m 

7.  lm 

2.  0  5m 

1  1.3m 

Nov 

8 

1 .  32m 

2.  5m 

8.  78m 

7.  lm 

2.  06  m 

1  1.  3m 

Nov 

1  5 

1 ,  37m 

2.  7m 

8.  6  3m 

7.1m 

2.  02m 

10.  6  m 

Nov 

22 

1.61m 

2.  6m 

8.  8  7m 

7.1m 

2 .  0  7m 

10,  6  m 

Nov 

30 

1.  73m 

4.  2m 

8.  7  lm 

8.3m 

2.  06  ir. 

11.  3m 

Dec 

7 

1.  62m 

— 

8.  70m 

7.  lm 

2 .  lorn 

F 

Dec 

2  1 

1  .  S9m 

2.  8m 

8.  79m 

7.6m 

2.15m 

F 

1978 
Jan  3 

2.  0  1m 

3.  5m 

8.6  7m 

2.  20m: 

F 

Jan 

10 

2.  09m 

2.  8m 

8.  6  6m 

6.7m 

2.  24m 

F 

Jan 

25 

2.  C8m 

2.  8m 

8.  8  7m 

7.  1m 

2.  23  m 

F 

Feb 

13 

2.  19a 

F 

8 .  £  7  m 

7.1m 

2.  22m 

V 

-W. 

Mar 

8 

2.  22m 

F 

8.  90m 

7.  lm 

2.  2  lm 

f 

Apr 

26 

2.02m 

F 

8.8  7m 

7.  lm 

2.16m 

F 

Jun 

2  7 

1.  5  7m 

3.  5m 

8.8  7m 

7.  In: 

2.  0  7  m 

- 

J  ul 

25 

1.65m 

3.  5m 

8.  66m 

7.  lm 

2.  10m 

— 

depth  belcw  qround  surface 
2:  height  above  qround  surface 
F:  water  in  piezometer  was  frozen 


. 

. 


Depth  (meters) 
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TEST  HOLE  PIEZOMETER  IB 

77-3 


PIEZOMETER  1A 


> 

:  ?' 

v  : 

^cy/zs/* 

A  ' 

■< 

6.1  - 

•V 

0.  \  ' 

L_ _ 

Li 

»l'l 

3*.l 

///^/// 

3  8  cm  PVC  pipe 


■  cement 
-|4  cm  hole 
•bentonite  plug 

20-40  froc  sand 


Figure  3H  Piezometer  details,  nest  #1 . 


Depth  (meters) 


180 


TEST  HOLE  PIEZOMETER  2B 

77-6 


PIEZOMETER  2A 


0 *a>. 


3.0- 


24.4- 


32.0- 


55.0- 

574- 


66.8- 


XGi 


1 

iiiilc 


82.3- 


yord  hydrant 


fill 

cement 
Bizanni  plug 


14  0  cm  steel 
cement 


175  cm  hole 


13  0  cm  hole 


114  cm  steel 


torch  slotted 


casing 


casing 


6.1- 

9.1- 


r  - ; 


■  3  8  cm  P VC  pipe 


•  cement 
-bentonite  plug 

20-40  free  sand 
14  cm  hole 


Figure  32.  Piezometer  details,  nest  §2 


1  8 1 


TEST  HOLE  PIEZOMETER  3B 

77-2 


PIEZOMETER  3A 


Figure  33.  Piezometer  details,  nest  #3- 
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Aquifer  lest:  Piezometer  IB 
LSD  4-  I  5-6 7-24-W4M 


Date  test  started:  July  29,  1977 

Pre-test  water  level:  flowing;  static  is  1.67nt  above  ground 

Test-hole  depth:  29.0m 

Average  test  rate:  14  L/min 

Depth  to  top  of  aquifer:  17.4m 

Screen  set  at:  20.4m  to  22.9m 

Pumping  interval:  51  min 

Datum  point:  top  of  casing;  (40.43m) 


pumping  interval 

Time  (t)  since  pumping  Depth  to  water  Drawdown  (s) 

started 


1 

min 

0.9  4m 

2.  6  1m 

2 

min 

1.52m 

3.  19m 

3 

min 

2.  06 mi 

3.  73m 

4 

min 

2 .  6  2m 

4 . 2  Sm 

5 

min 

3.19m 

4.86m 

6 

min 

3,44  m 

5.  1  1m 

7 

mi  n 

3. 6bm 

5.33m 

8 

min 

4.40m 

6.  0  7m 

9 

min 

5.11m 

6.78m 

1  0 

min 

5.36m 

7.05m 

15 

min 

8.  09  m 

9.7  6m 

20 

min 

10.37m 

12.04m 

2  5 

mm 

13.  08m 

14. 75m 

30 

min 

15.  75m 

17.42m 

40 

min 

20. 93m 

22. 60m 

5  1 

min 

- 

- 

1 


• 

recovery  interval 
piezometer  IE 


lime  (t) 
pumping  . 

since 
s  top  fed 

Depth  to 
#a  ter 

Residua 1 
ardtado¥n  (s’) 

Ratio  t/t 

1 

min 

15.51m 

21. 18m 

52.  00 

2 

min 

18.  39m 

2  0.  06m 

26 .  50 

3 

do  in 

1 7. 50m 

19.  1  7  m 

18.00 

4 

min 

16. 74m 

18.4  1m 

13.75 

5 

min 

16.  ISm 

1 7.86m 

11.20 

6 

min 

1  5.  6  7m 

17.34m 

Q  ^  1 

7 

min 

15.  19m 

16. 86m 

8.29 

8 

mm 

14.63m 

16.30m 

7.  38 

9 

min 

14.  10m 

15. 77m 

6.  68 

10 

min 

1 3. 56m 

15. 25m 

6.10 

15 

min 

1  1.76m 

13.  4 3n: 

4.40 

20 

min 

9.91m 

1  1 . 5  8  m 

3.  55 

25 

min 

- 

- 

— 

30 

mm 

7.  15m 

8.82m 

2.70 

40 

min 

4.  3  3m 

6.00m 

2.28 

50 

min 

3.  05m 

4.72m 

2.20 

6  0 

min 

1.  7  8  in 

3.  4  5m 

1.85 

75 

min 

0.  3  7m 

2.  0 4m 

1.68 

■ 

BAPTISTE  LAKE  OBSERVATION  WELL  IB 


i. _ I _ I _ I _ I _ 1 _ I _ I _ I _ L 

U|  (S)  U«Op«OJQ  |0npi*8y  puo  (t)  U«Op»OJQ 


Figure  3^,  Aquifer  test,  piezometer  IB. 
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Aquifer  lest:  Piezometer  2B 
LSD  1-2  1-66-24- *4M 


Date  test  started:  July  28,  1977 

Pre-test  water  level:  flowing;  static  is  6.8m 

Test-hole  depth:  82.3m 

Average  test  rate:  35  L/min 

Depth  to  top  ot  aquifer:  55.  £m 

Screen  set  at:  55.0m  to  51.4m 

Pumping  interval:  20  min 

Datum  point:  top  of  casing;  {  +  0.3m) 


pumping  interval 
test  #  1 

Time  (t)  since  pumping  Depth  to  water 

started 


1 

mi  n 

2.  9 2m 

2 

min 

4. 5  7m, 

3 

min 

6 .  45m 

4 

min 

8.  38m 

5 

min 

10.  08m: 

6 

min 

11.68m 

7 

mi  n 

1  3. 26m 

8 

min 

14. 68m 

9 

min 

16 . 08m 

10 

min 

1 7. 22m 

1  5 

min 

22.  73m 

20 

mi  n 

- 

185 


above  ground 


Drawdown  (s) 

7.  72m 

1  1.  3  7m 

13.  85m 

15.  18m 

16. 88c 
1 8. 48m 
20. 06m 

2  1.46© 

22. 88c 
2  4.  Q2m 
29. 53m 


. 
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recovery 

interval 

piezometer 

1 B ,  test  #1 

Time  (t) 

since 

Depth  to 

Residua  1 

h  a  t  i  o  t  / 1  * 

pu  raping 

stop  fed 

water 

drawdown  (s’) 

1 

min 

2  1,  56m 

28 . 3ora 

2  1.00 

2 

min 

2  0.  2  7m 

2  7. 0  7m 

11.00 

3 

min 

19.  2 5m 

26 . 05m 

7.6  6 

4 

min 

18.  16m 

24.96m 

6.00 

5 

min 

1 6. 94m 

23. 74m 

5.00 

6 

min 

1 5.  77m 

22.57m 

4.  33 

7 

rain 

1 5. 04m 

2  1.84m 

3.86 

8 

min 

1  4.  1  5m 

20.  9  5m 

3.50 

9 

min 

13. 2  1m 

20.01m 

3.  22 

1  0 

min 

12.  32m 

19.  !2m 

3.00 

1  5 

min 

8.  0  6  ra 

1 4, 86m 

2.33 

20 

rain 

4.  3  2ra 

11. i 2m 

2.00 

25 

min 

1.  98m 

8.78m 

1. 80 

30 

min 

0.  OOra 

6.80m 

1.6  7 

187 


Figure  35.  Aquifer  test  no.  1,  piezometer  2B. 


Aquifer  lest:  Piezometer  2B 
LSD  1-2  1-66-24- 


188 


Date  test  started:  July  28,  1977 

Pre-test  water  level:  flowing;  static  is  6.8m  above  qround 

Test- hole  depth:  82.3m 

Average  test  rate:  18  1/ min 

Depth  to  top  of  aquifer:  55.8m 

Screen  set  at:  55.0m  to  5  1.4m 

Pumping  interval:  90  min 

Datum  point:  top  of  casing;  (  +  0.3m) 


pumping  interval 
piezometer  2B,  test  #2 


since  pumping 
started 

Depth  to  water 

Drawdown  (s) 

1  mm 

i.  52m 

8.  32m 

2  min 

2.31m 

9.  1  lm 

3  min 

3.  3  9m 

10.  19m 

4  mi  d 

4.  3  4m 

1  1.  14m 

5  min 

5.  23m 

12. 03m 

6  min 

6.  1  0m 

12.  90m 

7  min 

6.7  8m 

13. 58m 

8  min 

7.  6  2m 

14.  42m 

9  min 

8 .  3  6  n; 

15.  16m 

10  min 

8.  9  9m 

1 5. 79m 

15  min 

1  1.  92m 

1  8.  72m. 

20  min 

13.99m 

20.  79m 

25  min 

15.  7  9m 

22.  59m 

3  0  min 

1  7.  6  0m 

24. 40m 

40  min 

19. 20m 

26. COm 

50  min 

20. 22m 

2  7.  02m 

75  mm 

19. 18m 

2  5.  9  6q; 

90  mi  n 

24. 08m 

3  0.  6  8m. 

■ 
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recovery  interval 
piezometer  2B,  test  #2 


time  (t)  since 
pumping  stopped 

Depth  to 
water 

Residual 
drawdown  (s') 

Ratio  t/t ' 

1 

mi  n 

2  1.  72m 

28.  5 2 mi 

9  1.00 

2 

min 

20.  3  7m 

2  7,  17m 

46.  00 

3 

min 

19.  20m 

26. 00m 

3  1.00 

4 

min 

17. 98m 

24. 78m 

23.50 

5 

min 

16. 89m 

23.  69m 

19.00 

6 

min 

15.  83m 

22 .  63  m 

16. 00 

7 

min 

1 4. 83m 

21.  63m 

13.86 

8 

min 

- 

- 

- 

9 

min 

1 2 .  9  6  ir, 

19, 76m 

11.00 

10 

min 

12. 09m 

1 8, 89m 

10.  00 

15 

min 

7.  70m 

1  4.  50m 

7.00 

20 

min 

4 . 60m 

1  i.  60m 

5.50 

25 

min 

2.32m 

9.  12m 

4.6  0 

30 

min 

0.  02m 

6.82  m 

4. 00 

40 

min 

0 .  00m 

6,  80ro 

3.2  5 
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Figure  30.  Aquifer  test  no.  2,  piezometer  2B. 
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-HOLES 


. 


Lithology  Logs  of  the  Shallow  Test-Holes 


SHE  #  i 


SE  of  LSD  1-26-66-24  - k 4  ?i 


Depth 

(E) 

Sedi ment 

Description 

0.0- 

1-5 

silt 

sandy,  storey,  light  brown. 

1.5- 

3.  7 

till 

clayey,  dark  brown. 

3.7- 

4.6 

till 

clayey,  moist,  aark  grey-brown. 

4.6  - 

6 .  1 

till 

as  above  nut  with  greater 
proportions  of  clay. 

6.  1  - 

7.  6 

till 

sandy,  light  grey,  water  at  6.4m 

SITE  #2 

S  W 

of  LSD  3-26-6 6- 24-W4J1 

Depth 

(m) 

Sediment 

Description 

0.0  - 

3.  7 

clay 

moist,  mottled,  light  yellow  to 
medium  brown. 

3.  1  - 

5.2 

till 

very  clayey,  some  stones,  moist, 
dark  grey. 

5.2  - 

6,  1 

ti  11 

sandy,  light  grey,  water  at  5. 2m 

SITE  #3 

3  Vi 

of  LSD  1  6-2 7-66- 24- W4M 

Depth 

(m) 

Sedi men t 

Description 

0.0- 

0.  9 

sand 

clayey,  rusty-brown. 

0.9- 

2.  1 

clay 

moist,  grey- brown. 

2.  I  - 

4.  6 

sand 

medium  grained,  well  sorted. 

oxidized  where  saturated, 
r usty- brown,  water  at  2.6m. 


SITE  #4 

SW  of  LSD  6-22-66-24-W4M 


Depth 

(m) 

Sedi men t 

Description 

0.0- 

0.  6 

till 

very  sandy,  medium  grained, 
orange- brown . 

0.  6  - 

1.  5 

till 

mottled,  yellow  and  brown  bands. 

1.5  - 

3.0 

till 

clayey,  grey-brown,  water  at  1.5 

3.0- 

4. 6 

till 

same  as  above  but  dark-grey. 

' 


■ 
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SITE  *5 


SE  of  LSD  5-22-66-2 4- Vi 4 M 


Depth 

Cm) 

Sedi  merit 

Description 

0.0- 

0.9 

till 

very  clayey,  moist. 

dark  grey 

-brown 

0.9- 

7.6 

till 

very  clayey,  moist, 
at  3.7m,  dark  grey. 

wet  sand 
water  at 

layer 

6 .  4  a,. 

SITE  #6 


NE  of  LSD  7-2  I-66-24-W4H 


Depth  (m) 

Sedi men  t 

Description 

0.0-  1.5 

till 

silty,  moist,  grey-brown, 
water  at  0.9m. 

1.5-  3.0 

till 

sandy,  grey. 

SITE  #7 

S  W  of 

LSD  10-2  1-66-24-W4M 

Depth  {m) 

Sediment 

Desc  rip t ion 

0.0-  1.5 

till 

very  sandy,  layers  of  rusty -brown, 
grey,  purple  colours,  water  at  0.6m 

1.5-  2.  4 

sand 

fine  to  medium  grained,  brown. 

2.4-  3.0 

sand 

fine  to  medium  grained,  dark  grey. 

SITE  #8 

S W  of 

LSD  12-2  1-66-24-K4M 

Depth  (m) 

Sedi men  t 

Description 

0,0-  2.1 

sand 

medium  to  coarse  grained,  very 
stoney,  moist,  rusty-orange. 

2.  1  -  3.  4 

sand 

fine  tc  medium  grained,  moist, 
dull  brown. 

3.4-  6.1 

sand 

medium  grained,  well  sorted,  moist, 
rusty-orange,  water  at  4.6m. 

SITE  #9 


SE  of  LSD  2-28-66-24-K4E 


Depth 

(m) 

Sed i ment 

Description 

0.0- 

2.  7 

silt 

very 

some 

mottled,  clay-loam  texture, 
oxidized  areas,  dark  grey. 

2.7- 

4.6 

sand 

fine 

water 

to  medium  grained, 
at  2.7m. 

buff- yellow. 

■ 

SITE  #10 


SV  of  LSD  4-28-66-24-K4M 


Dept  h 

{&) 

Sediment 

Description 

0.  3  - 

2.  1 

sand 

medium  grained,  well  sorted, 
orange-yellow. 

2.  1  - 

7.6 

sand 

medium  grained,  well  sorted, 
buff-yellow,  water  at  6.7m. 

7,6  - 

9.  1 

sand 

medium  qrained,  sliqhtly  clayey, 
buff-yellow . 

SITE  #  1  ! 

SW 

of  LSD  4-33-6  6-24— k’  4  K 

Depth 

{m) 

Sediment 

Description 

3.0  - 

7.  6 

till 

very  clayey  and  sandy,  moist, 
qr ey-brown. 

7.6  - 

12.2 

till 

very  clayey,  dark  qrey,  water  at 

1  1.6m. 

12.2- 

13.  7 

till 

very  sandy,  light  qrey. 

SITE  #12 


N  W 

of  LSD  S - 2 8- 6 6- 2 4- w 4 H 

Depth 

(m) 

Sedi men  t 

Descri p tion 

0.0- 

1 .  8 

sand 

medium  qrained,  well 
rusty-orange. 

sorted,  dry 

1,8- 

3.  0 

sand 

medium  qrained,  well 
rusty-brown,  water  at 

sorted, 

1.8m. 

SITE  #13 

S  w  of  LSD  d-33-66-2i--«4M 


Depth 

(m) 

Sediment 

Description 

o.o- 

1.5 

silt 

oxidized,  sandy,  yellow-brown. 

1.5- 

4 . 0 

clay 

grey-brown. 

4.0- 

12.  2 

till 

very  clayey,  dark  grey,  water  at 

11.0m. 

' 
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SITE  #14 


N£  of 

LSD  1-  3-67-24-H4M 

Dept  h 

(E) 

Sedi me  c  t 

Description 

0.0  — 

7.  0 

till 

very  clayey,  sandy  near  the  surface, 
moist,  dark  brown. 

-j 

. 

o 

i 

1  0.  7 

till 

very  clayey,  dark  grey,  water  at 
9.5m. 

SITE  #15 

SW  of  LSD  16-  3-67-24-W4M 


Depth 

<rc) 

Sedi me  n  t 

Description 

0.0- 

0.6 

sa  nd 

light  brown. 

0.6  - 

2.  4 

till 

sandy,  wet,  white. 

2.4  - 

6 . 4 

till 

very  clayey,  oxidized,  very  wet, 
dark  brown. 

6.4  - 

9.  1 

till 

very  clayey,  dark  grey,  water  at 
7.  9m. 

SITE  #16 

NE  of 

LSD  1 5 -  1 0- 6  7- 2 4 - fc 4 M 

Depth  (  m)  Sediment 

Description 

0.0-  1.5  sand 

medium  grained,  well  sorted,  dry 

buf f-yel low . 

1.5-  4.6  sand 

medium  grained,  well  sorted,  wet 

sliqhtly  clayey,  buff-yeilow , 
water  at  3.0m. 


SITE  # 1 7 


NE  of  LSD  8-  15-6  7-24- W4M 


Depth 

(m) 

Sediment 

Description 

0.0- 

1.  5 

till 

sandy,  dry,  dark  brown. 

1.5- 

3.0 

till 

mottled,  oxidized  areas,  wet. 

dark  brown. 

3.0- 

4.  6 

till 

very  clayey,  wet,  darl  grey-brown. 

4.6  - 

12.  2 

till 

very  clayey,  wet,  dark  grey,  water 

at  1  I.  Cm. 


■ 


. 
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Depth 

0.0- 

1.5- 
6.1  - 
7.6  - 


Depth 

0.0- 

1.5- 
3.0- 

4.6- 


Depth 

0.0- 

2.7- 
4.6  - 


Depth 

0.0- 

1.8- 

4.6  - 
6.  I  - 
7.3  - 

7.6  - 


SITE  #18 

NR  of  LSD  10-11-67-24-W4H 


(IE) 

1 .  5 

Sediment 

till 

very 

Desc 

c layey , 

ription 
sana  y 

at  the  surface. 

6  .  1 

till 

grey 

very 

-brown, 
c  layey , 

moist. 

grey-brown. 

7.  6 

till 

very 

clayey. 

moist. 

dark  grey. 

10.  7 

till 

very 

clayey. 

moist. 

slightly  more 

slity  than  above,  dark  grey,  water 

at  10.1m. 


SITE  #15 

SR  of  LSD  4— 13-67-24- R 4 M 


(®) 

Sed i me nr 

Description 

1 . 5 

till 

sandy,  moist,  yellow-brown. 

3.  0 

till 

sandy,  wet,  yellow-brown. 

4.  6 

till 

sandy,  moist,  yellow-brown. 

6.  1 

till 

very  clayey,  dark  qrey-brcwn, 
water  at  4.6m. 

SITE  #20 

NR 

of  LSD  4-12-6 7- 24-14 M 

(m) 

Sediment 

Descr lption 

2.  7 

till 

clayey,  dry,  oxidized, calcareous 
deposits,  yellow-brown. 

4  .  6 

till 

very  clayey,  moist,  dark  grey-brown. 

15.2 

till 

very  clayey,  moist,  dark  grey-brown, 
water  at  14. 3m. 

SITE  #21 

NE 

of  LSD  I6-34-66-24-R4M 

(m) 

Sedi men t 

Descrip tion 

1.  8 

sand 

medium  gramea,  well  sorted, 
yellow-orange. 

4.6 

till 

very  clayey,  moist,  yellow-brown. 

6.  1 

till 

very  clayey,  moist,  grey-brown. 

7.  3 

till 

very  clayey,  moist,  dark  grey. 

7.6 

till 

very  clayey,  sandy,  dark  grey, 
water  at  7.3m. 

9.  1 

till 

sandy,  light  grey. 

SITE  #22 


NW  of  LSD  2-  2-66-24-W4M 


Depth  (m) 

Sediment 

Description 

0.0- 

1.  5 

sand 

medi urn 
water 

grained,  clayey,  light  grey 
at  0.  9m. 

1.5- 

3.  0 

sand 

clayey 

,  light  grey-brown. 

SITE  #23 


S 8  of  LSD  1  1-35-6 6- 24-8  4 M 


Depth 

(m) 

Sediment 

Description 

0.0- 

2.  1 

silt 

clayey,  dry,  orange-brown. 

2.  1  - 

5.2 

till 

mottled,  moist,  stoney,  dark  crown 

5.2- 

6.  1 

till 

very  clayey,  moist,  stoney,  dark 
grey . 

6.1- 

9.  8 

till 

very  clayey,  mcist,  dark  grey. 

9.8  - 

12.  2 

till 

sandy,  moist,  grey-brown,  water  at 
10. 7m. 

SITE  #24 

NE  Ot 

LSD 

16-1 7-66-24-W4E 

Depth 

(m) 

Sediment 

Description 

0.0- 

6.  1 

till 

very 

clayey,  moist,  dark 

brow  n. 

6.  1  - 

15.  2 

till 

very 

clayey,  moist,  dark 

grey. 

*  eo  water  * 


SITE  #25 

NW  or  LSD  I  4-28-66- 2 4- W 4-?! 


Depth 

(ra) 

Sedi meat 

Desc  ri p tion 

0.0- 

3.0 

till 

some  stones,  alternating  moist  and 
dry  layers,  dark  brown. 

3.0- 

6.  7 

till 

as  above  but  no  stones. 

6.7- 

9.  1 

till 

dark  grey-brown,  water  at  8.2m. 

■ 

. 
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Lithology  Logs  of  the  Deep  Test-Holes 


TEST-HOLE  76-1 
LSD  4- 15-67-24- W4M 

Depth  (m)  Sediment  Description 


0. 

0  -  7.6 

silt 

sandy. 

few  pebbles,  dark  qreen. 

7. 

6  -  7.7 

gravel 

7. 

7  -  2  1.3 

clay 

sandy. 

few  pebbles,  dark  qreen. 

2  1. 

3  -  24.4 

till 

silty 

and  sandy,  dark  green. 

24. 

4  -  2  5.9 

till 

silty 

and  sandy,  taking  water. 

dark  green. 

25. 

9  -  29.8 

sand 

29. 

8  -  38.  1 

till 

sandy. 

dark  green. 

38. 

1  -  39.6 

sand 

39. 

6  -  42.1 

till 

sandy. 

with  pebbles,  dark  green. 

42. 

1  -  43.3 

gr  avel 

composed  or  shale  and  ironstone 
fragments,  some  clay,  quartz  pebbles 

fairly 

rounded,  5  ce  in  diameter. 

43. 

0  -  54.9 

shale 

dark  grey-green. 

Deptn  (m) 

TEST-HOLE  76-2 

LSD  1-23-66-24- K4M 

Sediment  Description 

3.3 

- 

29.  9 

sand 

oxidized  on  top,  few  pebbles,  coal 

29.9 

— 

44.  2 

till 

fragments,  green-grey. 

sandy,  coal  fragments,  very  clayey  m 

44 . 2 

5  7.9 

till 

top  1. 5®,  pebbles  and  thin  gravel 

layers,  green-grey. 

as  above  but  with  slightly  more 

57.9 

_ 

6  1.3 

till 

gravel . 

sandy,  contains  hard  shale  fragments 

6  1.3 

64.3 

till 

and  ironstone  fragments,  possible 
flowing  zone,  reddish-brown, 
as  above  but  harder  drilling. 

64.  3 

- 

73.8 

till 

contains  shale  and  ironstone 

73.8 

79.  6 

san  d 

fragments,  r eddish-brown. 

very  silty,  some  coal,  taking  water. 

79.6 

83.  5 

sand 

80. 5 

- 

85.3 

sand 

very  silty,  some  coal,  few  bedrock 

85.3 

9  1.4 

sand 

fragments,  taking  water. 

very  silty,  contains  coal  chips. 

9  1.4 

93.9 

till 

taking  water, 
sandy. 

93.  9 

- 

97.  8 

till 

very  clayey,  contains  some  fine 

97.  8 

98.  1 

till 

gravel. 

as  above  but  with  coal  chips. 

98.  1 

- 

102.  1 

shale 

dark  grey-green. 
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12.  2 

36.  6 
40.2 
4  1.1 


IES'1-HOLE  76-3 

LSD  1  3-  35-66-24-K4M 

(ffl) 

Sediment 

Desc  rpi tion 

12.  2 

silt 

oxidized,  very  sandy. 

guite  a 

pebbles,  medium  grained  sand. 

3  6.6 

till 

sand  y. 

40,  2 

sand 

silty. 

4  1.1 

till 

very  clayey,  contains 

peobles. 

56.  7 

shale 

dark  grey-green. 

Depth  {id) 
0.0-  5.5 

5.5  -  20.  1 

20.  1  -  36.6 


1EST-HOLE  76-4 

LSI:  12-26-66-24-fc4H 

Sediment  Description 

sand  medium  grained,  brown, 

till  sandy,  dark  qreen. 

shale  dack  grey-green. 


Depth 

(m) 

Sediment 

0.0- 

54.9 

till 

54.  9  - 

57.9 

till 

57. 9  - 

62.  5 

till 

62.5  - 

7  3.  2 

shale 

1E5T- HOLE  76-5 
LSI  1 2-36-66-24-K4M 

Descri ption 

sandy,  contains  pennies,  dark  qreen, 
as  above  Dut  with  coal  chips, 
sandy,  contains  pebbles, 
dark  grey-green. 


Depth  (m) 
0.0-  20. 4 
2  0.  4  -  3  0.5 


IESI-HOLE  76-8 

LSD  13-33-66-24-V<4E 

Sediment  Description 

till  sandy,  few  pebbles,  green-grey, 
shale  dark  grey-green. 


Depth  (m) 
0.0-  7.6 

7.6-  42.  7 


LEST- HOLE  76- 1 1 

LSD  14-  12  — 67-24-W4M 

Sediment  Description 

till  sandy,  with  pebbles, 
shale  dark  grey-green. 
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I  ESI- HOLE  77-  I 
LSD  6-2  1-66-24- W4K 


Depth 

(ID) 

Sediment 

Description 

o.o- 

1  .  5 

sand 

fine  grained,  some  coal  cnips, 
reddish-brown. 

1.5- 

5.8 

sand 

fine  grained,  some  coal  chips,  grey. 

5.8  - 

6.  7 

sand 

medium  grained,  black. 

6.  7  - 

9.5 

sand 

medium  grained,  coal  chaps,  grey. 

9.5  - 

1  1.  9 

clay 

contains  alot  of  organic  material, 
dark  grey. 

11.9- 

12.  5 

sand 

medium  grained,  lost  circulation, 
dark  grey. 

12.  5  - 

17.  1 

c  lay 

no  pebbles,  grey. 

17.1  - 

22.  8 

shale 

soft,  grey-green. 

1 ESI -HOLE  77-2 


LSD  6-2  1-66-24- W48 


Cep  t  h 

(m) 

Sediment 

Description 

0.0- 

2.  4 

s  and 

medium  to  coarse  grained,  oxidized. 

2.4- 

12.2 

sand 

fine  grained,  contains  wood,  shells 
and  coal  chips,  grey. 

12.2- 

13.  7 

clay 

dark  brown. 

13.7  - 

19.  2 

clay 

black. 

19.2  - 

4  1.1 

clay 

grey. 

4  1.1- 

4  4.  2 

till 

sandy. 

44.2- 

4  7.2 

gr  avel 

medium  grained,  salt  and  pepper 
texture,  taking  water. 

47.2- 

5  1.8 

shale 

dark  grey -green. 

lESl-HOLE  77-3 
LSD  4-  1  5-67-24-H4H 


Dept  h 

(m) 

Sediment 

Description 

0.0- 

2.  1 

s  ilt 

sandy,  no  pebbles,  grey. 

2.  1  - 

4.  3 

sand 

medium  grained,  grey. 

4.3  - 

1  7.4 

clay 

no  pebbles,  grey. 

17.4- 

26 . 8 

sand 

fine  to  medium  grained,  grey, 
taking  water. 

26.8- 

29.0 

till 

sandy,  grey. 

' 


. 


lEST  HOLE  77-4 


LSD 

i  7-10-6  7-24- W48 

Depth  (in) 

Sediment 

Description 

0.0-  3.  1 

clay 

oxidized,  some  pebbles,  brown 

3.1  -  11.3 

till 

very  clayey,  grey. 

11.3  -  13.1 

till 

oxidized,  grey. 

13.1  -  1  S. 8 

shale 

dark  grey-green. 

IESI-HOLE  77-5 

LSD  8- 10-67-24-K4M 

Depth 

(  ®) 

Sediment 

Description 

0.0- 

3.  4 

clay 

oxidized,  some  pebbles,  biown 

3.4  - 

7 . 9 

clay 

grey. 

7.9  - 

8.  2 

gravel 

very  silty. 

8.2  - 

14.  G 

till 

very  clayey,  grey. 

14.  0  - 

14.6 

gravel 

silty,  some  coal  fragments. 

14.6- 

27.4 

shale 

darx  grey-green. 

1ES1-HCLE  77-6 
LSC  1-20—66-24- W4M 


Depth 

(m) 

Sediment 

Descn  ption 

0.0- 

30.  5 

sand 

30.5- 

55.3 

till 

55.  8  - 

59.4 

gravel 

silty,  feldspar  and 
fragments. 

5  S .  4  - 

82.  3 

shale 

dark  grey -green. 

. 
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K  [JTTEB I  A  N  BLETHfiE  N  COLONY  KLLL 
SE-36-66-24-W  UM 


Depth 

(E) 

Sediment 

Description 

0.  0 

— 

3.  1 

clay 

oxidized. 

3.  1 

— 

7.  6 

clay 

grey . 

7.6 

— 

12.2 

till 

soft,  grey. 

12.  2 

— 

12.8 

clay 

oxidized. 

12.8 

- 

19.  8 

till 

grey. 

15.8 

- 

2  1.3 

till 

silty,  grey. 

2  1.3 

- 

3  6.6 

sil  t 

calcareous,  grey. 

3  6.6 

- 

38.7 

sand 

interbedded  with  silt 

38.  7 

- 

42.  7 

silt 

calcareous,  grey. 

42.  7 

- 

45.  1 

clay 

hard ,  non-ca lcareous , 

45.  1 

- 

5  0.0 

sil  t 

soft,  grey. 

50.  0 

— 

52.4 

clay 

52.  4 

- 

55.  2 

sand 

55. 2 

- 

6  7.1 

shale 

*  drilled  and  logged  ny  Elk  Point  Drilling  Ltd  on  €9/09/06. 

*  static  water  level  is  31.5  e  below  ground  surract. 


, 
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8 recipi tation ,  Athabasca  2  meteorlogy  station 


(all 

values 

in  c  m ) 

Month 

1970 

197  1 

1  972 

1973 

1974 

1  975 

1976 

1977 

J  an 

2.  4  1 

5.  8  9 

2.7  7 

0.  5  1 

5.  6  1 

2.64 

1.  42 

1.8  6 

Fee 

1.  7  3 

0.  9  1 

5.3  1 

2.  77 

3,  53 

1.  04 

1.85 

0,  89 

Mar 

5.  49 

2.29 

3.6  1 

0.36 

3.  87 

1.  73 

2.  0  1 

2.53 

Apr 

0.  23 

2.  49 

2. 3  1 

3.  12 

4.96 

2.4  1 

0.  4  1 

1 .  73 

May 

1.  79 

1.  68 

2.  72 

8.33 

2.  3  I 

2.  54 

3.  35 

1  1.25 

J  an 

15,32 

12.  75 

11.40 

9.40 

2.  90 

9.  78 

€.86 

6.  1  0 

J  ui 

12.  73 

10.  69 

5.  76 

5.  18 

1  7.42 

13.  3  1 

8.  79 

19.  34 

Aug 

4.  98 

6.  15 

3.58 

7.  09 

5.  06 

1  1  .  74 

9.  50 

5.  4  6 

Sep 

3.  53 

1.60 

3.30 

7.  26 

6.  10 

2.34 

2.72 

5.3  7 

Oct 

3.6  3 

1.  19 

1.19 

4.75 

1.  04 

3. 65 

3.  1  5 

1.5  0 

Nov 

1.  19 

2.90 

2,  03 

6.  12 

1.  06 

0.6  1 

0.  56 

0,  9  7 

Dec 

3.45 

3.48 

2.08 

4.  83 

1.9  1 

4.  34 

4.24 

1.73 

Temperature, 

A thabas 

ca  2 

Meteorology  S 

tation 

(all  values  in 

°C) 

Month 

1  9  70 

1  97  1 

1972 

1373 

1974 

1975 

1976 

1977 

J  an 

-  19.  1 

-22.  2 

-20.5 

12.9 

-21.1 

-  13.6 

-  12.  3 

-14.4 

Feb 

-7.  3 

-9.  6 

-  1  7.  7  - 

1  0.  9 

-  10.  7 

-  14.  7 

-9.  3 

-0.  9 

Mar 

-8.  1 

-7.  6 

-5.  3 

-2,  1 

-11,1 

-8.  8 

-5,0 

-3.  3 

Apr 

3.  4 

3.  4 

1.  1 

3.  1 

4.  0 

0.  7 

6.  5 

6.6 

May 

10.0 

12.  6 

1  1.  8 

11.6 

8.  1 

9,6 

1  1.  3 

1  1.  5 

Jun 

16,6 

14.  1 

14.  2 

13.  6 

14.8 

13,2 

12.  8 

14.  8 

Jul 

16.  1 

15.  7 

1  4.  0 

16.2 

14.  9 

1  7.  8 

16.  9 

14.  4 

Aug 

15.  1 

16.6 

1  6.  3 

14.6 

12.4 

12.5 

16.  7 

12.  4 

Sep 

9.  2 

8.8 

5.  5 

9.4 

8.  3 

10.9 

12.  4 

9.  5 

Oct 

2.  9 

3.  8 

).  8 

4.  1 

b ,  5 

3,  6 

3 .  5 

5.  2 

Nov 

-9.5 

-6,  4 

-4,9  - 

14.  1 

-3.  7 

-5.  2 

-2.  1 

-6.  9 

Dec 

-20.3 

-18.6 

-  1 7.  7  - 

1  4.  8 

-6.0 

-13.3 

-10.6 

—  17.8 

' 

. 
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Potential  Evaporation,  I  horn th waits ’ s  Method 

(all  values  in  cm) 


Mon  th 

19  7  0 

15  7  1 

J  an 

0.  00 

0.  00 

Feb 

0.  00 

0 .  0  0 

Mar 

0.  00 

0.0  0 

Apr 

2.  60 

2.6  0 

May 

8.  02 

9.  75 

Jun 

12.  6 8 

10.  97 

J  ul 

12.  65 

12.  30 

Aug 

10.  80 

1  1.  69 

Sep 

5.8  0 

5.  5  1 

Oct 

1.77 

2.  1  7 

Nov 

0.  00 

0.  00 

Dec 

0.  GO 

0.  00 

1872 

1873 

1974 

0.00 

0.  00 

0.  00 

0.00 

0. 00 

0.  00 

0.30 

0.  00 

0.  00 

1.03 

2.46 

3.25 

9.63 

9.25 

7.  0  1 

11.40 

10.  73 

1  1.  84 

11.53 

12.  77 

12.  13 

11.84 

10.  56 

9.  42 

3.93 

6.  00 

5.  60 

1 .32 

2.40 

3.  77 

0.00 

0.  00 

0.  00 

0.00 

0.  00 

0.  00 

1  9  75 

19  76 

197  7 

0.00 

0.  00 

0.00 

0. 00 

0.  0  0 

0.00 

u.  00 

0.  00 

0.00 

6.97 

4.  4  1 

4.  47 

7.  99 

8.  7  1 

9.  18 

10.  60 

9.  89 

11.59 

14.  0  1 

1  2.  6b 

11.56 

9.  25 

11.62 

9.  16 

6.  90 

7.  34 

6.  04 

2.  25 

1.93 

2.99 

0. 00 

0.00 

0.  00 

0.  00 

0.  00 

0.00 

• 

1 

